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PULSED ELECTROMAGNETIC GAS ACCELERATION: 
ACCELERATION PROCESSES IN QUASI-STEADY ARCS 

R. G. Jahn, W. Von Jaskowsky, and K. E. Clark 
Department of Aerospace and. Mechanical Sciences 
Princeton University, Princeton, New Jersey 

(NGL 31-001-005} 


A series of basic experiments has demonstrated that pulsed 
plasma accelerators can attain quasi-steady operation in a few 
tens of microseconds when provided with suitably protracted and 
synchronized inputs of driving current and mass flow. 1 ^ This 
quasi-steady mode is characterized by stationary, diffuse dis- 
charge current .patterns, constant terminal voltage, and constant 
exhaust velocity in a range encouraging propulsion application. 

Such devices are thus of interest both as pulsed plasma thrusters, 
and as laboratory simulators of the corresponding steady flow ac- 
celerators in power ranges where the latter cannot be conveniently 
studied in steady operation. 

For example. Fig. 1 shows a schematic diagram of a quasi- 
steady accelerator employing a coaxial electrode configuration 
like that of the conventional magnetoplasmadynamic arc. Mass flow 
is provided by a gasdynamic shock tube which can very rapidly pro- 
duce a high pressure reservoir behind suitable injector orifices. 
Arc current is derived from an array of capacitors and inductors 
which, depending on the particular configuration, provide rectan- 
gular pulses from 4000 amp x 600 jusec to 140,000 amp x 20 psec. “ 
The corresponding range of arc power — 100 kW to 100 MW — is com- 
mensurate with relatively large arc chamber dimensions, which in 
turn permits more detailed study of the discharge structure than 
would be possible in smaller scale facilities. 

A considerable variety of diagnostic studies have been 
performed in this facility and in other related apparatus involving 
Kerr-cell photography, electric and magnetic probing, velocity 
field mapping, spectroscopy, laser interferometry, and terminal 
measurements of current, voltage and mass flow. From these an in- 
structive-, yet incomplete, picture of the arc operation begins to 
emerge. Figure 2 shows a series of Kerr-cell photographs of arc 
luminosity at a given current for three different mass flows: one 

matched to the current value in terms of the classical MPD arc cri- 
teria, the others representing conditions of mass "starvation" and 
"overfeed, " respectively. The importance of matching the mass 
flow to the arc current is demonstrated here and in many other ex- 
periments, 2 most notably in a tendency for underfed arcs to erode 
electrode and insulator material. 

Figure 3a displays the distribution of current within a 
17,500 amp argon arc as determined by magnetic probing. Two par- 
ticularly significant features are the participation of the entire 
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cathode surface in the discharge, which should minimize the ero- 
sion of that surface, and the rather small protraction of the 
current pattern downstream of the anode orifice in contrast to 
the experience with external field MPD accelerators „9-12 Fig- 
ure 3b shows contours of floating potential for the same arc op- 
eration. Roughly, this pattern consists of four domains of sig- 
nificant potential drop: (a) a 10“20-volt anode fall, (b) a 20- 

volt cathode fall, (c) a drop roughly equal to the ionization po- 
tential across the injected gas jets, and (d) a 100-volt envelope 
around the cathode. Of these, it is the last which dominates the 
potential pattern, and probably holds the hey to the performance 
range and efficiency of this type of accelerator. 


Time-of-f light velocity measurements made downstream with 
biased double probes yield values inconsistent with any simple 
electrostatic acceleration model based on the potential pat- 
terns. p or the sample case of fig. 3, centerline ion speeds 
of the order of 25,000 m/sec are found, which imply ion energies 
of some 150 eV, and are therefore far above the difference between 
the anode and downstream plasma potentials, far above the Alfv&n 
critical value, ^ and, in fact, comparable to the full interelec- 
trode potential. Apparently some complex but effective electro- 
magnetic-electrothermal conversion process functions within or 
just downstream of the cathode potential envelope. Small wedges 
placed in this flow indicate that it is substantially supersonic, 
and that the ion temperature is of the order of 10 eV. 


While the full import of these various observations is not 
clear, the prognosis for efficient gas acceleration by this quasi- 
steady technique seems favorable. Repetitive pulse operation of 
thrusters of this type should combine the advantages of the high- 
power arc operation with tolerable mean power requirements and 
variable thrust capability via duty-cycle adjustment. ^ Systems 
aspects of this concept, including the major problems of power- 
conditioning and mass injection are also under basic study in 
this program, along with the development of suitable techniques 
for assessing the performance of thrusters of this type. 16 
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CURRENT AND POTENTIAL PATTERNS IN 
QUASI -STEADY MPD ARC 
J = 17.5 kA 
m = 5.9 g/sec 


FIGURE 3 


DYNAMIC CURRENT- VOLTAGE CHARACTERISTICS OF A PULSED 
MEGAWATT MPD-ARC 

by Charles J. Michels, NASA Lewis Research Center 

The investigation of the voltage-current (V-I) characteristics of a 
pulsed, megawatt, MPD-arc is described for various amplitudes of 
auxiliary magnetic field (0.5 to 2.0 T). The pulsed arc is powered 
by a crowbarred capacitor bank and the data is dynamically gathered 
for a hundred microseconds after the voltage and current traces have 
stabilized. The tungsten cathode is heated to emitting condition 
just prior to powering the arc. A puff of argon or nitrogen gas Is 
admitted into the evacuated arc chamber and allowed to distribute 
prior to powering the arc. From cold flow transient pressure traces, 
the time to power the arc is determined, as well as the chamber pres- 
sure at firing time. The experiments were operated for two different 
weight flows through the arc annulus, 6.0 g/sec , and 1.5 g/sec. 

The gross physical processes, such as the V-l characteristics, imped- 
ances, and plume characteristics described in this paper, are part of 
the means of identifying these arcs in relation to the varied forms 
of arcs in the literature. This paper extends the knowledge to the 
0,2-2. 5 megawatt level. Spoking instabilities are not noted for these 
power levels and time durations. 

The V-l characteristic obtained for each shot shows a positive-sloped, 
curvilinear trend. A series of these characteristics, for currents 
from 2 to 14 kA, shows them to be dynamic characteristics (as defined 
in Ref. 1) and not connected curves. A static characteristi c curve 
can be estimated from the family of dynamic characteristic curves. As 
the magnetic field is increased, the dynamic curves, although still 
curvilinear, tend to be more identifiable as part of a connected curve. 
For instance, at 20 kG, the dynamic curves are almost part of the esti- 
mated linear static character i'st i c curves. That is, the increased field 
has made the arc system almost res istor- 1 i ke and the short duration 
pulsed operation is sufficient to obtain static V-i characteristics. 
Positive sloped static character istics for higher arc current cases are 
observed as suggested in Ref. 1. 


REFERENCES 
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EXPLORATORY ELECTROMAGNETIC THRUSTER RESEARCH 


Adriano C. Ducati 

Plasmadyne, a division of Geotel, Inc. 
Santa Ana, California 

(NAS 1-9297) 


fl 2 31 

Recent work v ’ 5 ' has been concerned with obtaining an improved 
understanding of MPD thruster operation including pulsed (quasi- steady) (4, 5) 
as well as continuous arc devices. A large fraction of the effort has fo- 
cussed on evaluating the effect of ambient conditions on the measured per- 
formance. It has been demonstrated that ambient gas is entrained in a 
region of the plume that forms a part of the arc path and can therefore be 
subject to electromagnetic forces. 

Using a large insulating vacuum tank it was found to be possible to 
detect the flow of charged particles in the plume utilizing a large coil 
wrapped around the test chamber to introduce changes in the ambient mag- 
netic field strength. The field deflects portions of the plume which carry 
charged particles (as distinguished from excited particles). The phenome- 
non was observed in the insulated tank as well as a nonmagnetic metal tank, 
indicating that large scale current flow paths in the plume do not depend on 
a return flow path through the tank. Other observations (w show evidence 
of a flow of electrons returning in an annular path to the anode with suffi- 
cient intensity to cause severe heating of the exterior anode face. These 
results provide evidence that a significant fraction of the arc current fol- 
lows a path in the luminous plume outside of the thruster. 

To evaluate the rate of entrainment of ambient gas in the excited 
portion of the plume, tests were made with hydrogen as the working fluid 
and with an additional flow of hydrogen or some other- gas introduced di- 
rectly into the test chamber. It was found that the arc voltage is strongly 
dependent on both the ambient pressure and the type of ambient gas used 
(although a constant arc chamber pressure is maintained, and the propel- 
lant gas is not changed). Also, the color of the plume depends on the en- 
vironmental gas (apparently even at the center of the jet). 

These observations suggest that the jet mixes very rapidly with the 
surrounding gas, and strongly indicate that a large fraction of the arc en- 
ergy is expended on ambient gas mixed with nropellant in the plume. The 
rapid mixing is confirmed by another study' (°) in which detailed surveys 
were made to determine the flow pattern in the vicinity of a small jet in a 
vacuum chamber. This work also indicates that more rapid entrainment 
rates can be expected as pressure is reduced as long as continuum flow 
conditions predominate. A clear cut evaluation of entrainment effects ap- 
pears to require very precise measurement of thrust for a range of ambi- 
ent pressures in the molecular flow regime. For steady operation this is 
only possible with very small thrusters; and has not yet been accomplished, 
even though a sizable effort has been expended toward this end. It is con- 
cluded that we are not yet able to predict the performance of continuous 
(steady) MPD thrusters in space from laboratory, tests. 
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On the other hand, quasi-steady MPD thrusters appear to offer a 
better opportunity for obtaining performance data in ambient conditions ap- 
proaching those in space. The test chamber can be thoroughly evacuated 
before triggering a pulse. No interference will be felt until the pulse front 
reaches the chamber wall and is reflected back to the vicinity of the thruster. 
Interference free periods hundreds of microseconds long can be obtained 
in this manner in test chambers of moderate size. During this interval 
the 'thrust level can be high enough to produce an impulse that can be read- 
ily measured. 

Another important reason for investigating performance in the 
quasi-steady mode is that MPD thrusters show progressively better per- 
formance as the power level is increased. Since much higher power levels 
can be handled in momentary. pulses than would be feasible with steady 
operation, strong interest has developed in determining the performance 
attainable with quasi-steady MPD thrusters. However, performance has 
not yet been determined accurately enough to establish the degree of corre- 
lation with theoretical models or with steady state performance. 

Work has been initiated for determining the performance of quasi- 
steady arc jets operating at energies as high as 20,000 joules per pulse. 
These devices typically operate for a large fraction of a millisecond and at 
power levels measured in megawatts during a pulse. The investigation is 
concentrating on the use of solid and liquid propellants which are vaporized 
by the action of the arc. A critical design problem which is being actively 
studied is the choice of a feed system that will provide propellant at the 
proper rate to match the power pulse, and avoid wasteful expenditure of 
propellant. The problem is complex because it involves the power pulse 
shape and the electrode design as well as the propellant properties and the 
propellant feed system itself. 

Some refined instrumentation is required to accurately determine 
the performance of quasi-steady thrusters. For precisely measuring small 
values of thrust, a torsional ballistic pendulum is being developed. The 
complete pulse generating circuit is mounted on the thrust platform so that 
the connecting leads can be reduced to instrumentation and control lines 
and light power leads for recharging the condensers. Other equipment used 
includes a high speed motion picture camera to record the growth and tail- 
off of the plume, a spectroscope which is used to obtain data on the material 
in the plume by repeatedly exposing the same film during successive pulses 
until an adequate record is obtained, and various probes for determining 
conditions in the plume. 

An effort will be made to establish trade-off factors which will per- 
mit a reasonable choice to be made of design variables such as thruster 
geometry, size, duty cycle, and propellant; and tests will be made at a 
range of vacuum levels to see if performance is affected by ambient pres- 
sure to a measurable degree. 
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DESIGN OF A QUASI-STEADY MPD THRUSTER AND 
SPECTROSCOPIC VELOCITY MEASUREMENTS 

J. M. Hoell, J, Burlock, and 0. Jarrett 
NASA, Langley Research Center 
Hampton, Va. 


Performance data obtained for d.c. MPD thrusters during the past few 
years have indicated an efficiency increase with power. To verify and 
exploit this trend, the power applied to d.c. arcs has been extended to the 
limit of conventional electrode capability and, even more important, beyond 
the capacity of forseeable space-rated power supplies. The short pulse mode 
of operation (i.e. pulses of 1 to 10 ys duration) which has been studied by 
a number of laboratories solves the electrode problems and when repetitively 
pulsed permits use of average powers well below space-rated power supply 
capabilities. However, this mode of operation has had problems concerning 
proper matching of the propellant distribution with the short current pulse. 

An alternate mode of operation which takes advantage of the increasing 
efficiency with power associated with the true steady state MPD arc and 
with repetitive pulsing provides the low average input power condition of 
the short pulsed mode is a mode of operation, designated as "quasi-steady" 
(Q.S.). In this mode a high current pulse and a propellant pulse are applied 
for periods of time which are long enough to allow important parameters to 
reach steady state, i.e., v l msec. This note will describe research which has 

been underway at Langley during the past few months to develop a Q.S. MPD 

thruster and to spectroscopically measure the velocity distribution in a Q.S. 
thruster. 

The MPD arc used during these tests is a hollow cathode device which has 

been operated at Langley for a number of years. The cathode is a barium impreg- 

nated tungsten tube 2 inches long by one— quarter inch O.D. by one— eighth inch 
I.D. The anode is a copper- tungsten ring with a 1 inch I.D. and 32-1/16 inch 
O.D. gas-feed holes around the circumferance. A block diagram of the elec- 
trical circuit is shown in figure 1. The Q.S. power supply is coupled to the 
electrodes through an impedance matching transformer designed to work into a 
load of 0.14 ft. Operating into a properly matched load this power supply is 
capable of a square wave output of 1 megawatt for 1 msec with a rise and fall 
time of approximately 50 ys. Using the circuit shown in figure 1, the arc 
can be operated with or without a D.C. current and with or without an external 
B-field before applying a Q.S. current pulse. The time control block consists 
of several time-delay generators which are used to synchronize the Q.S. current 
pulse with the spectrograph shutter, the D.C. arc, and the gas pulse. 

In developing a Q.S. thruster, the first component considered was the 
Q.S. power supply. Of equal importance, however, is the gas-feed system and 
the region of gas injection. Using the present MPD arc, gas can be injected 
from three different positions: (1) through the hollow cathode, (2) through 
the anode, and (3) between the electrodes. Tests have been made using the 
first two methods. A fast-acting gas value capable of producing a gas pulse 
with a rise time of v 100 ys and a duration of several millisec. was used to 
introduce the gas into the arc chamber* The rise time of gas pulse for both 
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cathode and anode injection was approximately 500 ps. The light intensity 
produced by the 4806.017 £ line of argon was used as a measure of the 
effectiveness of gas injection from a particular position. It was found 
that when the propellant was injected through the anode, the light intensity 
from the 4806.017 S. line increased by approximately a factor of 5 over that 
for cathode injection. This suggests that for the latter case most of the 
gas was missing the discharge. A similar test is planned for injection 
between the electrodes. 

It Is of some concern in the development of a Q.S. thruster whether or 
not "quasi-steady" thermionic emission from the cathode exists during the 
current pulse as this relates to cathode life time, efficiency, and dupli- 
cation of the steady MPD mode. As a check on the emission from the cathode, 
experiments were made in which the pulsed current and voltage was measured 
for first a hot cathode and then a cold cathode. For the hot cathode case 
the Q.S. pulse was applied ^ 20 millisec. after a d.c. arc of 40 amps used 
to heat the cathode was extinguished. Figure 2 shows the current and voltage 
traces for these too conditions as well as for the condition of hot and cold 
cathode with B = 620 g. It can be seen that even for a pulse current of only 
2000 A there is no measurable difference between a hot and cold cathode for 
the • same external field; there is, however, considerable difference between 
the case of B ^ 0 and B = 0. 

A major part of the Q.S. program at Langley consists of spectroscopically 
measuring the exhaust velocity from a Q.S. thruster. For d.c. operation 
accurate measurements are relatively straightforward using a spectrograph with 
a resolution of 0.03 S or better. ^ For pulsed operation the major obstacle 
is insufficient light for good exposure for each pulse; consequently, for the 
pulsed powers that have been used to date, it has been necessary to use from 
three to six pulsed exposures per plate. The preliminary results that have 
been obtained have qualitatively demonstrated that just as the region of pro- 
pellant injection is important (see above) so is the uniformity of gas distri- 
bution. For example, figure 3 shows a radial distribution of the axial velocity 
which exhibits a radial gradient in axial velocity. The assymmetry in the 
axial velocity corresponds to an assymmetry that has been measured in the pulsed 
gas distribution. The unfavorable propellant distribution is due to the use of 
a single valve to inject gas from one side of the anode. To eliminate this 
nonuniformity in the velocity distribution a more uniform gas injection system 
is being designed; extensive surveys of velocity distribution will be made for 
this system. 
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Figure 1. - Schematic of Q. S. thruster system 
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a. Hot cathode, B = 0 k- Cold cathode, B - 0 


Upper trace : 100V/cm 
Lower trace : 2000 A/ cm 
Time sweep : 1 ms/ cm 



c. Hot cathode, B = 620 g 


d. Cold cathode, B = 620 g 


Figure 2. - Pulsed current and voltage 
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SPECTROSCOPIC STUDY OF ION-NEUTRAL 
COUPLING IN PLASMA ACCELERATION* 

A. Co Malliaris and D. R„ Libby 

Applied Technology Division, Aveo Corporation, 
Wilmington, Massachusetts 


The acceleration of a partially ionized gas by electromagnetic 
means might he quite inefficient, if the coupling between the ions 
and the . neutrals, for momentum transfer, is weak. Poor efficiency 
may arise from large velocity disparities between ions and neutrals 
in the same plasma flow. In the work reported here the absence or 
presence of a sufficiently strong coupling is generally revealed by 
velocity disparities determined spectroscopically, for several species, 
both ionic and neutral. Velocities are determined by Doppler shift 
measurements of selected spectral lines. Such lines are obtained not 
only from species which occur naturally in the main propellant, but also 
from any desired additives. 

Results regarding two different situations are being obtained, 
from experiments where two basically similar axi symmetric configurations 
operate under significantly different conditions. In the first case, a 
large disparity in ion-neutral velocities is observed, while in the 
second case, the said disparity disappears. In this second case, the 
common ion-neutral velocities are in good agreement with the center -of - 
mass velocity, as determined independently from measurements of motion- 
ally induced potentials in the flow, as well as from measurements, of 
the local momentum and mass fluxes. 

A satisfactory interpretation of the results is obtained when the 
calculated mean free path for momentum exchange is compared with a 
typical dimension of the electromagnetic acceleration region. 

A detailed paper will he presented at the AIAA 8th Aerospace 
Sciences Meeting, New York, N. Y. , January 1970. 


Work supported by NASA/LaRC Contract NAS 1-9298 
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NEUTRAL DENSITY MEASUREMENTS IN AN NH^ MPD ARC* 

Per Gloersen, General Electric-Space Sciences Lab, King of Prussia, Penna. 


The program described in this paper was to determine experimentally the 
range of applicability of a vacuum ultraviolet absorption spectrometry tech- 
nique to detecting neutral atom and molecule constituents in the exhaust stream 
of an MPD arc jet. Specifically, emphasis was on measuring the nitrogen atom 
densities in an ammonia-fueled MPD arc jet exhaust. To date, no other tech- 
niques have been established for measuring this quantity, possibly important for 
the corroboration of arc jet propulsion efficiency figures obtained from thrust 
and input mass flow measurements. An analysis is presented for the reduction 
of transmissivity data in terms of nitrogen atom flow rates. The main purpose 
of tnis program was to establish the experimental technique for this new appli- 
cation; more systematic use of this technique to help determine the operating 
characteristics of the NH 3 MPD arc would be undertaken in subsequent programs. 


The light source used in these studies consisted of an air-cooled quartz 
capillary through which nitrogen flowed. A 250 micron slit was placed at the 
vacuum chamber end of this tube to permit higher operating pressures in the 
tube than are maintained in the vacuum test chamber. Excitation is achieved in 
a tuned cavity driven by a 2.145 KMC, 125 watt oscillator. The background light 
source produced a 1200 S. nitrogen atom line received in the spectrometer, re- 
flected from a freshly coated mirror', some 15 times brighter than that received 
directly from the arc in ammonia; thus, meaningful absorption measurements were 
readily accomplished. 


The arc was constructed according to plans supplied by Langley Research 
Center and is identical with one of the configurations used in a number of 
experiments. Absorption measurements were made at a downstream position in 
the arc plume about 52 cm from the external nozzle (see Figure 1) . The flow 
range covered was from 6-18 mg/sec. The arc was operated at 300 amperes and 
63-80 volts dc, depending on the flow rate. The magnetic field was maintained 
at 5000 gauss for all of the runs. Observed transmissivity (I/I 0 ) of the arc 
under these conditions is indicated in Figure 2, along with similar data for 
the cold ammonia flow with the arc power turned off. 


Under the conditions of the present experiment, the background emission 
line is much narrower than the absorption line in the arc plume, principally 
because of the symmetrical Doppler broadening due to the arc plume spreading. 
Therefore, the only useful part of the absorption coefficient, k (\^) , is at the 
line center, which we shall call k Q . Thus, one can immediately write the ex- 
pression for the integrated absorption as 

I/I = e" k o nL ( 1 ) 

o 


= (1/k L) In (I /I) 

0 o 


Solving for the unknown. 


n 


( 2 ) 
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where, for a Doppler-broadened line, 

k = Tfr e cfA o /vtan<, 
o 


and 

r = 2.82 x ICT 13 cm, 
e 

10 

c = 3 x 10 cm/sec, 

f = 0.35, X = 1200 9., K. is the half angle of spreading in the exhaust 
plume and v is the average axial velocity of the absorbing species in the 
plume. The total rate of flow of the absorbing species is given by 

N = nAv = (7T*/16 r cf A ) Ltan°Cln (I /I) (3) 

e o o 

where the radius of the cross-sectional area A has been taken as L/4. K was taken 
to be 20°, as defined by the arc position and baffle aperture. It should be 
noted, however, that the visual appearance of the arc plume was indicative 
possibly of a somewhat smaller angle. Based on Kogelschatz' s results from his 
Doppler shift measurements in a similar arc," for the nitrogen atom is taken 
to be 3 x 105 cm/sec. over the operating range covered here. Two things 
should be noted in this regard. First, the range of operating conditions in 
these experiments extended to lower flow rates than those of Kogelschatz; it 
is assumed that the constancy of the velocity noted by him extends also to 
the lower flow rates. Second, Kogelschatz T s observations involve the -^ luminous 
nitrogen atoms whereas these experiments involve the ground state atoms. In 
view of the short radiative lifetime of the atomic states, it is assumed that 
the excited atoms observed by him resulted from re-excitation in the exhaust 
plume as a consequence of the extended currents from the arc. It therefore 
seems reasonable to assume that the excited and ground state atoms are moving 
with the same average velocity. 


Using these data and assumptions. Equation 3 can be evaluated for a number 
of different operating conditions. The results are shown in Figure 3. It can 
be seen that maximum fractional amount of ammonia observed in the plume was 
about 1.5% and that this fraction is almost independent of flow rate. The 
absolute accuracy of these results is no better than the + 50% available for 
the oscillator strength, f, and depends also heavily on the assumed appli- 
cability of the value used for v\ Thus, if the velocity of the non-luminous 
nitrogen atoms were significantly higher than observed by Kogelschatz for the 
luminous ones, the role of these consituents could be converted from minor to 
maj or . 


* 


Work performed under Contract No . NASI-8461 




TRANSMISSIVITY AT 1200 £ 



Figure 2* 




NITROGEN ATOMS IN THE PLUME 



Figure 3. Fraction of ammonia molecules at the inlet appearing as nitrogen atoms 
in the plume, 52 cm downstream from the nozzle. 
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SUMMARY OF LEWIS MPD ARC THRUSTER PROGRAMS 
by George R. Seikel, NASA Lewis Research Center 

This presentation will briefly review the various phases of Lewis 1 MPD 
arc thruster program including: (l) low power (100 to 1000 W) radia- 

tion cooled thruster; (2) 30 kW radiation cooled NH 3 thruster; (3) 
spectroscopic exhaust diagnostic measurements; and (4) MPD arc spoke 
phenomena. The scope of this program has been to investigate the 
feasibility of building lightweight completely radiation-cooled long- 
life thrusters and to simultaneously perform diagnostic studies of 
the plasma processes with the aim of improving thruster performance. 

A low power thruster has been operated in a hands off life test of 
638 hours, except for 11 intentional starts and stops. An automatic 
current control was used in this test, but with more recent cathode 
configurations no control appears necessary and no significant erosion 
has occurred during 100 hr. tests. The propellant yielding the highest 
performance is xenon but operation with many other propellants has been 
possible including a combination of H 2 and CO 2 . Recently a lightweight 
thruster was constructed. This thruster has very low power radiation- 
cooled magnetics and starter circuits. Presently parametric tests of 
this thruster are in progress. 

The present performance of the 30 kW NH 3 radiation cooled thruster is 
shown in an attached figure. This Lewis data shows both that the 
efficiency increases with arc power and that at constant power the 
specific impulse can be varied over a wide range. This thruster was 
developed and life tested for over 500 hrs. by McDonnell Douglas under 
contract to Lewis. It represents the first truly radiation cooled 
high power MPD arc thruster. The thruster is designed with adequate 
insulation between the anode and the magnet so that a radiation cooled 
permanent magnetic could be utilized. The anode heat loss (40 to 60% 
of the total power) is the major inefficiency in this thruster, and 
the requirement that this heat be radiated by the anode limits its 
power level to 40 kW. We are looking into the processes which affect 
the arc voltage particularly the way current attaches to the cathode, 
in the hope that by learning to control this we can raise the arc volt- 
age to lower the anode losses while maintaining the current distribu- 
tion that is presently yielding as thrust oyer 70 % of the net power put 
into the plasma. Why the conversion of the net plasma power into thrust 
can be so efficient has been shown by our spectroscopic measurements in 
the exhaust. These show both that the exhaust can be primarily high 
velocity neutrals and that the acceleration takes place in the magnetic 
nozzle outside the physical thruster. 

At present it is not clear whether the current spoke generally present 
in these thrusters is detrimental to performance. It does, however, 
greatly complicate theory and experiments to no clear advantage. There- 
fore, in addition to analytically and experimentally investigating spoke 
phenomena we are looking for ways to get rid of it. Preliminary tests 
with the low-power thruster for example indicate that use of inductors 
in series with an azimuthal ly segmented anode will eliminate the current 
fluctuations on the segments. 
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500 HOUR TEST OF A RADIATION COOLED MED ARC THRUSTER 


D. ¥. ESKER AMD J. C. KROUTIL 
MCDONNELL RESEARCH LABORATORIES 
ST. LOUIS, MISSOURI 63166 


The successful development of a MPD arc thruster to flight hardware 
status requires a demonstration of long term reliability. The work des- 
cribed represents an initial step towards this goal. Specifically, a 500 hour 
life test of a radiation cooled MED arc thruster is described. The thruster, 
(Fig. l), developed under contract to NASA, consisted of a radiation cooled 
anode- cathode- insulator assembly connected to a water cooled electromagnet. 

The design and performance characteristics of this thruster have been reported 
previously. 3 - The test was conducted in a three foot diameter vacuum chamber 
at a background pressure of approximately 9*0 x 10~3 porr with a propellant 
mass flow rate of ammonia of 30 mg/sec. Test instrumentation included arc 
voltage and current, magnet voltage and current, propellant flow rate, anode 
surface temperature, magnet heat transfer, and photographic observation of the 
interelectrode region. 

The test was initiated with a propellant flow rate of 30 mg/sec., arc 
current of 5^0 amperes, and axial magnetic field at the cathode tip of 0.14- 
tesla. These conditions produced a specific impulse of 1770 seconds with an 
efficiency of 17 . 1 percent. Rie thruster operated in a high voltage mode con- 
dition for the majority of the run with only several hours of low mode opera- 
tion at the start and end of the test. The arc voltage gradually increased 
from a value of near 55 volts early in the test to 64- volts after h-25 hours. 
Hiotographs of the cathode electrode indicated that the voltage increase 
correlated closely with the formation of a nodule or tit at the tip of the 
conical cathode. This nodule left the cathode at t» 4-28 hours with a subse- 
quent decrease in voltage, after which a new module began to form. 

The test was voluntarily shut down after 508. 5 hours for inspection of the 
thruster. An additional 4-5*5 hours was logged before the thruster was disassem- 
bled. The increasing yoltage characteristic served to increase the thrust so 
that the performance of the thruster varied over the course of the test, with a 
specific impulse of 1900 seconds being recorded prior to the shut down for in- 
spection. Erosion of the boron nitride insulator was the greatest at 4-1$ of 
its original weight. The cathode suffered a 2.9$ weight loss and the anode 
gained a slight amount of weight over its initial value. 

The results of this test suggest that the cathode electrode is the most 
significant component affecting long term operation rather than the anode. 

The control of voltage mode and the erosion of the cathode and insulator are 
items for future study. 


^Esker, D. W. , Kroutil, J. C. and Checkley, R. J. , "Design and Performance of 
Radiation Cooled Arc Thrusters", AIAA Paper 69 - 24-5, March 1969 . 
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FIGURE 1 MCDONNELL RADIATION COOLED MFD THRUSTER 
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EXPERIMENTS FOR PLASMA. DISTURBANCE IN MPD ARC 
AND LINEAR HALL ACCELERATOR 
IN THE PLASMA PHYSICS AND GAS LASER BRANCH 

0. Jarrett* F* Allario, B. D. Sidney and R. V. Hess 
NASA, Langley Research Center 
Hampton, Va. 


Onset of a rotating disturbance has been examined in two MPD accelera- 
tors, the linear Hall accelerator (LHA) and the MPD arc, as shown in figures 
1 and 2. This paper summarizes the results previously obtained-^-* 2 an( j 
presents results of recent experiments. 

The geometry of the MPD arc and the LHA is discussed in detail in 
references 1 and 2. The LHA has an approximately radial magnetic field 
and an axial electric field while the MPD arc has a coaxial electrode 
arrangement with a diverging axial magnetic field. The onset of the 
rotating disturbance in the two devices shows certain similarities despite 
the difference in power range; the LHA operating typically at a power of 
0.1 kw and the MPD arc as high as 30 kw. In figures 3 and 4, the critical 
magnetic field, B^, at onset of the instability is shown for the MPD arc 
as a function of mass flow and for the LHA as a function of pressure. Since 
pressure in the electrode region increases with mass flow in either device, 

B c increases with mass flow in each case. In addition, for a given mass 
flow, increases with molecular weight. 

Through variation of the radial magnetic field distribution near the 
cathode of the LHA (due to the magnetic field of the cathode heating current) 
it was determined that the total magnetic field near the cathode determined 
the onset of the instability in the entire electrode region. Because of the 
similarity in behavior of the rotating disturbance in the two accelerators, 
it was of interest to examine the MPD arc to see if the onset of the disturb- 
ance also occurred in the region near the cathode. 

In order to lower the power range of the MPD arc closer to the LHA, 
a hollow instead of a solid cathode was used. Stable operation of the MPD 
arc with a hollow cathode could be maintained from below 1 kw to 25 kw. 
Furthermore, a second solenoid with an independent power supply was in- 
stalled downstream from the anode to permit changes of the magnetic field 
configuration in the electrode region. The tip of the cathode was 
repositioned 1/4 inch behind the anode so that superposition of the mag- 
netic fields of the two coils permitted a wide variation of magnetic field 
in the electrode region. By comparing results from five different magnetic 
field configurations obtained at onset, it was determined that the axial 
magnetic field in a region close to the center plane of the anode determined 
the onset of the rotating instability. This result appears to be different 
from that obtained in the LHA where the magnetic field near the cathode 
determines the onset. The following considerations suggest, however, a 
common basis for the phenomena. 
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Previous studies of the MPD arc indicate a downstream bulging of 
the current, with formation of a cathode jet which is closer to cathode 
than anode potential. 3 It is reasonable to expect then that electrons 
emitted from the cathode follow the axial magnetic field lines for some 
distance downstream of the cathode, eventually crossing the magnetic field 
lines where the combination of crossed electric and magnetic field 
energetically allow it. The observed onset of the instability for the 
MPD arc near the center plane of the anode could coincide with the region 
where the mobile electrons begin to cross the magnetic field. This 
criterion appears to he satisfied also for the LHA where the magnetic field 
is perpendicular to the electric field everywhere including the immediate 
vicinity of the cathode in contrast to the MPD arc geometries, therefore, 
in the LHA emitted electrons are forced to cross the magnetic field within 
a short distance from the cathode. More sensitive magnetic field shaping 
and diagnostics are planned to gain a better understanding of the onset 
of rotating instabilities in an MPD arc. 
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ANALYSIS OF PLASMA EXPANSION IN A MAGNETIC NOZZLE ' 
by E. L. Walker, NASA Lewis Research Center 

An analysis is made of the expansion of a plasma through an axisymmetr i c 
magnetic nozzle produced by Helmholtz coils. The plasma is assumed to 
be inviscid and to consist of hot electrons and'cold ions. The steady 
state mass, momentum, and energy conservation equations along with 
Maxwell's equations are solved numerically near the 'axis of symmetry 
for the ratio of the ion exhaust energy to the initial electron thermal 
energy 

"YmZt - 

The solutions (see figure) are given by a family of curves which are 
described by two parameters: (1) the location of the on-axis sonic 

point x = j^-/a where a is the coil radius; and (2) a dimensionless 
magnetic field ^ which involves an effective electron Hall parameter. 
The solutions show that the inclusion of electron heat flow gives 
values of -T^ax which are significantly higher than that for an adia- 
batic expansion. These results can be used to calculate an upper limit 
on the efficiency vs. specific impulse of a device such as- the small 
MPD arc. 
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PERFORMANCE LIMITS OF ION THRUSTERS AND LOW POWER MPD ARC THRUSTERS 
by George R. Seikel, NASA Lewis Research Center 

By examining a power balance for these DC discharges an unders tandi ng 
of the plasma processes that limit their performance can be obtained. 

If only unavoidable discharge losses are considered then the limiting 
performance is determined. The two major unavoidable losses i n a DC 
discharge are the power carried to the anode by the current and the 
power consumed in the ionization and competing excitation processes 
in the plasma volume. 

/ 

If only the minimum anode power and ionization power losses are con- 
sidered then the limiting performance of a DC bombardment ion thruster 
can be evaluated. The limiting efficiency vs. specific impulse is a 
function of only the plasma discharge energy expended per AMU of beam 
ion produced (i.e. eV/ion-AMU). if the discharge has a constant elec- 
tron kinetic temperature, then for a given propellant, discharge volt- 
age, and anode work function, the discharge energy per beam ion AMU is 
only a function of electron temperature. For any discharge voltage 
there is an optimum discharge electron temperature. The calculated 
minimum values of eV/ion indicate that present 36 V Hg ion sources may 
be approaching a fundamental limit. This analysis also indicates per- 
formance gains possible through raising the discharge voltage and the 
potential performance of thrusters using propellants other than Hg. 

In a similar manner a power balance for the low power MPD arc thruster 
results in predicting the limiting efficiency vs. specific impulse. 

This analysis makes use of the results of a plasma expansion analysis 
(to be presented by Walker) and considers the minimum anode and volume 
ionization losses. The theoretical performance limit so calculated is 
compared with present NASA-Lewis low power xenon MPD arc thruster data. 
The discrepancy between the theoretical limit and present performance 
indicates that substantial performance improvement may be possible. 

In making these theoretical limit calculations both the magnetic field 
and sonic point location were optimized. If the sonic point were 
actually further downstream then the theoretical limit would be lower. 
This point brings up one of the big questions yet to be answered for 
such flows. Where does the sonic point occur and why? 
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EXPERIMENTAL INVESTIGATION OF CURRENT ROTATIONS IN A MPD ARC 
ROBERT A. COCHRAN, Research Assistant, and JAMES A. FAY, Professor, 
Department of Mechanical Engineering, Massachusetts Institute of 
Technology, Cambridge, Mass. 

Recent experimental studies Indicate that under certain conditions 
the current in the inner electrode region of a MPD arc is not uniformly 
distributed throughout the inner electrode region but is confined to an 
azimuthally localized region or spoke which rotates in the Amperian 
direction at steady frequencies on the order of 100 Khz.. This study 
experimentally determines the conditions under which this highly non- 
uniform current distribution occurs and investigates the variation in 
the frequency of the conducting region as well as the electrode voltage 
of the MPD arc operating in both the diffuse and spoke mode. The 
experiment was conducted in a pulsed quasi-steady MPD arc of 500 ^sec 
duration. The electrode geometry consisted of two concentric copper rings 
of 1cm axial length having a 1cm radial gap. The inner ring had an 
outer radius of 2cm. The primary propellant used was argon and the arc 
behavior was investigated for large variations in m(. 056-4,0 gm/sec), 
1(100-5000 Amp), and B (0. 1-2.0 telsa) . Additional measurements were also 
made using xenon, helium, and hydrogen as propellants. 

It was found that the single rotating current spoke occurs at inter- 
mediate current levels above a certain magnetic field. strength. The 
discharge is very uniform and diffuse at low current levels and at high 
currents it consists of many small fluctuations each corresponding 
to a small percent of the total current delivered to the arc. A comparison 
of the spoke frequency measurements with the actuator disk theory of Fay 
and Cochran (1) shows that the actuator disk model gives a satisfactory 
prediction of the behavior of the spoke frequency until a limiting frequency 
is reached after which the spoke approaches a constant rotational speed 
which is the same for all gases. The point at which the actuator disk 
model becomes invalid corresponds to a magnetic Reynolds number based on 
the spoke speed of order one for all gases investigated. The variation in 
the electrode voltage was found to obey the correlation of Patrick and 
Schneiderman (2) when the arc was operated in the diffuse mode at low 
current levels. In the spoke mode voltages several times this prediction 
were measured and can be explained in terms of the EMF generated by the 
rotating spoke. 

1. Fay, J.A. , and Cochran, R.A., "An Actuator-Disk Model for Azimuthally 
Nonuniform MPD Arcs," AIAA Journal , Vol, 7, No. 9, Sept. 1969, pp, 1688- 
1692. 

2. Patrick, R.M., and Schneiderman, A.M., "Performance Characteristics of 
a Magnetic Annular Arc," AIAA Journal , Vol. 5, No. 2, Feb. 1966, 

pp. 283-290. ‘ 



JXB ACCELERATOR WORK AT AMES 
James R. Jedlicka 

MSA Ames Research Center, Moffett Field, California 


The present steady-state accelerator, which was installed two 
years ago, receives plasma from a constricted arc. This combination 
offers high plasma entering velocity and electrical conductivity 
without the need for seed. The flow through the accelerator 
approximates fully developed pipe flow, a consequence of its high 
temperature, and results in severe circulating currents upon appli- 
cation of a transverse magnetic field. For many combinations of J 
and B, no net acceleration of the stream occurs. 

The most favorable combination of operating conditions has 
produced a velocity increase of times that for which J=B=0, or 
4 times that before J was applied but with non-zero B. The absolute 
velocity for these conditions, however, was only 8 km/s ; which is 
not spectacular when compared to the 6.3 km/s available from the 
constricted arc alone under other conditions. 


Possible means to improve performance are discussed. 



THE DESIGN OF THE LANGLEY TWENTY MEGAWATT LINEAR 
PLASMA ACCELERATOR FACILITY 

ArleaJE*. Carter 
NASA, Langley Research Center 


A linear, crossed-field plasma accelerator facility for high-speed 
aerodynamic testing and basic research in high-speed magnetohydrodynamics 
has been designed. Use was made of the experience and results obtained 
with previous smaller, lower speed accelerators together with a quasi- 
one- dimensional ^theoretical analysis. General design philosophy centered 
around three criteria: (1) make J x B as large as practicable; (2) operate 

with approximately constant current density; (3) maintain a constant Hall 
potential gradients. Convergent and divergent channel geometries were 
considered and rejected because of, among other things, lack of uniform gas 
properties across the channel cross section and the difficulty of properly 
matching the device to an existing power supply. The final design has a 
constant-area channel designed for constant electrode potential difference 
with the magnetic field tailored to achieve this condition. A compromise 
in the magnetic field strength at the upstream end of the accelerator to 
help assure proper operation of the arc heater resulted in a small decrease 
in the electrode potential difference and an increase in electrode current 
for the first four or five electrodes. Figures 1 and 2 are computer . 
prediction's of velocity and current density based on the modified magnetic 
field distribution shown in figure 3. 

The matching of the accelerator to the existing 10 megawatt power .supply 
is difficult. The power supply has eight separate modules which may be 
connected in any series-parallel combination. The number of accelerator 
electrodes greatly exceeds the number of modules available for power. With- 
out individual power supplies for each set of electrodes an extensive resist- 
ance network is needed to properly divide the current among the many electrodes 
since the applied potential distribution along the anode and cathode walls 
must match the internally generated Hall potential distribution or circulating 
Hall currents will result, degrading the accelerator performance. A special 
resistor network was designed to handle the high power dissipation. 

Figure 4 is a schematic diagram of the facility. The arc heater is shown 
at the left and consists of a cathode, constrictor, anode, plenum, and super- 
sonic nozzle. Nozzle exit conditions are a velocity of approximately 3800 m/s, 
enthalpy of 1.6 x 10^ J/kg, mass flow of- 18 g/s. The working gas is nitrogen 
seeded with 2 percent mole fraction of cesium vapor. One 10 megawatt power 
supply is used to power the arc heater and a second 10 megawatt supply is 
used for the accelerator which mounts between the poles of a large c-magnet. 

The accelerator channel is 6.3 cm square and 50 cm long; it contains 36 pairs 
of individually water-cooled electrodes. The sidewalls are boron nitride 
plates backed by water-cooled copper, coated with beryllium oxide for elec- 
trical insulation. Each water and electrical connection is of a "quick- 
disconnect" type, greatly facilitating the attachment of the 156 water leads 
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and 72 electrical, leads. The accelerator forms its own environmental 
chamber eliminating the need for mounting it within a vacuum tank. 
Immediately downstream of the accelerator is a test section used for flow 
observation and diagnostics. The flow is exhausted into a five-foot 
diameter duct which is cooled by spraying water on its outer surface. 

A four-stage steam ejector provides a vacuum source. The expected accel- 
erator exit conditions are: Velocity of 13,000 m/s, enthalpy of 

11.8 x 10 / J/kg, Mach number of 7.0, density corresponding to 58 km 
altitude, stagnation point heating rate of 7.2 kW/cm^ to a hemispherical 
body 2.54 cm in diameter. 
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OPERATING CHARACTERISTICS AND INITIAL TESTS OF THE 
LANGLEY 20 MW LINEAR PLASMA ACCELERATOR 

W. R. Weaver 

NASA, Langley Research Center 
Hampton, Va. 

The Langley Research Center’s 20 MW Linear Plasma Accelerator Facility 
has been made operational. The initial problems with the arc heater and 
the improvements necessary to make it a suitable plasma source are descri e 
The initial accelerator tests consisted of powering a small number of the 
upstream electrodes, adjusting the resistor network until the Hall currents 
were eliminated, and then connecting more electrodes and adjusting their 
resistors until 30 electrode pairs were connected and the desired current 
through each pair was obtained. The final current distribution, electrode 
potential difference, and the electrode-wall potential distributions are 
given. Pitot pressure measured at the accelerator exit was 0.45 atmosphere, 
which indicated a velocity of 80% of the design velocity of 12,000 meters 
per second. 



ELECTRODE POTENTIAL DISTRIBUTIONS AND POWER SUPPLY CONNECTIONS 
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PRESENT STATUS OF THE EXPERIMENTAL MHD POWER GENERATION PROGRAM 
by R. J. Sovie and L. D. Nichols 
NASA Lewis Research Center 

The NASA-Lewis closed loop MHD facility (Fig. 1) has been operated sev- 
eral times at 3200° F in the past year with no facility problems. Con- 
sequently, we have been able to concentrate on the plasma and generator 
problems. In order to increase the plasma conductivity it was necessary 
to install an impurity purge system and a new cesium injection system. 
The method of pre- ioni zat ion was also changed. The significance of 
these design changes will be discussed briefly. As far as the generator 
data is concerned open circuit voltages and short circuit currents have 
been measured. The effects of internal current leakages on the measured 
generator performance have been ascertained and will be discussed. It 
was also found that a coating is formed on the generator electrodes that 
seriously inhibits the current flow in the external load circuit. Vari- 
ous starting techniques are presently being investigated in order to 
eliminate this coating and these results will also be presented. 
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Figure 1 - Bedplate showing heater, nozzle MHD duct, magnet, and reheater. 
The flow distributor at the heater inlet is not shown. 
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THE EFFECTS OF ELECTROTHERMAL INSTABILITIES ON SEVERAL 
NON-EQUILIBRIUM PLASMAS 

by Allan R. Bishop, NASA Lewis Research Center 

The effects of electrothermal instabilities on several different non-equ? 
librium plasmas is examined. In argon seeded with cesium the effect of 
a molecular impurity magnifies the instability losses. A one-tenth per- 
cent carbon monoxide impurity combined with the instability reduces the 
maximum power density to one-eighth of its ideal value. This maximum 
value occurs at a load parameter of about 0.2. 

The effect of the instability on the output power of a specific generator 
using neon seeded with cesium, is also examined. With generator shape, 
Mach number, pressure drop, and magnetic field specified, the output 
power loss as a function of stagnation temperature is determined. The 
1'oss varies from a factor of- four at 1800° K down to zero at about 2500° 
K, where the instability disappears. 

The feasibility of using helium seeded with uranium as a working fluid 
in generators driven by liquid or gas core nuclear reactors 'is discussed. 
The use of a helium-uranium mixture appears to be a possibility al-though 
the power density is lower than with more conventional working fluids. 

The optimum mixture ratio is about six parts uranium per thousand parts 
helium. Even with the electrothermal instability losses included, the 
non-equilibrium conductivity is higher than the equilibrium conductivity. 
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"HOT" ELECTRODE SHEATH DROPS IN A MHD GENERATOR 
by M. A. Mantenieks, NASA Lewis Research Center 

Electrode sheath drops in an arc heated argon, supersonic MHD generator 
were measured. The Mach 3.5 generator has carbonized 2% Th-W electrodes 
and is capable of generating currents of up to 34 amps and open circuit 
voltages of about 20 volts with a magnetic field of .55 Tesla depending 
on the electrode gap. The other parameters measured In the experiment 
included gas pressure, gas temperature (measured calorimetrical ly) and 
electrode temperature as a function of generated current through the 
electrode. Because of the highly non-equilibrium state of the plasma 
(due to the arc heater and frozen flow) zero current sheaths at the 
anode and cathode are of the order of one volt. The cathode sheath 
drop increases from this value, as the current increases, to 6,5-10.3 
volts depending on the initial cathode temperature and plasma condi- 
tions. Once this level Is reached the sheath drop becomes relatively 
independent of the current. The anode voltage decreases with increas- 
ing current and becomes negative at higher currents (- 4.6 volts at 30 
amps ) . 

The cathode temperature was found to first decrease with current and 
then increase linearly at higher currents. (For example: Temp, drops 

from 2250° K at ! = 0 amps to 2232° K at f =3-5 amps and i ncreases to 
23 10° K at 1 = 22 amps). The anode temperature increases linearly with 
increasing current. A current balance at the electrodes following the 
work of Lyubimov' and an energy balance at the electrode surface give 
qualitative agreement with the experimental results except for the high 
current regime. 


REFERENCES 

1 G. A. Lyubimov, "Layers of Abrupt Change of Potential on 'Hot' Elec- 
trodes," translated from Teplofizika Vysokikh Temperatur, vol . 4, 
no. 1, pp. 120-132, January-February 1966 . 
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ANALYTICAL AND EXPERIMENTAL STUDIES OF MHD GENERATOR 
CATHODES EMITTING IN A "SPOT" MODE 
by Lester D. Nichols and Maris A. Mantenieks, 

NASA Lewis Research Center 

The mode of electron emission from an MHD generator cathode can affect 
the performance of the generator. Under certain conditions the emis- 
sion of the electrons is confined to a small spot on the cathode. These 
spots are similar to the spots observed in' any thermionic arc ( J_) . 

There are many theoretical studies of these arc spots - and they have 
in common the assumption that the spot is formed by Ion energy transfer 
to the cathode. The ions are formed in a region near the cathode as a 
result of electron collisions with gas atoms. The differences among 
the previous studies involve their treatment of the ionization region 
and the number of unknown parameters which must be specified in order 
to use the model to predict current voltage characteristics. Ecker (2) 
determines the ion current from the solution which gives the minimum 
voltage as a function of a "contraction parameter" rather than consider- 
ing the energy conducted from the spot. Bauer (3) and Lee and Greenwood 
(4) leave the ion current as a free parameter to be determined experi- 
mentally, Adams and Robinson (5.) consider a rather complete set of heat 
transfer mechanisms, but require the specification of the angle between 
the current stream tube and the normal to the cathode before comparison 
between their model and experiment can be made. The aim of this study 
is to present a model which requires ho empirical constants. Such a 
theory is presented by Bade and Yos (6 ) . But they make many assumptions 
about the nature of the phenomena to develop a simple theory. The 
present report is an attempt to improve the Bade and Yos theory by using 
fewer assumptions but retaining simplicity. Bade and Yos assume that 
the ions and electrons are at the same temperature, all the atoms are 
ionized, and that the electron current carries no energy from the ion- 
ization region into the discharge. We offer a means of calculating 
electron temperature, ion density, and the energy carried by the elec- 
trons from the ionization region to the discharge volume. In all other 
respects these two models are the same. The results for given gas con- 
ditions and electrode work function predict a sheath- voltage which 
decreases with both increasing current for a fixed undisturbed (i.e., 
far from the spot) cathode temperature, and increasing undisturbed 
cathode temperature at a fixed current. Experimental data is shown 
which supports the conclusions predicted by the model. 


REFERENCES 
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pp. 64-68. 

2 G. Ecker, "Die Stabi 1 isierung des Lichtbogens for Anode und Kathode," 
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USE OF A HIGH-INTENSITY R-F PLASMA RADIANT ENERGY 
SOURCE FOR RESEARCH ON GASEOUS-CORE NUCLEAR REACTORS* 

John S. Kendall, Supervisor Fluid Dynamics 
United Aircraft Research Laboratories 
East Hartford, Connecticut 

ABSTRACT 

The nuclear light "bulb rochet engine concept is based on the transfer 
of energy by thermal radiation from gaseous nuclear fuel through an inter- 
nally cooled transparent wall to seeded hydrogen propellant (Fig. l). In a 
reference engine design, the effective radiating temperature of the nuclear 
fuel is 15,000 R and the hydrogen propellant exit temperature is 12,000 R. 

For this design, the calculated specific impulse is 1870 sec, the engine 
thrust is 92,000 lh, and the engine weight is 70,000 lb. The gaseous nuclear 
fuel is kept away from the transparent wall by a vortex created by the tan- 
gential injection of neon near the inside surface of the transparent wall. 

The neon discharging from the cavity contains gaseous nuclear fuel which is 
then condensed, centrifugally separated from the neon, and pumped hack into 
the fuel-containment region of the vortex. The neon is also pumped hack 
into the cavity. The closed-cycle fuel system leads to the possibility of 
providing perfect containment of the gaseous nuclear fuel and fission 
products . 

Experiments were conducted to simulate the thermal environment in a unit 
cavity of a nuclear light bulb engine . These experiments include develop- 
ment of a high- intensity r-f induction-heated light source, transparent-wall 
model tests, and propellant heating tests. 

The objective of the high- intensity r-f light source research is to 
develop a non-nuclear source having a radiant energy flux in the range of 
interest for the full-scale engine (equivalent black-body radiating tempera- 
tures between 8000 and 15,000 R). This research was conducted using the 
UARL 1.2-megw r-f induction heater to deposit energy in argon plasmas (Fig. 2). 
The configuration used in these experiments consists of two concentric fused 
silica tubes located within the coils of the induction heater; argon is 
injected at one end to produce a vortex flow and is exhausted through ports 
located at the centers of the water-cooled copper end walls. Tests were con- 
ducted at pressures between 2 and l6 atm. 

For the maximum-power test condition achieved to date, a total of 216 kv? 
was deposited in the discharge (Fig. 3); of this, approximately 156 kw 


* This research was supported by the Joint AEC-NASA Space Nuclear Propulsion 
Office as part of Contracts NASw-8^7 and SNPC-70 and is reported in Refs. 1-3. 
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(75 percent) -was radiated outward through the peripheral wall of the test 
chamber. The corresponding radiant energy flux was 36.8 kw/in. 1 2 , which is 
equivalent to a black-body radiating temperature of 10,200 R . This is well 
into the range of interest for the engine; hence, the development has progres- 
sed to the point where useful experiments related to the thermal environment 
of the engine are being conducted. 

The objective of the transparent-wall model tests is to test models 
(Fig. 4) adjacent to the high- intensity light source to determine their 
operating limits . The heat deposition rate expected in the full-scale engine 
due to thermal radiation from the nuclear fuel at an equivalent black-body 
radiating temperature of 15,000 R is approximately I.65 kw/in .2 of wall sur- 
face area. Tests employing a 0.005-in • wall thickness model resulted in heat 
deposition rates on the inside surface of the model as high as 2.2 kw/in.^ 
(i.e., greater than that expected for the full-scale engine) with no apparent 
damage to the fused silica coolant tubes. 

The objectives of the propellant heating experiments are (l) to develop 
methods for injecting solid particle seeds into a simulated propellant stream, 
(2) to study the effectiveness of buffer layers for preventing wall coating, 
and (3) to conduct experiments in which a large percentage of the incident 
radiation from the light source is absorbed in a simulated propellant stream. 
This work is still in the initial stages. The propellant heating configura- 
tions used are generally similar to the geometries of the components expected 
to he employed in the nuclear light bulb engine. Argon seeded with micron- 
sized carbon particles (instead of hydrogen seeded with tungsten particles) 
was used as the simulated propellant. At the low radiant energy source power 
levels of up to 3 kw that were used in these initial tests, temperature rises 
up to 223 R were obtained. Further increases in simulated propellant tempera- 
ture rise will be obtained primarily ^y means of increased power, improved 
particle deagglomeration to increase absorption, and more effective buffer 
layers to reduce the coating of particles on the transparent walls . 

REFERENCES 
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SKETCH OF RADIANT ENERGY SOURCE CONFIGURATION USED IN TESTS 
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FIG. 3 

VARIATION OF RADIANT ENERGY FLUX WITH TOTAL R-F DISCHARGE POWER 
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URANIUM PLASMA RESEARCH AT THE UNIVERSITY OF FLORIDA 
R. T. Schneider, C. D. Kylstra and M. J. Ohanian 
University of Florida 
Gainesville, Florida 

The eventual successful operation of a nuclear energy source utilizing fuel 
in the plasma state, whether for propulsion or other applications, will require 
knowledge of the radiation emitted by the uranium plasma as a function of 
operating pressures, temperatures, wavelengths, and system geometry. Supple- 
mentary to this is the composition of the plasma, i.e., the partial pressure 
distributions of the various uranium Ionic species and the electrons, as well 
as that of any working fluid present. Calculations concerning criticality, 
radiative heat transfer to the working fluid, and system startup and control 
all need the above information. 

The object of the uranium plasma research program at the University of 
Florida is to acquire at least some of the radiation emission and composition 
data needed in addition to some nuclear data (Ref. 1-5). This' paper will be 
concerned with the preliminary experimental results of the measurement of 
temperature and partial pressure of a uranium plasma. (For other research 
efforts on radiation and composition of uranium plasmas see Ref. 6 and 7.) 

The experimental apparatus is a high pressure chamber capable of 100 
atmosphere pressure with water cooling to the thermal shield, tungsten cathode, 
and tungsten crucible with a uranium insert as the anode, the uranium is 
vaporized by a D.C. arc. 

The measurement method used for the arc temperatures and temperature 
distributions is the Boltzmann Plot technique and the single-line relative 
intensity method which can be used when the radial intensity distribution 
shows an off-axis maximum. 

The partial pressure distributions of the individual uranium ionic species 
can be obtained from a correlation of the temperature distribution from the 
Boltzmann Plot techniques and the temperature distributions from the single- 
line relative intensity method. Fig. 1 shows the result of a temperature deter- 
mination using Boltzmann Plots for singly ionized uranium. -Pig. 2 shows the 
partial pressure distribution for singly ionized uranium (UIl). 

The thermodynamic characteristics of uranium at high temperatures 
(T > 3000°K) are important for several applications, e.g., accident analysis. 

Due to temperature limitations little experimental data is available. One 
important property falling into this category is the vapor pressure. Recent 
work using the effusion method (Ref. 8 and 9), has been limited to a maximum 
temperature of approximately 2500°K. For higher temperatures , the vapor 
pressure must be extrapolated. The boiling point .temperature predicted by 
extrapolation is in excess of 4000^R at one atmosphere. For these conditions, 

uranium Is already a plasma. 

Nottingham (Ref. 10) has shown that for an arc, the voltage- current 
characteristic can be related to the anode boiling temperature. Thus, the 
boiling temperature can be determined by measuring the voltage- current charac- 
teristics . 
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Figure 3 illustrates a typical set of arc voltage- current characteristics. 
Figure 4 shows the uranium boiling temperatures determined in this experiment. 
There is good agreement between the experimentally determined boiling temperatures 
and the boiling points predicted by extrapolation from low temperature vapor 
pressure measurements. The determined boiling temperature is 3,837°K for 1.0 
atmosphere. 

Future plans for the research effort are the determination of the emission 
coefficient of a uranium plasma for a wide range of temperatures and pressures, 
the study of a fissioning plasma and the study of the interaction of fission 
products with a noble gas plasma. 
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COMBINED THEORETICAL AND EXPERIMENTAL RESULTS 
OF HEAT TRANSFER TO AN ANODE* 

P. F. Massier + and T. K. Bose^ 

Jet Propulsion Laboratory 
Pasadena, California 

ABSTRACT 


The purpose of this investigation is to evaluate the heat transfer to an anode 
in the presence of a transverse gas flow. The work was initiated because the 
heat transfer to some of the elements of plasma propulsion devices, especially 
to the anode, can be significantly higher than that encountered in chemical 
propulsion thrust chambers. Furthermore, impingement of electrons on the 
anode can produce an amount of heating greatly exceeding that caused by con- 
vection from the bulk motion of the plasma. Experimental heat transfer results 
are analyzed in conjunction with a theoretical prediction method for ionized 
gas flows over a segmented flat surface with an electric field applied to the 
upstream segment, which is the anode (Fig. 1). Under the influence of the 
transverse gas flow, the arc attaches near the downstream end of the anode, 
and the discharge is deflected from a vertical plasma column into an elongated 
U-shape by the motion of the argon gas giving a boundary layer type appearance 
over the anode and the segments downstream of the anode. Tests have been con- 
ducted in an argon atmosphere at a nominal pressure of 50 mm Hg and a current 
of approximately 40 amps over a gas velocity range of about 55 to 260 ft /sec. 
Each of the electrode segments are individually cooled and the heat transfer 
rate is determined from the coolant flow rate and the rise in the coolant 
temperature. In the absence of plasma free-stream temperature measurements, 
which are planned, the temperature is deduced from heat loads to the three 
segments downstream of the anode. By matching the heat load distribution of 
the downstream segments with the theoretical predictions by considering con- 
vection and diffusion only, it has been established that the origin of the 
thermal boundary layer for all tests occurred at a distance of about 0.8 of 
the segment length (0.20 in.) from the leading edge of the anode. It is also 
deduced, with the help of Fig. 2, that the free stream temperature associated 
with the heat transfer for one that was around 13,000°R, and was almost con- 
stant along the flow direction. From the experimental heat load to the anode, 
the theoretical heat load to the anode due to convection and diffusion is 
subtracted to obtain the contribution of the electric current, which was 0.247 
Btu/sec or 87 percent of the total heat load. This was then used in conjunction 
with the theoretical prediction as shown in Fig. 2 to estimate the combined 
anode fall potential and work function which was found to be 4.9 volts. This 
value seems quite reasonable compared with measurements of other investigations. 


This work represents the results of one phase of research carried out in the 
Propulsion Research and Advanced Concepts Section of the Jet Propulsion 
Laboratory, California Institute of Technology, under Contract NAS7-100, 
sponsored by the National Aeronautics and Space Administration. 

’^NASA-NRC Resident Research Associate 
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A STUDY OF THE EFFECT OF ION 
INTERACTIONS IN PLASMA ACCELERATORS 

Sol Aisenberg and K. W. Chang 
Space Sciences, Incorporated 
Waltham, Massachusetts 


As part of the study of plasma boundary interactions in magnetically 
accelerated plasmas, it was established that as much as 15 percent of the 
plasma driving force was transferred by the positive ions to the cathode. 

This ion current,drag at the cathode contributes to the cathode damage as 
well as to the reduced plasma propulsion efficiency. The study of the trans- 
lational cathode drag due to the positive ions (and to the neutral gas compo- 
nent as well) indicated that there is a threshold in the tangential drag force 
as well as a saturation effect. The threshold can be ascribed either to a 
threshold in the tangential ion velocity, in the ion current partitioning, or 
in the tangential accommodation coefficient. In studying the effect of ion 
currents on the cathode it was established that there is a need for an accu- 
rate measurement of the ion and plasma velocity. With this velocity meas- 
urement as a function of magnetic field and current it is possible to determine 
if there is, indeed, a threshold in the tangential velocity. At the same time, 
the measurement of the velocity together with the drag measurement can per- 
mit one to deduce the ion current partitioning at the cathode. The determi- 
nation of the ion current partitioning at the cathode spot will also be of help 
in the establishment of a model for the electron emission at arc cathodes. 

The recent effort, therefore, has been devoted to the determination of 
the plasma velocity in order to determine the 1) ion current fraction at the 
cathode spot, 2) any possible threshold in ion tangential velocity, and 3) the 
existence of a saturation effect in ion velocity. 

There are a number of different ways of measuring ion velocity: one 
is the back EMF method where the terminal voltage is measured as a function 
of transverse magnetic field for a constant arc current and gap spacing. A 
second method of determination of plasma velocity is by means of a time-of- 
flight measurement where the velocity of a moving plasma is measured by 
means of photo-multipliers, and is an unsatisfactory method of determining 
plasma velocity since it actually measures the velocity of propagation of the 
plasma phase rather than the velocity of the individual plasma components. 
This is shown by the existence of the phenomenon of retrograde motion. 

Measurements were made in a circular accelerator channel to compare 
the velocity deduced from back EMF with the velocity deduced by time-of- 
f light measurements for the identical J x B accelerator situation. The analysis 
of Alfven suggests the existence of an ionization limiting velocity. Measure- 
ments have shown that the back EMF determination of plasma velocity is much 



closer to the values expected from the ionization limiting velocity and that 
the time- of-f light method gives about two orders of magnitude of lower ve- 
locity and is, therefore, unsuited for velocity determination. The ionization 
limiting velocity as determined by the theory' of Alfven is predicted to be a 
constant independent of magnetic field. Some results are illustrated in 
Figure 1. The velocity deduced from the back EMF method is essentially a 
constant since it is related to the slope of terminal voltage with respect to 
magnetic field and it was observed that the slope is essentially a constant. 

The time -of- flight method indicates a slight increase of transit velocity as 
a function of magnetic field. The concept of ionization limiting velocity sug- 
gests that this velocity is independent of pressure while measurements of the 
velocity determined from back EMF and also from time -of-f light show a de- 
crease of velocity with increasing pressure. 

As part of this research program a technique has been developed for 
the direct measurement of plasma ions and atoms by means of a Doppler 
shift technique. The system is illustrated in Figure 2, Light from the plasma 
is sampled by means of fiber optics and can be used to measure velocities 
for both approaching and receding plasmas. The wave length shifts of the 
radiation is measured very accurately by the means of a Fabry- Perot inter- 
ferometer. By modulation of the interferometer mirror and by means of a phase 
sensitive detector it is possible' to close the feedback loop in order to op- 
erate the system as an optical discriminator where a resonant frequency of the 
interferometer is locked to the line center. In this way the Doppler shift sig- 
nal consists of a DC output voltage which is directly proportional to the ion 
or atom velocity. Preliminary measurements of the Doppler shift velocity 
indicate that it is comparable with the values determined from the back EMF 
method as well as with the values approximately expected from the ionization 
limiting velocity. Doppler shift measurements are being made for different 
gases and electrodes and as a function of gas pressure, arc current and trans- 
verse magnetic field. 

The role of positive ions in the cathode spot has also been studied. 

It is well known that in the arc discharge, as the gas pressure is reduced, the 
current density at the cathode decreases gradually until at a lower pressure, a 
transition phase is reached, after which the current density suddenly increases 
to very high values. This transition is believed to be a change from a thermionic 
cathode to a field type of cathode. The usual model of field emission, which is 
based on the electrostatic field of the average ion space charges in the bipolar 
sheath, is inadequate to explain the observed level of emission current den- 
sities. A new model based on a modified field emission process which takes 
into account the local electric field associated with each individual ion is 
presented. Preliminary calculations show that this added effect due to the 
intense local ionic force field can indeed provide an emission current large 
enough to account for the observed current densities. Experimental verification 
of this model requires an experimental determination of the positive ion current 
fraction at the cathode. 
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Despite recent advances in the area of arc and plasma research, 
there has been no satisfactory way of measuring the fraction of electron 
current at the cathode surface. Although this information is vital for the 
understanding of the emission mechanisms, only very crude estimates are 
available in the literature. These estimates are usually obtained from 
energy balance calculations, and they are inaccurate because of the un- 
certainty in various unknown parameters. The present results of the tan- 
gential drag measurements and the Doppler shift determination of velocity 
will permit evaluation of the ion current fraction at the cathode spot. An 
additional method for the determination of current partitioning at the cathode 
is proposed, based on the principle' of perpendicular force and momentum 
balance at the cathode surface. A unique model for various components of 
electrode forces is presented which is capable of explaining the known 
behavior of cathode forces at different pressure and current levels. The 
fraction of electron current at the cathode can be determined from the fol- 
lowing equation: 


I F + *> e kT e A - r evap m c A 

f= — s 1 - ■ — 

1 (I/e) (1 - V ) V Zm.e | 4> cd j 


Where I is the measured total arc current, I , the electron current at the 
cathode, tr , the normal momentum accommodation coefficient, e,the 
electronic charge, m^, the ionic mass, ^ cc j, the cathode drop, F, the 
measured resultant force on the cathode, N g , the electron number density 
of the arc spot, T e , the electron temperature of the arc spot,F , the 
flux of evaporating cathode material, T c , the cathode temperatur^f m c , the 
molecular mass of the cathode material. A, the cathode spot area and k, 
the Boltzmann constant. 

Because in the low pressure range of interest the conditions of 

F>>T \ / 2kT_ n and F>>nJdIhA are usually satisfied (as a simple, 

evap — — A e e 

c 

order of magnitude calculation indicates), the present method is believed 
to be effective and accurate especially in the low pressure range where the 
arc is supported by the ion microfield emission process. 
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MED BOUNDARY LAYERS WITH NONEQUILIBRIUM 
IONIZATION AND FINITE BATES 
Arthur Sherman^ Hsuan Yet?; and Eli Reshotko 2 


During a previous MSA research contract an analysis "was made of 
the nonequilibrium plasma boundary layer developing over an insulator 
wall including imbedded electrodes. Such a case was chosen to simulate 
the boundary layer in an MED channel with finite electrodes. The pri- 
mary difficulty -in that analysis proved to be the singularity at the 
initial and trailing edges of the electrode which made the numerical 
solution tedious. In the earlier analysis we had assumed the flow to 
be in quasi-equilibrium such that the electron density could be calcu- 
lated from the Saha Law using the electron temperature. The primary 
objective of tbe present analysis is to extend the previous study to the 
situation where ionization and recombination rates are slow enough to cause 
significant departures from quasi -equilibrium,, In other words, we wish to 
consider finite rates and study flows which are to some extent frozen. 

In the earlier work, the electron temperature was taken into account 
by considering the partial differential equation governing the conservation 
of electron energy. In order to provide the wall boundary condition the. 
plasma sheath was considered in a simple manner. For the present investi- 
gation it is necessary to include one additional equation governing the 
electron density, which is the partial differential equation for the con- 
servation of electron species. Including this additional equation means 
that we are solving four highly nonlinear partial differential equations by 
the finite difference method we have developed. This has the effect of 
substantially slowing down an already tedious numerical calculation, since 
now k- X k matrices have to be inverted at each calculation point rather 
than 3 X 3 as before. 

To date, the new system of equations allowing for finite rates have 
"been formulated and their transformation completed, Tbe transformed equa- 
tions have also been reduced to finite difference format in a much more 
compact form than had been done before. The hope here is that this new 
form will permit a reduction in computer time. A computer program has 
been written to solve these equations in their finite difference form, and 
calculations are in progress. 


^Computer and Applied Sciences, Inc. 
2 

Case Western Reserve University 
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STRUCTURE OP A COLLISIONLESS SHEATH IN THE 
PRESENCE OF A MAGNETIC FIELD* 


H. A. Hassan 


North Carolina State University, Raleigh, North Carolina 27607 


The objective of this Grant is to study the surface-sheath-plasma interactions 
in devices being considered for plasma acceleration and MHD power generation. In 
such devices the electron temperature is different from the heavy particle temperature 


and one needs to consider the structure of the sheath near the electrode to ob- 
tain necessary boundary conditions. The description of the plasma in a given 
device is then obtained by matching the sheath solution and the solution of the 
conservation equations of the various plasma constituents at the sheath edge. 

The appropriate transport properties for a plasma in a magnetic field and 
the influence of the surface parameters on the sheath structure were studied 
under this Grant^ - ^. At present, we are in the process of integrating these 
studies so that we may obtain an understanding of the coupling between surface, 
sheath and the bulk of the plasma. The work is in the computational stage. 

The work carried out in References 2-b assumes that the sheath thickness 
is less than the mean free paths and the Larmor radii. When the restriction 
on the Larmor radii is relaxed, the equations governing the collisionless sheath 
are the Vlasov equation and the appropriate Maxwell's equations. These equations 
can be expressed as 


3f . 
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B - V x A , V A = - y Q J 


( 2 ) 


E = - V $ , 



(3) 


with 

J = E e. J f . df. , P' =E e. n. = l e, 1 f . df. (4) 

0 J J J - 0 J 3 3 0 


and , e. and m. are the velocity, charge and mass of species jj E and B are the 
electric ^.nd magnetic fields j <f> and are the electric and magnetic potentials, and 
y and s are the permeability and permittivity of free space. For the case of 
an infinite planar electrode with uniform work function, the electric and magnetic 
potentials and the distribution functions are functions of the coordinates normal 
to the electrode. Solution of the Lagrangian subsidiary equations gives the constants 
of the motion which can be written as 


* Work Supported, in part, by Grant NGR 3^-002-0^8. 
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where U , V and W are the velocity components in the x, y and z directions. If 
one asstimes that particles far from the electrode have no velocity components paral- 
lel to the electrode, then the distribution function may "be chosen as 
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where subscript » designates conditions at the sheath edge, a. is constant pro- 
portional to the number density, and 6 is the Dirac delta function. 

Examination of the solution shows that in order for the collected particles, 

5=1 5 to reach the electrode, their energy must be greater than e , where 

2 21 

1 ^1 p 

£in = max [ 2 — (A y - A y (»)) + e;L *] (7) 


The quantity is obtained as part of the solution of the two coupled nonlinear 
equations 




u J 
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( 8 ) 


It should be noted that, in this case, one needs not invoke the Bohn $heath criterion 
because the solution is always stable. 

The results of the above analysis were employed to study the anode losses in 
plasma device where the dominant energy transfer to the anode is that due to the 
electron current. It was shown that the heat flux to the anode is given by 

kT e - e J 

q Z j [ + 4> A + W + e] , e = (9) 

where <j> A is the anode drop, W is the work function and j is the current density. 

The results of the computation show that e, which depends on the magnetic field 
strength, is much less than <J> A . Hence, in those “devices where the d omi n an t energy 
transfer to the anode is that due to the electron current, the anode losses are 
independent of the magnetic field . 
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ANALYSIS OF NOISY LANGMUIR PROBE CHARACTERISTICS 
by Roman Krawec, NASA Lewis Research Center 

The theory of current collection by a probe immersed in a plasma was 
developed in by Langmuir and Mott-Smith for the case of a quiescent 

plasma. Unfortunately, plasmas of interest in fusion research are sel- 
dom, if ever, quiescent. 

A static current-voltage characteristic is. generally obtained by filter- 
ing out any noise, and the resulting curve is then analyzed in the normal 
manner to obtain electron temperature. 

It is not generally realized that the variation in noise amplitude as 
probe voltage is changed may possibly be used as a means of estimating 
the electron temperature. In some cases it may even be possible to 
determine whether the electron temperature, number density, or plasma 
potential is the time varying quantity. This note will investigate the 
effect of a sinusoidal variation of the above quantities on peak-to- 
peak probe noise. 

Consider a plane probe and let the electron current to i t be given by 

I = %» A [Sirf* [" ' %T] (V< o) (la) 

I- h nA (•I?*'/ 1 ' ( v *°) (,b) 

where V is the probe voltage measured with respect to plasma potential, 
and the other quantities have their customary meanings. Let the time 
variation in electron temperature, density, and plasma potential be of 
the respective forms: . 

r„ (// Tfc , T, (i f (3 r cos» t) 

K- Vi. (n-n'/n v cosut) • f), (l+(3 h 

V" 'J?* Ot ri/vf . coswt ) -- v pa c i + ^ cts w v 

This allows us to write eqn. (l) in the form a 

i= w[- 

'" he ; e p bVkt ; t^/kT 

A dimensionless peak-to-peak noise is now defined as: 
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where ! max and J m ; n are the maximum and minimum values of eqn . (2) 
respectively, for a given value of probe voltage. 

In this note we look at the peak-to-peak noise for the case when only 
a single parameter is time varying and the perturbation, p , is held 
constant as probe voltage is varied. Three typical cases are shown 
in Fig. 1. Both a time varying density or plasma potential have the 
same effect on the noise in the electron retarding region until probe 
voltage approaches plasma potential. Noise caused by variation of 
density remains constant in the electron accelerating region, while 
that due to variation in plasma potential goes to zero. The slope 
of the against probe voltage gives the correct value of the 

electron temperature in both cases, A variation in plasma tempera- 
ture manifests itself by having a non-linear portion for the elec- 
tron retarding region, and by taking on a fixed value within the 
accelerating region. 

After the particular time-varying parameter is determined, its ampli- 
tude may be found from a measurement of noise amplitude at plasma 
potential as shown in Fig. 2. 
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THE EFFECT OF SPACE CHARGE ON THE OPERATION OF 
ELECTRO-STATIC ENERGY ANALYZERS 

by Richard R. "Woo lie tt 
NASA-Lewis Research Center 


When electrostatic energy analyzers are operated in plasmas, space 
charge inside the analyzer may lead to incorrect ion and electron 
temperature measurements. These space charge effects may become im- 
portant when: (a) the plasma density is high; (b) the. magnetic field 
is high; (c) or the spacing between the analyzer grids is large. The 
potential of the analyzer aperture plate, which is in contact with the 
plasma, also may affect the ion temperature measurements. This can 
also be related to space charge effects. 

According to theory, space charge effects should become less pro- 
nounced as the grid spacing is reduced. To look for this effect, two 
energy analyzers were tested, each with a different grid spacing as 
shown on the analyzer schematic in figure 1. The analyzers were im- 
mersed in a 2 inch diameter plasma of about 10 cm " density . A 
1,000 gauss field was alined with the analyzer axis. The plasma space 
potential was about -15v and the probe floating potential about -40v. 
The analyzer aperture plate was used as a planar Langmuir probe. A 
cylindrical Langmuir probe was also used to determine the plasma 
characteristics . 

The plasma electron temperature, as determined from both the planar 
and the cylindrical Langmuir probes, was 8.6 ev. The cylindrical 
Langmuir probe gave the same electron temperature when its location 
was varied from 3 to 25 mm ahead of the aperture plate. Changing the 
potential of the aperture plate did not change the temperature deter- 
mined by the cylindrical probe in any of its positions. 

The electron temperature was also measured with the analyzers. Ana- 
lyzer #2 measured the same electron temperature as the Langmuir probes. 
The #1 analyzer gave an electron temperature of 14 ev or roughly 1% 
times higher than the Langmuir probes . Since the #2 analyzer had a 
grid spacing of roughly a third of the #1 probe, we conclude that de- 
creasing the grid spacing improved the accuracy in electron temper- 
ature measurement. This trend is consistent with the space charge 
effects . 

Ion temperatures determined by both analyzers are shown on figure 2 
as a function of aperture plate voltage. The effect of varying the 
aperture plate voltage below the space potential is quite dramatic, 
especially for analyzer #1. For analyzer #1, the measured ion temper- 
ature varies by a factor of five over the range of aperture plate 
voltages shown. This variation is reduced to a factor of two in 
analyzer #2. Again we see a result that can be shown to be consis- 
tent with space charge effects, namely, that reducing the grid spacing 
reduces the variation in ion temperature measurement. 
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In all of the tests care was taken to insure cleanliness of all inter 
nal surfaces of the analyzers. If these surfaces are contaminated, 
especially the grids, a large hysteresis effect results. This is 
probably due to charge accumulation on the dirty surfaces. In this 
condition the analyzer would be unreliable. 

There are two well-known one dimensional space charge models that pre 
diet the above mentioned trends. In both models we assume that only 
one species of charged particles is present. In model A the charged 
particles have no initial drift velocity, only thermal motion. In 
model B the particles have an initial uniform drift velocity, but no 
thermal motion. In both models, the potential variation between two 
adjacent grids is determined from the Chi Id -Langmuir space charge ex- 
pressions. Between adjacent grids the charge flow is considered as 
a current flowing between two plates. The ensuing space charge sets 
up a potential barrier that can affect the collected currents . 

To approximate a one-dimensional energy analyzer the aperture plate 
would have to be replaced by a screen. And then if the analyzer were 
operated in a strong magnetic field, the charged particles would be 
tied tightly to the field lines, yielding a good approximation to a 
one-dimensional model . 

In figure 3, results are presented for model A along with results for 
a model without space charge. Figure 3 is a plot of current to the 
gate versus gate voltage. The charged particles had a temperature of 
7 ev and no initial drift velocity. A space charge will be set up 
between the entrance screen and the gate. If an exponential is fit 
to the model A curve we get a temperature -of 40 ev. Therefore, space 
charge results in a temperature that is almost 6 times greater than 
the true value. In our experiments analyzer #1 gave electron temper- 
atures that were about 1% times the Langmuir probe results . Since an 
aperture produces a pencil beam, the one -dimensional approximation 
used in figure 3 is probably too severe. 

The results for model B are given in figure 4. In this case the 
plasma had zero temperature. Grid A was set at 60 volts to separate 
the ions and electrons . The electrons reaching grid A therefore had 
zero temperature, and a uniform drift velocity equivalent to 60 ev. 
Consequently & grid A we have electrons flowing toward the gate. The 
electron density at grid A was assumed to be 10^ enf 3 . The current 
collected on the gate is presented in figure 4 as a function of gate 
voltage. If space charge were not present, the full electron current 
would reach the gate for zero gate voltage as shown for the no-space- 
charge curve. But since space charge does exist between grid A and 
the gate, the knee is displaced. If an exponential is fit to this 
curve a temperature of 38 ev is obtained. The displacement of the 
knee of the curve is a characteristic observed experimentally. 

We have derived expressions for a new one-dimensional model where 
both thermal and drift velocities are included. This was done for 
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both a single species and a two species model- The resulting equation 
requires a solution to an integral non-linear differential equation 
with two point boundary conditions. A solution was obtained by con- 
verting the problem to a single point boundary condition. The single 
boundary condition was the minimum potential point where the electric 
field is zero. Figure 5 compares the results with model A for which 
the drift velocity is zero. We see that the results for the new 
model are considerably different than those of model A. The' differ- 
ence between the curves is large enough to show that the simplifie 
model does not predict the space potential sufficiently accurate to 
yield a correct current versus voltage curve for an electrostatic 
analyzer. Similar results are obtained when a comparison is made 
with model B, 
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CALCULATIONS OF ION-CYCLOTRON WAVE PROPERTIES IN A HOT PLASMA 
by Donald R. Srgman, NASA Lewis Research Center 

In the past, a Fourier integral technique was used to calculate electric 
and magnetic wave fields excited in an infinitely long plasma column by 
the current of a finite length radiofrequency coil. To solve the Fourier 
integral using the residue method, it was necessary to find the roots of 
a complicated function of the axial wave number, • This could be 

done with little dJjLficulty when the plasma was cold with no damping, so 
that the roots of Ilk) = 0 were real. 

However, in a hot plasma (or one dominated by collisions) the terms mak- 
ing up J p(k) ate complex and the roots are likewise complex. In this 
case it Is often difficult to find a numerical technique to locate the 
roots quickly and easily. To avoid this difficulty with roots, one can 
use a Fourier series approach. For a finite length system of length 2L, 
the fields can be written as a Fourier series as follows 

£ Cr, 2,-fc) r 2 § Cr ^)n f ( L >*k € Cost 

n-=o 

The function is either a sine or cosine function of argument 

«rr £/u or •depending on whether the current distribution in the 

coil is asymmetric or symmetric about Z = 0 and on whether open circuit 
or short circuit boundary conditions are assumed at the ends of the 
plasma column (Z - + L) . j:C T '»*Oin is the Fourier transform of the 
electric field evaluated at \\ = mr/ L or k = nir/3 l again depending on 
the end boundary conditions. 

The solutions of the Fourier series and integral techniques become iden- 
tical when the damping is sufficient to damp the normal modes in a dis- 
tance short compared to L. But the important point is that roots of 
T 1 lk)= 0 do not have to be found when using the series method. 

The series method has been used to calculate fields and power transfer 
to the plasma for conditions favorable for the generation of ion-cyclo- 
tron waves (1 .e. to SL^i ). The coil used was a Stix coil and the 
effects of varying ion and electron temperatures were studied. It has 
been found that whenever ion or electron temperature is increased the 
efficiency of transferring energy to the plasma waves is ‘reduced. This 
effect becomes appreciable when the damping length of the natural modes 
is of the order of, or less than, the coil length. When temperatures 
are high the electron density for which optimum coupling occurs is 
higher. In addition, the use of shorter Stix coil wavelengths reduces 
the coupling, while increased ion temperature anisotropy (Tj.i fXt\[ ) 
improves coupling for a fixed Ti*£ . The effect of reduced efficiency 
with increasing temperature may explain why efficiency goes down with 
increasing coil current in experiments at the Lewis Research Center. 



MODEL FOR THE CALCULATIONS 




PEAK POWER TRANSFER TO PLASMA VERSUS ELECTRON 
DENSITY WITH ION-TEMPERATURE AS A PARAMETER 




POWER TRANSFER TO PLASMA VERSUS 
ION-CYCLOTRON FREQUENCY WITH > ION < 
TEMPERATURE AS A PARAMETER 
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THERMAL I ZATI ON OF ION-CYCLOTRON WAVES IN A MAGNETIC BEACH 
by Clyde C. Swett* NASA Lewis Research Center 

An experimental program is being conducted to study the ion-cyclotron- 
wave method of producing high ion temperatures in a hydrogen or deute- 
rium plasma. In this method waves are propagated in the plasma along 
an axial magnetic field into a region of slightly lower magnetic field 
(called the magnetic-beach region) where the wave energy is converted 
into thermal energy of the ions. This method is capable of heating 
ions to the kilo-electron-volt range. 

The magnetic-field and vacuum-chamber configuration used is shown in 
Figure 1. It consists of an rf center-section with two side-arms in 
a magnetic-mi rror geometry. A kS kW rf transmitter* operated at 6.5 
megahertz* supplies a constant current to a Stix coil. This coi 1 is 
wound so that ion-cyclotron waves propagate out both of .its ends. The 
plasma* although initially started by a hot-cathode discharge* is main- 
tained by the rf alone when metal grids are located at the peaks of the 
mirrors. There are three ways of operating the magnetic field: flat 

field* near beach* or far beach. In the beach regions the wave fre- 
quency will match the ion.-cyclotron frequency and the wave will be 
thermal? zed. . 

In previous experiments the capability of transferring rf power effi- 
ciently to the wave has been shown* the wave has been identified by 
magnetic probing* and the wave has been damped in a beach. The purpose 
of the recent experiments has been to measure the ion temperature result- 
ing from- damping of the wave. ion temperatures were determined by means 
of coils located in the near-beach region which measured the diamagnetism 
of the plasma when the transmitter was unkeyed. 

Ion temperature measurements are shown in Figure 2 for an electron den- 
sity of 2x1012 cm - 3 * deuterium gas* and maximum power. The wave is 
propagating from left to right. The magnetic field of the near beach 
is shown in the lower curve and the temperatures measured by 10 diamag- 
netic coils are shown by the points in the upper curve. The points at 
which the operating frequency is equal to the ion-cyclotron frequency 
is shown as SL - 1 . As the wave propagates into the beach it becomes 
damped and should be completely damped at the position^. = 1. Because 
of the local mirrors created by the beach, the ions become trapped. The 
data show that the highest transverse ion temperature (■'*"'500 ev) is at 
the position of minimum field. The neutron flux measured at the beach 
was compared with theoretical predictions to verify the temperatures 
measured. 
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OBSERVATION OF ION HEATING IN A PLASMA SUBJECT TO THE 
CONTINUITY EQUATION OSCILLATION 
by J. Reece Roth, NASA Lewis Research Center 

A photograph of the modified Penning discharge used in this experiment 
is shown in Fig. 1. The plasma is confined in a superconducting mag- 
netic mirror apparatus with a mirror, ratio of 2.5:1, and a maximum 
field on the axis that can be varied up to two Tesla. The plasma is 
15 cm in diameter at the midplane, and can be operated in the steady 
state for several hours.^ When a high DC potential is applied to the 
anode ring, thermalized ions of high kinetic temperature are observed 
coming out through the magnetic mirrors ^ These ions have a quite 
accurately Maxwellian energy distribution along a radius in velocity 
space. 

These high ion temperatures were investigated further, with the results 
shown in Figure 2. The ion kinetic temperature is shown plotted as a 
function of direct-current anode voltage for four background pressures 
of deuterium gas. The kinetic temperature is approximately linearly 
proportional to the anode voltage up to the limit of the power supply 
used. This suggests that ions of any desired energy can be produced 
by going to a sufficiently high anode voltage. The only other steady- 
state heating system to produce such high ion temperatures is the Burn- 
out-V device at Oak Ridge National Laboratory. 

A retarding potential energy analyzer was used to measure the efflux of 
ions as a function of time. Results are shown in Figure 3. The zero 
level is the baseline between pulses. The charged-particle efflux dis- 
plays pulsar- like sharp, narrow peaks between broad, flat minima. For 
different plasma conditions, these pulses display fine structure, 
periodic interpulses and other interesting detailed features. The 
particle efflux, the light output, and the electrostatic fluctuations 
were in phase at the same frequency. The light output was monitored 
at two axial and azimuthal stations with a pair of photomultipliers, 
and showed that the plasma was oscillating in unison itfith a single 
phase and frequency. Parametric variation of the plasma properties 
showed that the frequency of this oscillation depended on both the 
electron and neutral number densities, but not on the magnetic field. 

These characteristics suggested that the oscillation in question is a 
new phenomenon, and this subsequently proved to be the case. This 
oscillation was labeled the continuity-equation oscillation, because it 
arises as a consequence of periodic solutions to the coupled pair of 
continuity equations for neutrals and charged particles in a slightly 
ionized gas. The frequency is proportional to the square root of 
the product of electron and neutral number densities, and is independent 
of magnetic field. 3 

The oscillation frequency was systematically studied as a function of 
the electron and neutral number densities . Figure 4 shows the results . 
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The observed frequency is plotted as a function of the product of rela- 
tive electron and neutral number density. The square root dependence 
predicted theoretically holds over a range of 5 orders of magnitude in 
the density product.- 5 

The following table summarizes the best temperatures, densities, and 
confinement times simultaneously observed in this experiment thus far: 


Ion Kinetic Temperature, T^, kev 5 

Electron Kinetic Temperature, T e , eV 200 

Plasma Number Density, n e , particle [vr 2x10^ 

Ion Confinement Time, sec 27 


The ion temperatures listed are gratifying, but the densities and con- 
finement times are quite low. These are now limited by the rapid loss 
of ions out the mirrors. Confinement time and density can be improved 
with a closed configuration such as the bumpy torus. This consists of 
12 magnetic mirrors end-to-end in a toroidal array. Such a supercon- 
ducting magnet facility is presently under consideration for future 
use in this experiment. The modified Penning discharge scheme will 
be applied to this configuration. It is one of the few steady- state 
plasma injection and heating schemes capable of being applied to a 
toroidal geometry. A reasonable goal for the bumpy torus would be an 
increase of a factor of 10 in confinement time and number density. 
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FINITE BETA INSTABILITIES PROPAGATING PERPENDICULAR 
TO A MAGNETIC FIELD 

by Lewis C. Himmell, NASA Lewis Research Center 

A theoretical analysis was performed to determine the effect of finite 
beta on electromagnetic modes propagating perpendicular to a uniform 
magnetic field in a fully ionized plasma described by a loss-cone 
distribution. The theory includes the effects of all components of 
the four vector potential, ft was found that finite beta creates new 
resonant mode instabilities which propagate at frequencies that are 
very close to the ion cyclotron frequency and its harmonics. The 
theory was applied to the Dory-Guest distribution functions which are 
of the form 

I 

ii 

and predicted unstable waves with growth rates as large as a f.ew percent 
of the ion cyclotron frequency for a value of beta of one-half. 




The dispersion relation describing resonant modes near the ion cyclotron 
frequency for example is of the following form: 

•a 1 _ a r <*« 0 o) + not) i 

c x a\L x-i (x-i)Pt*)J (2) 

whsr6 ^ ^ ^ 2. / -v 

(> _ ktyt n A — W 4n<A X ~ ci ) 

The propagation vector is denoted by k. whi le tO and C refer to 

the ion plasma frequency, the ion cyclotron frequency, the complex fre- 
quency of the wave and the speed of light respectively. The term involv- 
ing O^ooCO’S the usual electrostatic contribution to the dispersion 
relation while the second term on the right side of the equation gives 
the finite beta effect. Since the functions DC*) and HC*) are second 
and thi-rd degree polynomials in x respectively, Eqn . (2) can be solved 
analytica 1 ly for x. 


An analysis was made of the distributions corresponding to the first 
four values of j for 0.5. When " 0*5 for example, the j = 1 and 
j = h distributions are stable while the j = 2 distribution is unstable. 
Growth rates of the order of a few percent were found for $ - 0.5. The 
j = 3 distribution was found to be unstable for A 4 1.0 X 10“3. Larger 
values of ^ make this distribution stable, V 

Table I lists growth rates in units of Mel as a function of ft for the 
j - 2 distribution. The "resonant" approximation is seen to be a good 
one since is so close to unity. 



TABLE I 


GROWTH RATES FOR j = 2 DISTRIBUTION IN UNITS OF ION 



CYCLOTRON FREQUENCY 


f 5 

1 - ^ /cod 

b-'l/aJ c t ' 

0.5 

- 5.2 X lo' 2 

3.12 X 10" 2 

0 . 1 

- 6.2 X 10“ 3 

1.84 X 10 -2 

0.01 

- 1.53 X 10“ 4 

2.04 X 10“ 3 

0.001 

- 9.85 X 10~ 6 

2.05 x 10 "^ 

0.0001 

- 9-35 X 10" 7 

2.05 X 10-5 

0.00001 

- 9.30 x 10" 8 

2.05 X 10-6 

0.000001 

- 9.30 X io“ 9 

2.99 X 10-7 
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INSTABILITIES AND PLASMA TURBULENCE STUDIED USING CORRELATION TECHNIQUES. 

Jack H. Noon., Rensselaer Polytechnic Institute, Troy, N. Y. 

A hollow cathode gas-fed arc facility has been used at Rensselaer to 
carry out correlation measurements on low frequency instabilities in a 
highly ionized argon plasma. Study of the higher frequency oscillations 
which are not of a coherent nature is planned. 


A schematic of the arc and associated instrumentation is shown in Fig. 
1. The actual construction is easily varied for different experiments. The 
discharge vessel is fabricated of 6 inch inner diameter pyrex pipe sections. 
The cathode is a hollow -barfated tungsten cylinder. Gas can be introduced 
through cathode or anode. The electrode separation can be adjusted over a 
range of one meter. Axial and radially movable Langmuir probes measure 
steady state plasma parameters and perturbations in charged particle density 
and plasma potential. An axial magnetic field of up to 2000 gauss is sup- 
plied using 12" inner diameter solenoid coils. By powering the coil in the 
neighborhood of the cathode independently, the magnetic field in this region 

B can be varied with respect to the main magnetic field B„. 
c h 


Fig. 2 shows how the amplitude (V) and frequency (f) of the dominant 
coherent instability (an ion acoustic wave) are affected by B and p. This 


wave can be suppressed entirely by increasing p or by increasing B rela 
tive to B h - * ' " "" ' ‘ ' ' C 

at Langley. 


A similar effect has been reported in the MPD arc by workers 
Fig. 3 shows changes in floating potential and in plasma 
potential V (calculated knowing the electron temperature; which occur when 
the wave is^suppressed. The plasma potential must be positive with respect 
to the cathode for onset of the instability. The enhanced loss of electrons 
from the core could be due to either increased axial losses or radial losses, 
but all our data is consistent with an increased rate of radial plasma 
transport . 


Another coherent instability (a drift wave) close in frequency to the 
ion acoustic wave appears to have little effect on plasma losses. This can 
be explained in terms of the phase difference 0 between the perturbed 
charged particle density and the fluctuations in plasma potential. A dif- 
ference of k/ 2 is found between the two cases. A theoretical dispersion 
relation has been developed to predict the growth rate and frequency depen- 
dence of these oscillations as a function of B and p and to explain the 
effect of gas pressure on 0. 


Bohm or anomalous diffusion is usually attributed to fluctuating elec- 
tric fields which are not of a coherent nature. To estab l ish the effect of 
general plasma turbulence requires a wider range of pressure and magnetic 
field than those studied up to the present, and the use of both auto-corre- 
lation and cross- correlation measurements of axially separated probes. It 
is hoped to establish by these measurements, correlating the spectra of 
density fluctuations and electrostatic field fluctuations, whether a cas- 
cade in wave number space leads to a universal turbulent spectrum, irrespec- 
tive of the means used to excite turbulence in the plasma medium. A tech- 
nique for measuring the form of the energy distribution of the charged par- 
ticles is also under investigation. 
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THEORY FOR PLASMA INSTABILITIES IN MPD ACCELERATORS 
AND BARIUM CLOUDS 


John J. Kim 

NASA, Langley Research Center 
Hampton, Va. 


MPD Accelerators .- This study is concerned with a theoretical evalua- 
tion of the experiments for onset at a critical magnetic field of rotating 
disturbances in the Linear Hall Current Accelerator-'- (LHA) and the MPD Arc.^’^ 
The evaluation mainly concerns the theories for the E x B instability^*^ and 
the electrothermal instability.® The applicability of other mechanisms ^ 5 ® 
and acoustic wave instabilities have not yet been explored. The E x B 
instability is attractive since in the LHA where measurements inside the 
electrode region could be performed the onset of the instability occurs in a 
region of large density gradient (near the cathode) . It is also important 
that the necessary condition, E • dn /dx > 0, for the onset of the E x B 
instability seems to be satisfied? although further measurements are 
required to obtain full assurance of the electric field E near the cathode. 


On the other hand the frequency characteristics of the rotating disturbance 
as .a function of the applied magnetic field seems to be in better agreement 
with the electrothermal instability.® The variation of onset with mass flow 


(or pressure) and molecular weight is approximately satisfied for both types 
of instabilities. Another criterion for the distinction between the 


two instabilities is the observation of a helical type of disturbance in the 


LHA. Sufficiently detailed measurements inside the electrode region of the 
MPD arc would have to he performed to identify a comparable helical shape. 

The theories for the E x B instability^*® do not explicitly predict a helix, 
although an intuitive nonlinear generalization of the analysis appears to be 
consistent with this result. The electrothermal instability theory as applied 
to the MPD arc,® however, does predict a helical instability in the following 
manner. The theory® deduces two disturbance modes, one purely rotating and 
the other a mixed mode which combines a rotational and a radial wave in the 


direction of the electric field. The analysis of the mixed mode is not 
pursued in the analysis® to its fullest extent since the experiment for the 
interelectrode region of the MPD arc only shows the rotational disturbance. 
(Note: The helical disturbance observed in the exhaust of the MPD arc has 


its axis in the flow direction hut not in the radial direction.) The helical 
disturbance observed in the LHA can be interpreted as a mixed mode consisting 
of a rotational wave and axial wave in the electric field direction. It is 
also of importance that only the mixed mode® analysis predicts the onset of 
instability, whereas the single rotational mode analysis does not yield 
an onset. The experimental results from the LHA-*- thus permit a more detailed 
check of the theory for electrothermal instability .6 The theory® gives the 
correct pitch of the helix. Experiments and electrothermal instability 
theories for striations observed in rectangular geometries^ (related to MHD 
generators) need to be further analyzed to evaluate their relevance to the 
experiments in the LHA and the MPD are. 


*NRC-NASA Resident Research Associate 
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Another means for comparison between the two instabilities is tinder 
investigation which concerns their growth rates and motion as predicted 
by linear analysis. The growth rates and motion of disturbances will be 
compared with the results of pulsed experiments. The nonlinear post onset 
conditions, involving arc columns or spokes, of the E x B and electro- 
thermal and other instabilities must be evaluated in detail to permit a 
better comparison with experiments. 


Barium Clouds .- Attempts have been made to explain the striation phenom- 
enon in barium clouds released in the upper atmosphere on the basis of the 
dissipative E x B instability as developed by Simon ^ which involves collisions 
with neutrals. An application of the E x B instability to the barium cloud 
has recently been reported. ^ Striations have, however, also been observed at 
very high altitudes (12 to 13 earth radii) in a recent experiment. It is 
suggested that a possible origin of these high altitude striations could be 
the drift wave instabilities^ which do not depend on collisions with neutrals. 
These instabilities have been thoroughly studied experimentally and theoreti- 
cally in the fusion program. The flute-like disturbances due to these drift 
wave instabilities could be used for interpretation of the striations in the 
barium clouds. It is then of special interest that the wave fronts of drift 
wave instabilities may be alined with the magnetic field or have an angle to 
it. The experimental observations need to be refined to evaluate their angle 
with respect to the magnetic field. The extension of the E x B instability 
to the collisionless case (including inertial effects) will also be explored. 


The drift motion of the barium cloud is primarily used to study strength 
and direction of electric fields in the magnetosphere. In the cases where 
striations are found the cloud drift motion (and the electric fields) would 
most likely be judged by the motion of the striations. The angle of the 
striations with respect to the magnetic field is then of considerable import- 
ance for the interpretation of electric fields. The possible relation of 
drift wave instabilities to auroral striation phenomena will also be discussed. 
Acknowledgment is made of discussions with F. Allario , R. V. Hess, 0. Jarrett 
and B. D. Sidney. 


REFERENCES 


MED Accelerators: 

1. Sidney, B. D. ; Allario, F. ; Hess, R. V.: Onset of Rotating Disturbance 

in Linear Hall Current Accelerator. AIAA Paper N. 69—381. AIAA Seventh 
Electric Propulsion Conference, Williamsburg, Va. , Mar. 3-5, 1969. 

2. Allario, F. ; Jarrett, 0., Jr.; and Hess, R. V.: Onset of Rotating 

Disturbance in the Interelectrode Region and Exhaust Jet of an MPD Arc. 
AIAA Paper 69-232. AIAA Seventh Electric Propulsion Conference, 
Williamsburg, Va. , Mar. 3-5, 1969. 



114 


3. Jarrett, 0.; Allario, F. ; Sidney, B. D. and Hess, R. V* : Experiments 

for Onset of Rotating Disturbance in Linear Hall Current Accelerator 
and MPD Arc. 6th Intercenter and Contractor's Conference on Plasma 
Physics, NASA Langley Research Center, Dec. 8-10, 1969. 

4. Simon, A. : Instability of a Partially Ionized Plasma in Crossed Elec- 

tric and Magnetic Fields, Phys. of Fluids, vol. 6, no. 3, March 1963, 
pp. 382-388. 

5. Hassan, H. A. and Thompson, C. C. : Onset of Instabilities in Coaxial 

Hall Current Accelerators. AIAA Seventh Electric Propulsion Conference, 
Williamsburg, Va. , Mar. 3-5, 1969. 

6. Smith, J. M. : Electrothermal Instability - An Explanation of the MPD 

Arc-Thrustor Rotating- Spoke Phenomenon. NASA TM X-52535, Technical 
paper proposed for presentation at Seventh Electric Propulsion Conference, 
Williamsburg, Va. , Mar. 3-5, 1969. 

7. Fay, J. A. and Cochran, R. A. : An Actuator Disc Model for Azimuthally 

Nonuniform MPD Arcs. AIAA J., vol. 7, no. 9, September 1969. 

8. Roth, R. : Private communication through F. Allario. 

9. Nedospasov, A. Y. : Striations- Vol, 11, no. 2, Soviet Phys. USPEK HI 

(1968), pp. 174-187 

Velikhov, E. P.; Degtyarev, L. M, ; Samarskii, A. A.; and Favarsokii, A. P. : 
Calculating the Ionization Instability in Low-Temperature Magnetized 
Plasma. Sov. Phys. DOKLADY, vol., 14, no. 1, (1969), pp. 68-71. 

Lengyel, L. L.: On the Numerical Simulation of Ionization Instability. 

Preprint, Institut fur Plasmaphysik, Garching bei Munchen, F. R. of Germany. 

Barium Clouds: 

1. Haerendel, G. ; Lust, R. ; Rieger, E.; Volk, H. : Highly Irregular Artifi- 

cial Plasma Clouds in the Auroral Zone. Max Planck Institut fur Physik 
und Astrophysik. Muenchen. MPI-PAE/Extraterr 21, January 1969. 

2. Simon, A. : Instability of a Partially Ionized Plasma in Crossed Electric 

and Magnetic Fields. Phys. of Fluids, vol. 6, no. 1, pp. 382-388 (1963). 

3. Linson, L. M. and Workman, J. B. : Formation of Striations in- Ionospheric 

Plasma Clouds. Bull. Am. Phys. Soc. Abstract 7B7, vol. 14, no. 11, Nov. 1969. 

4. Private information to the author by David Adamson and Leo D. Staton of 

the Interplanetary and Magnetospheric Physics Section at NASA Langley 
Research Center. 



115 


5. Krall , N. A.: Drift Waves. Advances in Plasma Physics, vol. 1, 

pp. 153-199. 

6. McCormac, B. M. , Edit.: Aurora and Airglow, Rheinhold Publ. Co., 1966. 

Walt, M. Edit.: Auroral Phenomena. Stanford Univ. Press, Stanford, 
California 1965. 

Chamberlain, Joseph W. : Physics of the Aurora and Airglow. 

Academic Press, 1961. 



ii6 


GRAVITOHYDROMAGNETIC INSTABILITY IN CQNTRASTREAMING FLUIDS 

S . P . Talwar* 

NASA, Langley Research Center 
Hampton, Va, 

The hydroma gnetic stability of a multicomponent system composed of 
infinitely conducting gravitating fluids is investigated by the normal 
mode analysis. The effects of hydromagnetic and gravitational coupling 
between the components are studied in an initially static and a contra- 
streaming configuration. It is found that an N-component static nongravi- 
tating system allows, in general, (-N+1) mixed magnetosound modes in addition 
to a purely Alfven mode. The gravitational coupling results in an insta- 
bility leading to an incipient fragmentation, whose size depends profoundly 
on nonthermalization in the composite medium. Finally, the conditions of 
contras treaming instability in plasma-plasma and neutral gas-plasma 
gravitating streams are derived in view of their likely importance to star 
formation from primordial matter. 


* 

NRC-NASA Resident Research Associate 
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TRANSVERSE INSTABILITIES IN SELE-FOCUSING STREAMS 


by 


Willard H. Bennett 
North. Carolina State University 

Broadly stated, the objective in this program has been to investi- 
gate the instabilities in magnetically self -focusing streams in which 
the mean free paths of the stream electrons are more than several orders 
of magnitude longer than the stream diameter and magneto-hydrodynamic 
methods of analysis cannot be applied. Electron beams have been pro- 
duced which are electrically pinched in order to insure complete 
neutralization of the beam electron space charge, and the current and 
voltage were then to be increased until the conditions for the magnetic 
pinch were satisfied, whereupon the instabilities could be studied. 

A directional probe which has been developed recently a this Laboratory 
has been used for measuring the radial potential distribution in and 
around these pinches. 

Using the value of energy of transverse motion, B, as measured 
in high voltage beams, the current required for the- magnetic pinch to 
become dominant appears to be about the same as the maximum current 
which has been obtainable in these experiments with our present equip- 
ment. There have been several instances when it appeared that the beam 
pinched to a smaller radius than could be expected with an electric 
pinch alone. However , in order to make an adequate study of these 
beams and their instabilities, it will be necessary to make measure- 
ments at current values up to much more than the minimum critical 
pinch current. 
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CASCADE MECHANISM OF NONLINEAR INTERACTIONS BETWEEN 
MODES IN A TURBULENT PLASMA 


C. M. Tchen 

The City College of the City University of New York 


By means of a hydrodynamic model, valid for a plasma of hot elect- 
rons and cold ions, we have established the development of spectra of kinetic 
and electrostatic energies. The nonlinear interactions have been described by 
two transfer functions valid for the nonlinear transfer of modes across each 
spectrum, and by a production function serving the coupling between the two 
spectra. The closure problem, as characteristic to any nonlinear system, has 
of course to be solved, and in the present theory, the approximation is resort- 
ed to the cascade concept. The concept helped in the derivation of the equa- 
tions for the energy spectra, and in the formulation of a "turbulent relaxa- 
tion frequency". The latter frequency is very important in determining the 
structure of turbulent transport coefficients entering in all the above non- 
linear transport phenomena. Finally the equations of energy spectra are 
solved for the collisional and collisionless cases. Some of the theoretical 
results have been confirmed by experiments in gas plasmas. The spectra of 
kinetic and electrostatic energies are investigated for both the collisional 
and collisionless cases. In the former case, the spectrum of potential energy 
from the electrostatic potential fluctuations, as a function of the wave number 
k, and denoted by G(k) follows the laws G'-'k"'^^ and k“^/2 i n the inertial and 
dissipative subranges respectively. In the latter case, the dissipative law 
G' f '- / k“^ is derived. The solutions of the kinetic energy spectrum are also 
found. The anomalous diffusion, called the Bohm diffusion, is derived in the 
present theory, and is found to play a role of a sink of energy flux across a 
potential spectrum. 


£3 
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MEASUREMENTS OF THE COHERENT OSCILLATIONS AND TURBULENCE 
IN PLASMA USING RESISTIVELY- AND CAPAC 1TI VELY-COUPLED PROBES 
by John S. Serafini, NASA Lewis Research Center 

The knowledge of the diffusive and other transport properties of a plasma 
depends upon the experimental study of plasma turbulence and oscilla- 
tions. Fluctuating plasma parameters which need to be measured include 
the particle densities-} the potentials and/or fields. The writer's 
previous work (Refs. 1 and 2) used resistively-coup led electrostatic 
(Langmuir) probes to sense the plasma fluctuations. These fluctuations 
exhibited both coherent-1 i ke and random-like properties. The writer's 
present work ‘is a continuation of these studies which additionally 
employs an electrostatic probe capacitively coupled to the plasma (such 
as one described by J. A. Schmidt in Ref. 3). The fluctuations sensed 
by this type of probe are considered to be fluctuations in potential. 

The measurements have been obtained ?n a crossed-f ield plasma source 
(axial electric and radial magnetic fields) with a slightly ionized, 
low current, low-pressure argon discharge employing a hot cathode. The 
plasma had an average electron temperature of about 10 ev and an aver- 
age electron number density of about 5.0 X 10^0 cm“3 . 

Amplitude spectra from .1 kilohertz to 6 megahertz as well as the auto- 
correlations which have been obtained allow comparison of the density 
and potential fluctuations. The spectra obtained with both types of 
probes contain a peak indicative of the coherent oscillation and are 
otherwise generally similar. At the higher frequencies the spectra 
obtained with the two types of probes do not agree in slope. Possible 
reasons for the observed similarities and dissimilarities in the spec- 
tra are discussed. 

Previous work on the space-time correlations found that the coherent 
part of the current-density fluctuations were wel 1 -correlated over the 
measured range. The present work on correlations provides a like result 
for the potential fluctuations sensed by capacitively-coupled probes. 

In addition correlations of the potential fluctuations sensed by capac- 
itively-coupled probes with the current-density fluctuations sensed with 
a resist! vely-coup led probe indicate strong correlations near the fre- 
quency of the coherent oscillation. The effect of the turbulent-like 
part of the fluctuations is also observed to be present in all of the 
correlation results. It appears as minor peak in the curves of the 
correlation coefficient versus time-delay at values of time-delay near 
zero. 

Also presented are spectra of measurements of the fluctuations from two 
capacitively-coupled probes being sensed differentially. These diffei 
ential-mode spectra are compared to the spectra of the fluctuations from 
the individual probes. From these results is deduced that the diffei 
ential-mode spectra depend (I) on whether the fluctuation is due to a 
coherent oscillation or is turbulent-like in nature at the portion of 
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the spectrum -of interest and (2) on the temporal and spatial scales of 
the fluctuations. Such information prescribes necessary conditions if 
fluctuating electric fields are to be directly measured by use of such 
probe systems. 
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INVESTIGATION OF MICRO-ELECTRIC FIELDS IN PLASMA 


Ward Halverson 

NASA Electronics Research Center 
Cambridge, Massachusetts 


The interaction of the small scale electric fields with the 
charged particles of a plasma determines the basic physical properties 
which distinguish a plasma from an. un— ionized gas. .The micro-electric 
fields arise from CouLomb potentials around the plasma particles, 
density fluctuations (both thermodynamic and turbulent ) , and electro- 
magnetic radiation. This interaction can he studied by injecting 
charged particles into a plasma and observing the dynamics of these 
"test 11 particles. .The behavior of the test particles is representative 
of that of the "field" particles of the plasma and can be used for 
the computation of plasma transport properties, spectroscopic data, 
as well as the trajectories of energetic particles in geophysical or 
astrophysical plasmas. Measurements of test particle interactions 
with a plasma essentially provide a direct experimental test of 
plasma kinetic theory which has been formulated for many years and 
has been generally accepted without direct proof. 

An earlier experiment- 1 - indicated that the energy loss of 5 to 10 
KeV proton test particles losing energy to a fully ionized lithium 
plasma is less than predicted by 4o to 50$. Unfortunately, tbe 
experimental error was of the same order as the energy discrepancy, 
and a more precise measurement is necessary for a good comparison. 

An experimental program is now underway at EEC which will measure 
the energy loss of fast protons injected into a lithium plasma of 
known properties . The plasma is formed by a magnetically confined 
reflex arc with electron density of 5 x lO^- 2 cm~3 and electron 
temperature of about 1.5 eV. The test protons, with energies from 
1 to 20 KeV, interact mainly with the plasma electrons and lose 
energy at the rate of a few eV per cm of plasma traversed. The main 
experimental problem is the very small energy loss of the test particles 
in the plasma. 

A method has been devised to increase the effective path length 
of the test particles in the plasma. (See Figure 1.) Fast neutral 
hydrogen atoms are passed through the magnetically confined lithium 
arc and a fraction of the neutrals are ionized by charge exchange 
collisions in the plasma. The resultant protons are captured by the 
magnetic field and lose energy in the part of the gyro-orbits which 
is in the plasma column. The proton orbits drift tangentially around 
the arc column because of the radial gradient of the confining 
magnetic field, and the average proton energy will be related to the 
angular position of the orbit gyrocenters. The captured protons are 
lost mainly by charge exchange in the plasma, and an energy analysis 
of the particles lost from a narrow angular region of the arc column 
will show a series of peaks, (See Figure 2.) The energy difference 
between two consecutive peaks of the energy spectrum will represent 
the energy loss during a complete precession of the captured proton 
gyrocenters around the arc axis. The peak widths are related to 
the energy dispersion of the test particles in the plasma. 
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At present, all of the major experimental components have been 
constructed and are undergoing tests in the laboratory. A first 
test with the lithium arc, neutral particle accelerator, and neutral 
particle energy analyzer operating simultaneously was successfully 
conducted recently. Measurements of test particle energy loss will 
begin as soon as the microwave interferometer bridge for plasma 
density measurement and a system for measuring the captured proton 
precession rate are put into operation. 


1. W. Halverson, Energy Loss of Past Test Ions in a Plasma. Bull. 
Am. Phys. Soc. 13 , 318 ( 1968 ) and Association Euratom - C.E.A. 
EUR-CEA-FC-1+72 (1968). 
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INFLUENCE OF CONTROLLED TURBULENCE 
ON GASEOUS DISCHARGES 

by 

S. T. Demetriades, G. S. Argyropoulos, and G, Fonda-Bonardi 
STD Research Corporation, Pasadena, California 


ABSTRACT 


Tlie properties and behavior of an electrical discharge perpendicular 
to a plasma stream can be influenced by introducing gasdynamic turbulence 
of controlled properties in the flow. This study has shown analytically 
(Ref. 1) that fine-grained, intense turbulence can prevent the formation of 
arc filaments, and has verified experimentally that such turbulence can 
increase the conductance of the discharge channel. Four areas are currently 
under investigation: 

1. Measurement of the increase of cross -arc conductance S as a 
function of turbulence parameters under different operating conditions. 

The scale and intensity of introduced turbulence are being controlled by 
means of different wire screens placed in the plasma stream. Preliminary 
results shown on Fig. 1 illustrate the substantial increase of conductance 
that is obtained when turbulence is introduced. 

2. Measurement of the- relative erosion rate of the electrodes. 

This investigation was started very recently and only incomplete results 
are available. It appears that turbulence promotes a much more uniform 
erosion of the electrode (cathode) surface, and that the arc is not attached 
to a small area, but is either diffuse ox moves very rapidly and evenly 
over the entire surface of the electrode. Anode erosion is negligible with 
or without turbulence in the stream. The net amount of material removed 
from the electrodes appears larger in the presence of turbulence. This 
may be due to higher surface temperatures because of enhanced heat 
transfer to the electrodes through the turbulent boundary layer, since in 
the few tests made the flow of coolant water was maintained constant. 

Also, this effect depends on the operating conditions, in particular the 
chamber pressure and the plasma enthalpy and velocity, 

3. Measurement of turbulence properties (in particular, turbulent 
cell dimensions and density fluctuations) and correlation to the observed 
conductance enhancement. Turbulence properties are determined by 
means of a laser interferometer which measures the average (integrated) 
refractive index along a thin cylinder transverse to the plasma flow. A 
special arrangement permits simultaneous duplication of the measurement 
along two parallel paths which can be displaced with respect to each other 
along the axis of the plasma jet, from zero to 3 cm; this is equivalent to 
two identical, parallel interferometers each having an independent detector. 
Analysis of the detected signals permits computation of the decorrelation 
intervals of the turbulence parameters both as functions of time interval 

at a fixed distance, or as functions of distance at a fixed time interval. 

Figure 2 shows a schematic diagram of the double interferometer. 
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4. Analytical and theoretical studies relating to the interpretation 
of the experimental results and the prediction of model performance. The 
theory of the phenomenon and a general dispersion criterion have been 
developed (Ref. 1): Prevention and/or dispersion of arc filaments by tur- 
bulence is accomplished when the sufficient condition T >2 > 1 is satisfied, 
where D£ depends on the plasma properties, the operating conditions, 
and the turbulence scale and intensity. Computation of under the 
operating conditions of the experiments is obtained by utilizing computer 
programs that have been developed at STD Research Corporation in the 
past for the detailed analysis and prediction of plasma properties in MHD 
flows (Refs. 2 and 3). Results of these computations are shown on 
Pigs. 3 and 4 for two different levels of turbulence intensity. These 
results indicate that turbulence is more effective in dispersing arc fila- 
ments at lower arc pressures and electron temperatures and are in excellent 
agreement ■with the experimental measurements. 
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Figure 2. Layout schematic diagram of laser 

interferometer, with two parallel paths 
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ATTACHMENT OF ARTIFICIAL PLASMA CLOUDS TO THE 
MAGNETIC FIELD LINES OF THE EARTH 

Leo D. Staton 

NASA, Langley Research Center 
Hampton, Va. 


The release of an artificially produced., ionized barium, plasma cloud 
into the earth's magnetosphere initiates several transitory physical processes 
of both theoretical and practical Interest. The initial phases Include the 
ionization of the initially deposited neutral barium atoms and the expansion 
of the resulting plasma blob in a spherically symmetrical fashion until a 
size is reached at which there occurs a rough balance between the gas pressure 
within the plasma and the effective pressure of the magnetic field. During 
this period and for some time subsequent to it the plasma cloud as a whole 
is moving with respect to both the magnetic field and the ambient magneto- 
spheric plasma, because of the initial velocity imparted to the cloud by the 
space vehicle which released it. It is only after the cloud has come to rest 
with respect to the magneto spheric plasma that optical tracking methods can 
be used to measure the cloud velocity and from this velocity to allow the 
deduction of the ambient, steady, electric-field component perpendicular to 
the magnetic field. 

The time required for the synchronization of the transverse cloud velocity 
with that of the ambient plasma (i.e. the time required for "attachment" of 
the cloud to the field lines) is of practical interest because of its import- 
ance in designing such an experiment. The physical processes responsible for 
this synchronization are of theoretical interest for both fundamental plasma 
physics and magnetospheric physics. Theoretical effort to date in describing 
these processes has taken three approaches. Common to all three is the idea 
that the relative motions among the magnetic field lines, the plasma cloud, 
and the ambient magnetospheric plasma are capable of launching disturbances 
which may be propagated along and guided by the magnetic field lines. 

The first approach recognizes that the polarization electric field which 
is induced in the cloud by its relative motion with respect to the B field is 
capable of driving electric currents through that electric circuit composed 
of the magnetic lines passing through the opposite (charged) faces of the 
cloud and of a highly conducting path connecting these same lines within the 
ionosphere. These currents are carried down one of these lines and up the 
other by electrons of the ambient magnetospheric plasma, which have very 
high mobility along the lines. In this way the kinetic energy of the rela- 
tive cloud motion is dissipated in ionospheric heating. A central feature 
of the model is that the "transmission line" for these currents is filled 
with the magnetospheric plasma, which has an effective dielectric constant 
high enough to reduce the electromagnetic propagation velocity to the level 
of the Alfven velocity. 
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A second, more rigorous approach is to seek solutions of the MHD 
equation for a related, mathematically tractable, problem in an attempt to 
gage the validity of the model described above. It is possible to- obtain 
a rigorous, one-dimensional formulation of an MHD problem in which an 
infinitely long strip-like current source in the x-y plane (with the 
current flow across the strip) excites a pure MHD wave propagating along 
the z-axis. Unfortunately, the one-dimensional character of the solutions 
is destroyed if the source strip is truncated to finite size. This fact 
implies that a correction must be added to the first model above to account 
for the inevitable divergence of some energy away from the direction of the 
transmission line. 

The third approach, one which was originally proposed by workers at 
the Max Planck Institute, ignores the polarization electric field and instead 
regards the magnetic field lines as frozen into both the magnetospheric and 
cloud plasmas. As the cloud moves the lines are then distorted, resulting 
in a magnetic disturbance which propagates down the field lines with Alfven 
velocity. The magnetospheric plasma is then disturbed with the result t ha t 
the cloud relative momentum becomes shared with the particles of the magneto- 
spheric plasma. 

Continued work is required since, at present, theoretical values of the 
synchronization, or attachment, time for the computing models are very nearly 
equal, allowing no certain description of the physical processes involved. 
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COMPUTER SIMULATION OF EXPANSION OF IONIZED BARIUM CLOUD 
PERPENDICULAR TO THE EARTH'S MAGNETIC FIELD LINES 


F. C. Grant 

NASA, Langley Research Center 
Hampton, Va. 


The elongated form of the Ionized barium cloud suggests that the motion 
transverse to the Earth's magnetic field be modelled by the two-dimensional 
motion of a large number of infinitely long rods of charge with appropriate 
masses. A computer program based on the work of Hohl and Hockney has been 
prepared to simulate the cloud's transverse motion. A more precise treatment 
of the dynamics has been introduced so that the treatment of the electronic 
rods exactly parallels that of the ionic rods. 

Large amounts of computing time are required for the simulation, so a 
prime concern is to use the smallest number of rods consistent with good 
simulation. Another primary concern is to use the proper initial conditions. 
These two matters are subjects of parametric study. 

When the cloud is subdivided into a sufficiently large number of small 
charged masses and they are properly started, the subsequent transverse 
motion provides a measure (when compared to observed cloud motion) of the 
exciting magnetic field fluctuations, or turbulence, of the Earth's magnetic 
rield . The effect of magnetic field fluctuations on the transverse diffusion 
may he estimated theoretically. 
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RADIATION FROM AN ELECTRIC DIPOLE SOURCE IN 
ANISOTROPIC COMPRESSIBLE MEDIA 

A 

A. K. Sundaram 

NASA, Langley Research Center 
Hampton, Va. 


A detailed discussion on the radiation characteristics of an electric 
dipole embedded in an infinite homogeneous, anisotropic, compressible 
plasma is made. The cases treated include the effects of Hall and electron 
inertia terms in the generalized Ohm's law and finite pressure in the momen- 
tum equation. The investigation bridges the gap between the MHD and the 
conventional magnetoionic theories on radiation from electromagnetic sources. 
Mathematically, the far fields are obtained using a simple procedure out- 
lined in the analysis which deviates from the existing theory, 1 and it is 
shown that the method has some advantage in obtaining the radiation resistance. 
From physical point of view, the directional distribution of waves excited 
from the dipole source is discussed when Hall and electron inertia effects 
are included and the results thus derived hold good for all frequencies, 
particularly for frequencies comparable to ion cyclotron frequency. The 
validity of the theqry is also examined in detail for a special case when 
E4-VxB=0, as usually treated in the literature. 
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SIMULATION OF SOLAR WIND-EARTH MAGNETOSPHERE INTERACTION 
by D- L. Chubb, NASA Lewis Research Center 

A laboratory simulation experiment of the interaction of the solar wind 
and the earth's magnetic field is being conducted at NASA-Lewis. The 
plasma flow is produced by an MPD arc and a magnetic dipole is used to 
simulate the earth's magnetic field. 

Preliminary results were encouraging. Regions, of intense light were 
observed in locations similar to those of the earth's radiation belts. 

A short movie will be presented of these preliminary results. Although 
these results are encouraging significant improvements are necessary to 
obtain a proper simulation of the solar wind-earth magnetosphere inter- 
action. Density and velocity requirements necessary to obtain a proper 
simulation will be discussed. These requirements result from the follow 
i ng simulation criteria: (1) bow shock thickness less than shock stand- 

off distance; (2) stagnation magnetic field twice field for undisturbed 
dipole; (3) collisional mean free paths greater than dimensions of exper 
ment; and ( 4 ) hypersonic flow. 

In addition to the simulation criteria, diagnostic results of a colli- 
sionless stagnation flow will be presented. These diagnostics consist 
of plasma potential and magnetic field measurements. Finally, modifica- 
tions that are being made to produce a proper simulation of the whole 
solar wind-earth interaction will be discussed. 
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CURRENT AND PLANNED RESEARCH IN ARGON ION LASERS 
IN THE PLASMA PHYSICS AND GAS LASER BRANCH 

P. Brockman 

NASA, Langley Research Center 
Hampton, Va. 


The argon ion laser offers the possibility of continuous single 
mode operation at high power in the visible and medium power in the UV. 
There are a wide range of applications including underwater communica- 
tions, materials testing, and deep-space communications. Continuous 
argon ion lasers have been operated at output powers over 100 watts and 
at efficiencies up to jv 0.2 percent. Tube diameters are typically a few 
millimeters but an increase in laser power with tube diameters up to 
15 mm has been shown. 


The interest in argon ion lasers at Langley evolved from plasma accel- 
erator experiments involving various high-power nonequilibrium gas discharges 
at low pressures. To date the program has consisted of equipping a laboratory, 
designing and building several experiments, and beginning cooperative research 
with Old Dominion University. 


One experiment has consisted of looking for gain in the plume of the 
continuous MPD arc. At about the time these measurements were begun gain 
in a pulsed MPD arc was reported. ^ To date, both double pass balanced gain 
measurements (Fig. 1-A) and internal cavity measurements (Fig. 1-B) using a 
pulsed TRW laser have indicated less than measurable gains ( <_ 1%), across 
the plume. It is still possible that inversion does exist in the rotating 
spoke which is evident in the MPD arc for many operating conditions . Future 
experiments will consist of looking for gain in the spoke using a d.c. laser 
as a probe and of .operating spoke free at much. higher powers using the milli- 
second pulse power supply described at this meeting. ^ 


Additional experiments have consisted of constructing and operating 2.7 
and 7 millimeter internal diameter quartz lasers. These are static filled 
lasers in which the fill pressure can be varied. The 7 mm i.d. laser has 
been operated at a range of pressures and currents from 15 to 30 amperes. 

For these conditions performance is similar to that reported in reference 1. 

Quasi-steady operation at power inputs up to 1 MW will be studied using 
flat pulses of 1 raillisec. or longer, and the limitations due to radiation 
trapping will be investigated. In this manner the problems for high-power 
steady lasers will be studied without involving materials problems. Efforts 
will be made to obtain proper quasi-steady conditions through simulation of 
steady-state conditions of the lasing plasma and of the tube walls. For 
example, preheating the quasi-steady pulsed laser with a steady d.c. discharge 
will be used to investigate the effect of gas temperature on radiation trapping. 
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An additional laser is under construction. This is a hollow-cathode, 
hollow-anode design with metal constricting apertures,^ (Fig. 2). In this 
design studies of high-power operation will be made and investigation of 
potential drop along the axis as well as a study of applying varying 
potential to the apertures. For operation at high input powers large- 
diameter hollow cathodes which have been operated for steady operation of 
the MPD arc at currents up' to 900 amps will be used. 


The work under the Old 'Dominion University grant will be discussed in 
another talk. 
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SPECTROSCOPIC MEASUREMENT IN AN MPD ARC* 

N. M. Nerheim + 

Jet Propulsion Laboratory 
Pasadena, California 


ABSTRACT 


Spectroscopic studies of the emission spectra of the plasma produced by an 
MPD arc have revealed a highly non-equilibrium distribution of the argon ion 
excited states. For some operating conditions, the spectra suggested the 
possibility that population inversions may occur for some of the excited 
levels. Furthermore, the spontaneous emission of an argon ion laser^) was 
noted to be very similar to the spectra' of the MPD arc. Because of the, simi- 
larities, an investigation of the conditions necessary to produce population 
inversions in the MPD exhaust was started. The beam of a static argon ion 
laser was directed through the flowing plasma and monitored to detect any 
changes in the beam intensity. An enhancement of the beam would indicate 
stimulated emission and therefore an inversion of the upper and lower states 
of the laser line. The instrumentation was sufficiently sensitive to detect 
a 1% change in the beam intensity. The probing laser beam was directed 
through the MPD exhaust at right angles to the flow direction, and a multiple - 
pass mirror system was used to increase the path length in the plasma to about 
500 cm. The exhaust plume was probed as a function of position over a wide 
range of arc operating conditions, exhaust pressure, and beam strength. 

No evidence of stimulated emission was found. The method used was limited 
to the investigation of the normal A+ lasing lines that constitute the out- 
put of the probing beam, the strongest of which is the 4880 & line. Because 
the lower levels of these lines are the resonant levels, direct spectro- 
scopic measurements of the ratio of the upper and lower levels of the transi- 
tion could not be made. This ratio can be measured for transitions originating 
from the higher energy levels and the study revealed the possibility of popu- 
lation inversions for some of these lines. However, this speculation is 
strongly influenced by the accuracy of the transition probabilities used to 
calculate the population densities from the measured line intensities. Be- 
cause of the lack of experimental values for the lines originating from the 
higher energy level, the transition probabilities of a number of lines with 
4d upper levels were determined. The measurements were made by H. N. Olsen 
using techniques that were very similar to those used previously by 01 sen 
to determine argon A-values. The radiation source for these studies was the 
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plasma produced by a free-burning 400 ampere DC arc in argon at 1.1 atm. A 
plane -grating monochromator was used to isolate the individual line and fused 
silca optical components were used to allow observation of lines with wave- 
lengths in the near uv region. Side-on observations were corrected for the 
background continuum and the measured integrated intensities were inverted 
using a standard Abel inversion technique. The A-values obtained were lower 
than the theoretical values by nearly a factor of two, but agreed well with 
one set of experimental values. A comparison of the published sets of Aril 
transition probabilities show that there are significant differences in the 
available values. The comparison is given in terms of the temperature differ 
ences that will result from the use of different sets of A-values when they 
are used to determine plasma temperatures. Additional measurements of transi 
tion probabilities are planned. These measurements will include argon atom 
lines and possibly some xenon lines. 


(1) R. I. Rudko, C. L. Tang, J. Appl. Phys., 38, 4731 (1967). 

(2) H. N. Olsen,' JQSRT, 3, 305 (1963). 
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ULTRAVIOLET SPECTRUM OE LOW PRESSURE ARGON ION DISCHARGES 

G. Ofelt, J. Becher and E. Young 
Old Dominion University, Norfolk, Va. 


The objective of the study is to understand the basic excitation 
mechanisms of the argon ion laser in order to be able to increase the 
power output and efficiency. This study is being conducted in close 
cooperation with the Plasma Physica and Gas Laser Branch of the Aero- 
Physics Division at the Langley Research Center and is intended to be a 
joint 'effort utilizing the plasma technique experience of that group to- 
gether with the spectroscopic experience of the Optics and Spectroscopic 
Group at Old Dominion University. 

Previous investigations have been concerned with visible and near UV 
spectroscopic observations. However, we intend to examine resonance 
transitions for singly ionized argon which lie in the vacuum UV region at 
approximately 1000 A. Since radiation trapping is believed to be a major 
limitation of argon ion lasers, it seems important to examine spectral 
line intensities in this region which will lead to a -determination of popu- 
lations of the various states of the ion. This study was undertaken to 
improve the technology of gas lasers and diagnostics and to examine the 
lasing system with the aid of spectroscopic techniques in the UV region. 

It is hoped that information will also be obtained for establishing popula- 
tion inversions for multiply-ionized argon which may cause lasing action 
for wave lengths as low as 1850 &. Parametric dependence of the operating 
conditions will be investigated for creating population inversion among the 
various lasing levels. 

Laser configurations to be investigated are both continuous and pulsed 
time durations from microseconds to milliseconds. Variation of pressure, 
current density, and capillary diameter are to he made in the course of this 
investigation. 

We have conducted a preliminary investigation on a small diameter 
discharge (approximately 4 mm) for the purpose of calibrating our spectro- 
scopic equipment for the argon discharge application. Initial studies were 
conducted at relatively high pressures of approximately 1 Torr, as compared 
to the approximately 0.1 Torr laser working pressures at this capillary 
diameter. The intent was to obtain as strong as possible line intensities 
during the initial investigation to facilitate the initial orientation. The 
experimental configuration used for this preliminary study of the argon ion 
discharge consist of an approximately 5 cm long capillary discharge tube, 

4 mm in diameter, which was mounted to a 1/2 meter Seya Monochromator. Obser- 
vations were made from the visible region of the spectrum to approximately 
300 £ with the aid of a dosium salicitate scintilator in front of the photo- 
multiplier tube. The output signal was analyzed by a PAR box car integrator 
which enables improved signal- to-noise ratio when the electronics is operated 
in the pulsed mode operation. The discharge power supply is capable of 
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supplying up to 1000 amps at 1000 V, for a period of 50 secs. Examina- 
tions of argon ion discharges made at this laboratory show a production 
of a continuum from approximately 1000 A to the near UV. . This is the 
molecular transition which is caused by the close proximity of argon atoms 
while undergoing a pulsed excitation. Descrete lines are superimposed on 
this continuum. At present, we are attempting to establish the experimental 
parameters that influence the production of the continuum and their relation 
to the lasing lines. In addition, hydrogen is present as an impurity which 
should aid in the determination of the gas temperature. 



FEASIBILITY OF A MAGNETOGASDYNAMIC CROSS-FIELD LASER* 

G. R. Russell 

Jet Propulsion Laboratory 
Pasadena, California 

ABSTRACT 


An analytic study has been carried out to determine if population inversions 
of atomic excited states can be attained in the entrance and exit regions of 
a cross-field accelerator using argon as a working fluid. A computer program 
is used in the analysis which includes the simultaneous effects of electron 
heating, ionization and recombination rates, and the rates of population of 
atomic excited states. The generalized equations for a three component fluid 
with different specie temperatures are formulated using the theory of 
Zhdanov (Ref. 1), and .include a generalized Ohm's law and electron energy eq- 
uation. Electron thermal conductivity, thermal diffusion, and radiative 
transport are included. Optically thin free-free and free-bound radiation 
losses are calculated using the theory of Biberman, et al (Ref, 2); also 
included is the nonequilibrium thin line radiation resulting from transitions 
between excited states. The resonance transitions result in a leakage loss 
from the thin line wings plus the effect of diffusion of photons in the thick 
line core. 

A ten level hydrogen-like model is used to calculate the collisional-radiative 
recombination and ionization rates, and the populations of the first ten 
transition arrays of atomic argon. It is not assumed either for the case of 
ionization or recombination that the rates of change of excited state popula- 
tions tend to zero as is done by Bates and Kingston (Ref. 3) , Radiative 
transition probabilities are calculated using the theory of Bates and 
Damgaard (Ref. 4) . Collisional excitation and deexcitation cross sections 
are evaluated using the theory of Gryzinski (Ref. 5) for transitions permit- 
ted by the selection rules for dipole radiation. 

Solutions have been obtained for constant pressure accelerator flows and for 
flows accelerated in divergent ducts. Depending on the initial conditions 
of the flow prior to acceleration, population inversions are predicted for 
both the entrance and exit regions for atomic transitions in the infrared. 

At the accelerator entrance the sudden increase in electron temperature leads 
to an inversion of the 5p»3d transition array. The predicted gain for this 
transition is relatively small and would not allow oscillation in a cavity 
less than about one meter in length. In the exit region of the accelerator 
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where the electron temperature is rapidly decaying, electron capture associated 
with collisional-radiative recombination causes inversions of the 4d~5p tran- 
sition array. Depending on the final values of the electron temperature and 
density, the gain for this transition array is large enough to make it pos- 
sible to attain oscillation in a cavity of the order of 10 cm in length. 

On the basis of these predictions, a cross field laser has been constructed 
at JPL to attempt to produce lasing in a cavity placed downstream of the ac- 
celerator exit. In this device argon is pre -heated in four MPD arcs, ex- 
hausted into a water-cooled plenum chamber, and then expanded through an 
uncooled graphite Mach 3 nozzle. The gas then passes through a constant 
pressure cross -field accelerator containing twenty-one uncooled thoriated 
tungsten electrode pairs. Provision is made to vary the current density and 
potential along the accelerator axis to minimize the flow of Hall currents,. 

The laser cavity is located at a fixed point in a large vacuum tank that 
houses the entire experiment. The cavity is positioned relative to the ac- 
celerator exhaust by moving the accelerator and associated equipment, includ- 
ing the electromagnet, on a track within the vacuum tank during a test. 

Design of the complete system including the laser cavity and associated 
optics has recently been completed. Preliminary diagnostic studies of tbe 
gas properties in the accelerator have begun prior to attempting to provide 
the correct threshold excited state populations required to achieve oscil- 
lation in the laser cavity. 


(1) Zhdanov, V. M. (PMM Vol. 26, Mo. 2 , 1962, pp 280-288), pp 401-413, 
Applied Mathematics and Mechanics (Translation) . 

(2) Biberman, L. M. , Morman, G. E., and Ulyonov, K. M. , Opt, and Spect. 

10, pp 297-299, 1961. 

(3) Bates, D. R. and Kingston, A. E., Plant Space Sci., 1963, Vol. 11, 

pp 1-22. 

(4) Bates, D. R. and Damgaard, A., Phil. Trans. Roy. Soc. (London) 242A, 

14 (1949). 

(5) Gryzinski, M. , Phys. Rev., Vol. 138, No. 2A, pp 336-358, 19 April 1965. 
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MOLECULAR LASERS USING ELECTRICAL DISCHARGES AND HIGH 
FLOW RATES. CURRENT AND PLANNED RESEARCH IN THE 
PLASMA PHYSICS AND GAS LASER BRANCH 

F. Allario, R. Lucht and R. V. Hess 
NASA, Langley Research Center 
Hampton, Va. 


Within the past few years, considerable research on gasdynamic 
molecular lasers using mixtures of CO^ with various gases has been per- 
formed with the aim of obtaining high laser power. This system although 
attractive because of its large input capability, has proven to have 
extremely low efficiency. Recently, interest in other systems capable of 
producing high laser power has arisen because of their relatively high 
efficiency, although limited power inputs. Among these systems are^ ^ 3 
electrical discharge lasers^ogerating in relatively high flow rates ’ * 
and chemical laser systems. ’ In both systems, rapid progress is being 
made in increasing laser power per unit length and there is promise of 
even higher power input and higher efficiencies than presently available. 

In the last six months, members of this Branch in collaboration with 
'Prof. Hassan of North Carolina State University, have developed plans for 
research in molecular lasers using electrical discharges and high flow 
rates of CO^, and mixtures. Some preliminary studies have already 
been performed. 

Professor gassan will discuss aspects of the theoretical program at 
this conference and in this paper some preliminary experimental results 
and plans for future work will be discussed. In figure 1, a schematic of 
the apparatus used to study electrical discharge lasers with high flow rates 
is shown. As is indicated, a simple electrode configuration and gas injec- 
tion scheme were used, but results obtained indicated that an improved version 
would probably be more desirable for applying greater power to the electrical 
discharge. Furthermore, a convenient optical cavity was chosen but again 
calculations show that the choice of an improved optical cavity could lead to 
a substantial increase in efficiency of the system. Volume flow rates in 
the one- inch-diameter tube were varied between 5 and 150 liters/min, the 
upper limit being set by available flow meters. Electrical power was supplied 
by a 100 kw (10 kv, 10 A) d.c. power supply. The laser beam showed a multi- 
mode structure and power was recorded with a conventional power meter. A 
gain-measuring system has also been assembled with an auxiliary laser to 
measure laser gain parallel and perpendicular to the flow axis. 

A parametric study has been performed of the variation of laser output 
power with flow rates of the constituents (C0„ , Ng* H ) and has been re- 
lated to the voltage-current characteristics 01 the discharge. These prelimi- 
nary results show that almost an order of magnitude increase in power output 
can be expected by varying the flow rate considerably less than 2 orders of 
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magnitude. For example, laser power was increased from less than 10 
watts (efficiency ^ 5%) to approximately 60 watts (efficiency a. 10%) 
for a tube length of less than half a meter (18") by increasing the 
flow rate from 5 to 130 liters 1/min. These results by far exceed the 
maximum power outputs per meter for no- flow lasers and are consistent 
with recently published results in references 1 and 2, showing that 
convective cooling of the gas mixture by the high flow rate helps to 
maintain a low gas temperature leading to enhanced laser power. The 
comparatively low efficiencies obtained (v 5% and 10%) suggests that 
improvement in the optical cavity and flow uniformity are needed; at 
the lower flow rates cooling of the laser tube walls is desirable. 

However, the more pertinent result for the present high flow rate experi- 
ment is that it was possible to increase input power to the discharge by 
increasing the flow rate, which for a given efficiency of’ course corres- 
ponds to an increase in laser power. The inability in the present experi- 
ment to obtain further increase in laser power arises from the fact that 
the input power was limited by the flow rates and the difficulties in ob- 
taining higher current densities due to the transition from glow discharge 
to arc-discharge operation which prevents the attainment of further useful 
power for lasing. 

It was further observed in our experiment that the discharge appeared 
to strike in an annulus between the electrodes rather than distributed 
over the cross section. Thus, optimum use is not made of convective 
cooling and the current density is not uniformly spread over the cross 
section. An improved electrode, configuration and injection system should 
permit higher power inputs in the present, experimental setup and should 
also permit operation at higher pressures without transition from glow- to 
arc-discharge yielding higher laser powers. 


A research program is planned to investigate the operation of molecular 
lasers with electrical discharges and high flow rates, at comparatively high 
pressures (> 50 torr) and large tube diameters. Problems associated with 
cross-flow lasers^ 5 ^ will be investigated using various techniques including 
magnetic fields to prevent bowing of the discharge with the flow. Operation 
of glow discharges with flow rates at high pressures involves new research 
problems concerning the transition of the (abnormal) glow to arc discharge. 
Diagnostics of the distribution of electric field, particle density and 
electron temperature are planned, together with transverse gain measurements 
to relate the distribution of lasing action to the distribution of plasma 
properties in the discharge. Transverse gain measurements along the tube 
length should yield information on the efficiency of multiple excitation 
of N£ along the discharge (and cavity) length, and should allow a com- 
parison of long tube lasers and short transverse flow lasers. Fluid mixing 
schemes will also he studied. 


High power pulsed .operation of electrical discharges at high pressures 
and high flow rates will be investigated to determine the limits of current 
density and power to which the discharge laser can be extended since in 



155 


pulsed operation, heating of the gas and filling of the lower laser level 
can be kept small through the use of short pulses. ^ Diagnostics on pulsed 
discharges are also planned utilizing the powerful techniques of time 
sampling and signal averaging. 

A program is also in prog'ress to study the feasibility of excitation 
mechanisms other than the glow discharge for molecular and chemical lasers 
in cooperation with Prof. Richard Schneider from the University of Florida. 
As a preliminary experiment, the effect on a conventional 00^ laser of the 
high energy products released in the reaction of H 3 with thermal neutrons 
will be investigated. Experiments at higher particle densities than in 
conventional tube lasers are- ultimately planned for investigation. The 
effect of neutrons on HF chemical lasers using UF, as fluoride source will 
also be studied. l ' 
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THE ROLE OF ELECTRICAL DISCHARGES IN POPULATION 
INVERSION IN C0 2 FLOW LASERS* 

H. A. Hassan 

North Carolina State University, Raleigh, North Carolina 27607 


The objective of this research is to analyze CO^ laser systems in -which the 
population inversion is obtained by electrical discharges . This entails solving 
the governing equations for various laser parameters, such as mass flow rate, 
length and diameter of the tube laser , pres sure and composition. In such discharges 
^excitation of CO molecules to the upper laser level results mainly from energy 
transfer from vibrationally excited N 2 while the lower laser lever Is depopulated 
mainly by collisions with He. 

If one assumes that molecules in different quantum states are different species, 
then the governing equations are the conservation of species equations , the over- 
all momentum and energy equations, and an electron energy equation. These equations, 
together with the reactions given in Eq. (.5) below, will serve as the basis of our 
model. The considerations presented here are based on these equations. 

Examination of the conservation of species equation, which can be written as 

PY. 

P -p + pu • VI. + V . (PY.V.) = W. , Y. = Pi /p (1) 


where p is the density, u the mean velocity, V. and W. are the diffusion velocity 
and mass rate of production of species i, shows that £here are three characteristics 
times in a steady laser system. These are 


t 

o 




( 2 ) 


where L and d are the length and diameter of the tube, respectively, and U^, p^, 

D and W are characteristic velocity, density, diffusion coefficient and production 
rate, respectively For a successful laser operation 

t < t (3) 

r o 


For a given temperature, t is inversely proportional to the pressure; thus, for 
an increase in pressure an N 2 molecule can be excited and de-excited many times 
for a given tube length. When t > t the laser is of the gas-transport type; in 
this case the cooling is mainly due to the convective cooling m C AT. This can 
be seen from 


t . a 2 u ,, -lie 

a _ _ o L p 

t D L ir e EL 

o o 


(U) 


* 
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where is the Lewis number (of order unity) , m is the mass flow rate* C is 
the specific heat at constant pressure and K is the thermal conductivity.-^ The 
m C 


parameter 


K L 


can he thought of as the ratio of convective cooling to cooling 


hy conduction and diffusion through the wall. Thus to improve the laser power 


one needs to operate with high 


m C 

I 

K L 


For a given laser system t is obtained from consideration of the rates of 
vibrational energy transfer in ?he system. Based on the recent survey of Taylor 
and Bitterman 2 , the important relaxation precesses in a CO^ laser system are 


N 2 ( V = 0) + e + N 2 (V= 1) + e 


(i) 


C0 2 (o,! 1 ^) + M -+ C0 2 (0,0°,0) + M + 667 C m 1 


(ii) 


C0 Q (0,0° ,1) + N. (V = 0} + CO (0,0° ,0 ) + U (Y = l) + 18 C 

d d d d m 


-1 

(iii) 


C0 2 (0,0 ,0) + M £ C0 2 ( 0,m ,0) + M + \ r&ns _ (iv) 


(5) 


where M stands for any species in the system: C0^ 1? 2 , H , H 2 0, CO, etc. Re- 
action (ii) describes the exchange of energy between translational and vibrational 
modes while reactions (iii) and (iv)describe the vibration - vibration energy ex- 
change. Examination of the rates indicates that H is active in depopulating 
the lower laser, level and has little effect on the upper laser level. On the 
other hand is active in depopulating all laser levels. This indicates that 
there is an optimum H 2 0 concentration beyond which a loss in laser power will 
take place. 

The role of H in a C0 2 laser system is discussed next. It has been sug- 
gested that, in addition to depopulating the lower laser level, it has a high 
thermal conductivity; this results in a decrease in the gas temperature and this, 
in turn, increases the population inversion. However, increased conductivity re- 
sults in an increase in the wall losses and this results in a decreased- efficiency. 
Thus, the explanation for the improved efficiency, especially in the gas— transport 
lasers, must lie somewhere else. It is proposed here that the improved efficiency 
is a result of the high electron temperature in the discharge. This can be seen 
from the electron energy equation 


J -E = 
e 


3 k N m 
e e 


(T - T) S 
e 


6 v 
s es 

m 


( 6 ) 


where J is the electron current, E the electric field, T the electron temperature, 
T the gas temperature ¥ the electron number density, and m the particle mass, 6 
the energy-loss factor Ind v is the collision frequency o? species s. The 
quantity <S is unity for monatomic gases and large for heavy particles (about 18 for 
IT )3, thus S to maintain a high electron temperature one needs to employ an inert 
monatomic carrier gas. 
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PUMPING MECHANISM OP C0 2 LASER AND FORMATION 

JL 

RATE OE C0 2 PROM CO AND o" 

C. J. Chen + 

Jet Propulsion Laboratory 
Pasadena, California 


ABSTRACT 


The pumping mechanism of a high-current pulsed pure C0 2 laser has been inves- 
tigated. It was found that there is a time delay of the laser pulse behind 
the current pulse. For a fixed gas pressure, there is an upper bound dis- 
charge current and pulse width below which the time delay decreases as the 
discharge current is increased. Above the upper bound current, the time 
delay is independent of the discharge current and the pulse width. The delay 
time decreases very rapidly as the gas pressure is increased. The gas pres^ 
sure p and the delay time t exhibit the inversed cubic relation as t a l/p . 
The intensities of the atomic oxygen line 7771 S. and the band head of the 0-1 
band of CO 4835 8. increase as the magnitude and/or pulse width of the dis- 
charge current is increased. These intensities reach a plateau for certain 
upper bound values of the discharge current and current pulse width. After 
the current pulse is over, the line intensities dec^g gradually to zero. 

The electron density measured typically is about 10 cm and decreases to 
about 17o of this initial value when the laser pulse appears. The electron 
temperature drops from about 2ev to the gas temperature within SOjxs. The 
gas temperature, which varies from room temperature to 400°K by heating the 
gas with a temperature tape, was measured with a thermocouple in the center 
of the tube. 

These experimental results indicate that the pumping mechanism of the upper 
level of the C0 2 laser is due to the reaction CO + 0 -* C0 2 . During the 
pulse discharge, the C02 is decomposed into CO and 0. The saturation pheno- 
mena as shown in delay time as a function of discharge current and pulse width, 
and the plateau phenomena, as shown in line and band head radiation indicate 
the percentage of CO 2 , which decomposed into CO and 0 depends upon the dis- 
charge power input which is proportional to the magnitude of the current and 
the pulse width. For a sufficiently high power input, the C0 2 is totally 
decomposed into CO and 0. The saturation and plateau phenomena appear when 
CO 2 is totally decomposed. Aftgr the current pulse is over, CO and 0 start 
to recombine to form excited C0 2 molecules. Each recombination contributes 
one excited C0 2 molecule. Thus, the recombination rate of CO and 0 is equal 


*This work represents the results of one phase of research carried out in the 
Propulsion Research and Advanced Concepts Section of the Jet Propulsion 
Laboratory, California Institute of Technology, under Contract NAS7-100, 
sponsored by the National Aeronautics and Space Administration. 
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to the rate of pumping of the upper laser level. The delay time of the laser 
pulse behind the current pulse is the time required to build up the laser up- 
per state population to meet the threshold condition for oscillation. By 
equating the threshold population density of the upper laser level to the 
number of CO 2 molecules formed, the rate constant K can be determined: 

2.88 x 10 13 (kT) 3 6 -1 

K = — i — — cm sec 

pt 

The value of K can be obtained by measuring T, p.and t. The result for a 
pressure range from 4 to 10 mm Hg at a temperature of 300°K is 

K = 8.5 x 10 33 cm s ^ 

The temperature dependence of .the value of K is obtained by measuring the 
time delay at different gas temperatures which is accomplished by using a 
heating tape wrapped around the discharge tube and the tubing upstream of the 
flowing gas system. The temperature dependence is assumed to be in the fol- 
lowing form: 

_ B 

T 

K = A e 

Evaluation of the constants A and B yields the complete expression for the 
reaction constant : 

1670 

rr o in in’ 32 T 6-1 

K = 2.19 x 10 e cm s 

It should also be noted that the water vapor impurity in the tube, which is 
very difficult to entirely remove even with elaborate precautions, does not 
affect the measured reaction rate in the present case. Since the recombina- 
tion energy in the reaction with OH is carried away by the product H, the CO 2 
molecule formed by this process is in the ground state (both electronic and 
vibrational) . Therefore, there is no contribution to the population of the 
upper energy level of the laser. 
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CONTINUOUS -WAVE CHEMICAL LASER OPERATION 

WITHOUT SUSTAINING EXTERNAL ENERGY SOURCES* 

Terrill A. Cool 
Laboratory of Plasma Studies,, 

Cornell University, Ithaca, N. Y. 14-850 


C ont inuous -ware laser operation at 10. 6ji has been achieved in the 
BF-COg and HP-CO^ molecular systems by purely chemical means. Both lasers 

operated solely by the simple mixing of- bottled gases; no external sources 
of energy were required. 

A chain reaction mechanism has been found -which provides continuous 
pumping of the upper CO^ laser state by vibrational energy released from. 

chemical reaction. The reactions can be sustained- continuously without 
an external source of energy. The chemical reactions are: 

+ NO ^ ONE + E 

F + U 2 > (DF)* + D 

D + F 2 ^ (DF)* + F (l) 

The analogous set of reactions obtained by substituting hydrogen in 
place of deuterium has also been employed successfully. The vibrational- 
rotational energy of the vibrationally excited (DF)* or (EE)* product 
molecules of the above reactions is coupled to the upper COg laser level 
by means of intermolecular vibrational energy transfer processes. 

Experimental results concerning the performance capabilities of these 
first purely chemical lasers will be presented. These measurements show 
that the DF-C0 2 system is capable of high optical gain and high power output. 
More significantly, this laser offers a remarkably efficient means for 

the continuous direct conversion of chemical energy into laser output. 

/ 

The DF-C0 2 purely chemical laser offers unique possibilities for 
efficient self-contained energy storage and utilisation. 


^Supported by NASA under contract NGL 33-010-064- and by ONR under 
contract N 00014~67-A-0077-0006. 
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EXPERIMENTAL AMD THEORETICAL STUDY OF 
CHEMICAL INVERSION IN GAS LASERS 

Andrew R. Saunders 
NASA, Langley Research Center 
Hampton, Va. 


Production of critical inversion densities by chemical reactions in 
gas lasers has assumed increasing importance in recent years. Pulsed and 
e.w. gas lasers using chemical inversion were operated successfully quite 
recently uji^izing vibrationally excited hydrogen halides formed in chemical 
reactions. ’ Highest power levels were obtained in a continuous flow 
reactor where the vibrational energy of the hydro, gen halide was used to excite 
the 00°1 level in C0„ to obtain laser emission on the 00°1 -s- 10°0-5- and 
00°1 ■+ 02 2 0 transitions. ^ 

f j tfc 

In our recent studies we have investigated chemical inversion in CO. 
formed in low pressure CO - 0 2 flames at 50 - 80 torr. The radiative 2 
transitions for CO^ are shown in figure 1. It was known from earlier work~* 
that electronically or vibrationally excited CO^ can form in the reactions 

CO + 0 ->■ co 2 ( 3 B 2 ) 

co 2 ( 3 B 2 ) + co 2 ( l B 2 ) 

co 2 ( B 2 ) -+ co 2 (v 1# v 2 , v 3 ) 

Radiative transitions from th^ excited electronic level *B 2 take place 
into vibrational levels of the ground state of CQ„. By comparison of 

s 

experimental measurements of radiative emissions corresponding to the 
11 0 02 z 0, 10°0 01 1 0 and 00°1 -»■ 00°0 transitions with the theoretical 

values for equilibrium vibrational distribution as shown in figure 2 for the 
latter two, it was shown that for nearly equilibrium population of the 10°0 
symmetric stretching and 020 bending modes the population of the 00°1 assym- 
metric stretching mode is substantially less than the equilibrium value. 

In addition the data show a definite correlation between the electronic 
transitions populating the vibrational levels and the infrared emission inten- 
sities for these levels. Thie preferential filling of the bending and sym- 
metric stretching vibrational energy levels of CO. may be explained by the 
transition from the bent ^B 2 state to the linear configuration. In 

this transition the bending mode vibration would be preferentially excited 
and energy readily transferred into the symmetric stretching mode by Fermi 
resonance coupling whereas the uncoupled assymmetric stretching mode would be 
populated largely by collisional transfer. 

Comparing these results obtained by chemical excitation on CO„ with 
the results of electron impact excitation in a C0 2 electric discharge laser 
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cavity the role of electron impact excitation as the primary process in 
producing 00 °1 inversion becomes quite apparent. In the recent studies 
of Clough and Thrush^ on the formation of vibrationally excited N 2 0 in 
the reaction N + N0 2 -> ^0 (v^, v^) +0 at about 0.1 torr preferential 

excitation of the v^ and V£ modes was observed similarly to our results with 
C0_. Since ^0 is isoelectronic with C0„ and the distribution of the 
vibrational energy levels is very similar, these results strongly support our 
observations on the CO - 0 - CO^ system. 

At this point a substantial experimental and theoretical effort is 
required to evaluate the electronic, vibrational and rotational energy trans- 
fer processes controlling chemical inversion. Quantum mechanical calculations 
are required to predict the probabilities of electronic and vibrational 
excitation in a given reaction and the initial energy distribution of the 
species formed in the reaction. Using the rate parameters available for given 
systems we shall attempt to relate the internal electronic, vibrational and 
rotational energy transfer processes in the presence of chemical reactions to 
the inversion densities obtained in a laser cavity for the same system. Com- 
plexity of the chemical rate processes is indicated by the rate processes on 
CC >2 inversion as shown in figure 3. 

The initial effort will be concentrated on measurements of the rotational 
and vibrational distributions of the laser emission radiation in halogen-halide 
and other chemical lasers in pulsed and c.w. operation using high flow rates as 
a function of the reaction kinetics. The theoretical studies will be concen- 
trated on solution of the radiation coupled chemical kinetic rate equations 
with detailed description of the fluid dynamic and chemical rate processes. 
Mechanisms leading to higher inversion densities and total radiated power will 
he examined both experimentally and theoretically. 

References 

1. Cornell, Paul H. ; and Pimentel, George C. : Hydrogen-Chlorine Explosion 

Laser. II. DC1. J. Chem. Phys. , vol. 49, no. 3, Aug. 1968, pp. 1379-1386. 

2. Gross, R. W. F. : Chemically Pumpsed C0„ Laser. J. Chem. Phys. vol. 50, 

1969, pp. 1889-1890. 

3. Cool, Terrill A., Falk, Theodore J., and Stephens, Ronald R. : DF-CO^ and 

HF-CC^ Continuous-Wave Chemical Laser, Applied Physics Letters 1970. 

4. Koopmann, R. K. ; and Saunders, A. R. : Non-Equilibrium Vibrational Energy 

Distribution of CO^ in CO-O^ Flames at Reduced Pressure. Jour. Quant. 
Spectr. and Radiative Transfer. 1970. 

5. Callomon, J. H. ; and Gilby, A. C.: The Carbon Monoxide Flame Spectrum Under 
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Thermal Laser Excitation by Mixing in a Highly Convective Flow 

Wayne G-. Burwell 

United Aircraft Research Laboratories 

CW infrared laser emission has been achieved from the mixed flow region 
formed by injection of CO2 or U2O into a vibrationally excited N2 stream, formed 
via the rapid adiabatic expansion of thermally excited (1000-3000 °k) if 2 through 
a supersonic nozzle. Population inversion and laser emission have been achieved 
on the 001 100 transition in COg, brought about by near-resonant vibrational 

energy transfer during collisions with excited Eg. Results are presented showing 
the i nf luence of flow properties downstream, of the expansion nozzle on optical 
gain and laser power characteristics . An analysis of nonequilibrium flow/ optical 
field interactions is presented. Effects of penetration and mixing are represented 
as instantaneous changes in free-stream composition, specific enthalpy and Mach 
number. Subsequent vibrational and radiative energy transfer processes are coupled 
with fluid flow conservation equations. Qualitative agreement between the thermal 
mixing laser analysis and the experimental results has been found. 
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INVESTIGATIONS ON A PLASMA-EOCUS APPARATUS 


N. W. Jalufka and J. H. Lee 
NASA, Langley Research Center 
Hampton, Va. 


Investigations are under way to learn more about the plasma in a 
plasma-focus apparatus. This is a relatively newly-discovered method 
of producing extremely hot and very dense plasmas, of the order of 50 
million degrees Kelvin, 10^ particles per cubic centimeter, and a 
million atmospheres pressure. The method is of much interest in the 
scientific .community and of considerable potential in several fields of 
both research and practical application. 

Some of the results obtained to date are as follows: 

1. Neutron time-of-f light measurements have been successfully used to 
obtain the neutron energy distribution - with, however, lower resolution 
than is desired but much better than that obtained by others. The statis- 
tics is currently being improved by using large (15") scintillators. 

2. Some of the distributions had a rounded peak (like a Gaussian), 
indicating thermalization , but some had a sharp-pointed peak, indicating 
a beam model, or pure radial collapse. 

3. The width of the neutron spectrum was measured and compared with 
predictions of theoretical models of the plasma focus. It was found that 
the half width of neutron spectra is much greater than expected from thermal- 
ized ions at temperature of 5 keV. Therefore ei-ther the ions are not 
thermalized or T ^10 keV. 

4. The electron temperature was also measured with soft x-rays 

(bremsstrahlung) by using differential absorbers made of aluminum and 
copper. The temperature was found to be in excess of 5 keV. 

5. Other investigators had found soft x-rays, of order 10 keV. Some 

had found up to 300 keV. We were the first to report 500 keV and greater. 

6. The x-rays are from a point source rather than a rod-shaped or 

moving lines, giving sharp pictures of x-rays penetrating more than 1 inch 
thick aluminum sheet , or equivalent . 

7. We have pinhole-camera x-ray pictures that show the shape of the 
plasma near the central electrode. 

The big unsolved problem at the time we started this experiment and 
still at the present time is whether or not the plasma is thermalized. We 
have recently exposed nuclear-emulsion stacks to measure the x-ray energies 
or spectrum in various directions. These have been sent to Vanderbilt and 
Brookhaven for analysis and will be used to derive the electron velocity 
distribution. 


*Vanderbilt University 
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With nuclear emulsion, refined measurements are planned for Cl) 
angular distribution of hard x-rays, (2) calibration of emulsions with 
radio isotopes and an electron accelerator, and (3) electron and proton 
spectra with a magnetic spectrometer. 

A preliminary spectroscopic investigation in the region between 500 2 
and 1400 2 has been carried out. A one meter, normal incidence, vacuum 
spectrograph was used to record the spectra photographically. 

The (time integrated) spectra in the region above 1000 2 was dominated 
by copper line radiation (CuII,CuIII) while the spectra below 1000 2 was 
dominated by line radiation from ionized oxygen (OIII, OIV, OV, 0VI) . 

Spectra were also recorded with a center electrode cap of aluminum and 
carbon. These spectra were predominated by line radiation from highly 
ionized stages of the electrode material. Some copper lines were still 
observed but were very weak indicating that some electrode material is 
carried by the current sheet. The use of the different electrode materials 
had no apparent effect on the production of hard x-rays. 

An investigation of the radiation in the soft x-ray region is planned. 
Spectra of the radiation in this region will be, obtained with a bent crystal 
spectrograph. Line radiation in this region is due to transitions in highly 
ionized atoms and it will be possible therefore to identify and classify 
these lines. An attempt will be made to determine the shape and peak of 
the bremsstrahlung continuum. The peak of this continuum radiation is a 
strong function of temperature (electron) and these measurements should give 
some indication of the temperature obtained in the plasma focus. 

An investigation of the effects of different gases will be made and 
spectra of the highly ionized atoms of these gases will be recorded and 
classified. Also the effects of different electrode materials will be 
investigated further. 
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ATTENUATION OF X-RAYS FROM PLASMA FOCUS 
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AN INVESTIGATION OF HARD X-RAYS FROM A PLASMA FOCUS 
J. H. Lee, D. S. Loebbaka, and C. E. Roos 
Vanderbilt University, Nashville, Tennessee 
Presented at the Sixth NASA Intercenter and Contractor Conference 
on Plasma Physics at Langley Research Center, 

December 8 , 9 , and 10 , 19^9 


INTRODUCTION 

A dense plasma focus produced by fast collapse of current sheets 
at the end of a coaxial electrode is a copious source of neutrons (from d-d 
or d-t reaction) and x-rays. Despite extensive studies of this device by 
many workers, the mechanism of neutron and x-ray productions is not well 
understood. A few models have been proposed but none of which is 
universally accepted. . The boiler model is based on the assumption that 
ions and electrons are thermalized in the plasma focus at a kT of a few 
keV and anisotropy of the neutron yield and spectra are results of plasma 
CM motion. The. beam-target model considers a high current of charged 
particles accelerated in an electric or magnetic field giving the production 
mechanism as well as the anisotropy. Accurate spectral distributions of 
neutrons and x-rays would distinguish between these models. With the plasma 
focus apparatus at the Langley Research Center a maximum photon energy of 
greater than 500 keV has been observed by means of differential absorber 
method.-^- However, the determination of the spectral distribution of x-rays 
was difficult since the absorption coefficients for most elements are near 
a minimum for the energy range in consideration. This investigation is an 
effort to obtain the spectrum of hard x-rays from the plasma focus using 
electron sensitive nuclear emulsions to clarify corresponding electron distri- 
bution in the plasma supplementing other diagnostics currently used to 
determine the nature of a plasma focus at the Center. This report covers the 
first phase of the work initiated in February, 19^9 • 

EXPERIMENT 

The plasma focus apparatus used in this experiment is described in 
a separate presentation of this conference. Hard x-rays are emitted simul- 
taneously with neutrons at the maximum compression of the plasma. However, 
it is found that no correlation exists between the intensities of neutrons 
an d x-rays , excluding the possibility of significant neutron induced gamma- 
ray emission. To determine the adequate dose of exposures a trial set of 
13 emulsion (UQO micron thick Ilford K— 5) stacks were placed around the 
source at varying distances from 15 to 90 cm. end exposed to 100 runs of 
plasma focus formation. The second series of experiments were carried out 
to determine the spectral distribution of electron tracks produced in the 
emulsion by either photoelectric or Compton effect of x-rays from the plasma 
focus. For this run 600 micron thick Ilford emulsions obtained by a new 
order were used. The emulsions, stacked behind a 1.5 mm thick lead absorber, 
were placed at a distance of 90 cm from the source on the axis of the 
electrode. Discharges in helium were also used for comparison with the 
deuterium plasma. The electron ranges were measured under a Leitz scanning 
microscope equipped with an 18 micron field of view reticle grid. 
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RESULTS AND DISCUSSION 

The histogram in Figure 1. is a result of electron range measurements 
in 3 plates with 3807 tracks (longer than 18 microns) observed over background. 
The ranges were converted to electron energies using the relation given by 
Barkas.^ A few tracks extended more than 600 microns indicating electron 
energies of greater than 500 'keV. In the figure the dotted curves represent 
the expected electron track distributions, arbitrarily normalized, for three 
different models of the source. The curves were obtained by taking account 
of the attenuation through the lead absorber and the efficiency of nuclear 
emulsion. A kT = 10 keV is chosen for the curve for Maxwellian electron distri- 
bution since it is maximum temperature one could reasonably assign for the 
plasma focus . It is clear that a Maxwellian distribution alone cannot account 
for the high energy spectrum observed. The curve for thick target brems strahlung 
is drawn for an electron beam with- kinetic energy of 500 keV. This model, by 
itself, cannot fit the observed spectrum. The power law curve is derived from 
the electron distribution, f(E)cc E“I and the best fit value of y = k ± 1. 

(The spectrum of thick target bremsstrahlung corresponds roughly to y = 0.5- ) 

The present results are in agreement with thei power law distribution reported 
for the energy range up to iO keV by Meskan, et al. although the accelerating 
mechanism (turbulent heating) is yet questionable. 3 It is interesting to note 
that such' distribution of electrons is proposed for astrophysical x-ray sources , 
such as SCO-X-l.^ After determination of response function of emulsion to 
jnonoenergetic ga mm a-rays the experimental results will be analyzed in detail 
and comparison experiment with a theta pinch will be carried out. Beside 
the scanning of the electron tracks, the recoil protons produced by neutrons 
were also studied. The results show that the maximum neutron energy reaches 
3 MeV, confirming the measurement made by the time-of-f light method, and 
anisotropy of neutron yields in. axial and radial direction is 1.6 ± 0.^, which 
is higher than the expected value (l.25) for the boiler model. 

CONCLUSIONS 

There is a strong Indication that the hard x-ray spectrum observed 
with nuclear emulsions are produced by a power law distribution of electrons 
in the plasma focus. The detailed study of this dependence should give 
decisive information on the proposed heating mechanisms in the plasma. The 
authors are grateful for the support and encouragement on this investigation 
by G. P. Wood (NASA), J. P. Barach (VU), and colleagues associated with the 
plasma focus apparatus at the Langley Research Center. We also wish to 
express our thanks to Mrs. T. Davis, Mr. S. Thatphithakkul and Mrs. U. Ahuja 
for their excellent scanning efforts. 
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TIME- INTEGRATED NEUTRON FLUX MEASUREMENTS ON A 
PLASMA FOCUS BY GAMMA SPECTROMETRY 

* 

L. P. Shomo and K. H. Kim 
NASA, Langley Research Center 
Hampton, Va. 


The plasma focus is produced by the discharge of a 20 kilovolt - 
25 kilojoule capacitor bank between the electrodes of a coaxial plasma 
gun filled with deuterium under pressures of a few Torr. The 1 mega- 
ampere discharge yields up to 10 1 neutrons in a burst of approximately 
100 nanoseconds half-width. The extremely short duration of the neutron 
production time prevents direct measurement of the neutron flux. Instead, 
the decay products of a silver target activated by the neutron burst from 
the plasma focus are counted. Figure 1 shows the experimental apparatus. 

A brief summary of the sequence of operation is: (1) A burst of 2.45 MeV 

neutrons is produced ; (2) some of the neutrons enter the polyethylene 
moderator and are thermalized; (3) the silver becoyg^ activated by the 
thermalized neutrons through the two processes Ag (n, y) Agl^® and 
Ag^- 99 (n, y) AgHO; (4) the activated silver decays by beta and gamma 
emission; (5) the gammas are detected by a three by three inch Nal (T&) 
detector; (6) the multichannel pulse-height analyzer provides the energy 
resolved gamma intensity spectrum, and (7) the intensity of the 656 keV 
gamma peak of Ag - ^ with 24 second half-life is measured. 

The use of a discrete gamma peak eliminates some of the corrections 
and uncertainties inherent in the beta counting method. The dead- time 
correction due to the 125 microsecond characteristic dead-time of a four- 
geiger-tube beta counting system can be neglected due to tbe 1 microsecond 
dead-time characteristic of Nal (T&) detectors. The beta counting method 
nondis criminat'ely counts betas from decaying Ag-^ and Ag^ having half- 
lives of 24 seconds and 2.3 minutes, respectively. This introduces errors 
due to the uncertainties of two neutron capture cross sections , and compli- 
cates the calibration procedure due to the two different half-lives. The 
greater penetrating power of the 656 keV gammas allows the use of a thick 
silver target thereby eliminating any disadvantage suffered by the gamma 
counting method due to the relative weakness of the 656 keV gamma peak. 

Thus the use of the 656 keV gamma peak of Ag-^® provides a significant in- 
crease in accuracy over the normally used beta counting technique. Yields 
of 10 9 to 10l° neutrons per burst from the plasma focus device have been 
measured by this method with an estimated uncertainty of +16 percent, which 
includes the +5 percent error contributed by the calibration source. 

This improved method is presently being applied to neutron flux aniso- 
tropy measurements on the plasma focus device. It is hoped that these im- 
proved measurements will shed some light on the validity of either the 
boiler or target models of the plasma focus mechanism. 

* 

North Carolina College, Durham, N. C. 
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FIGURE 1. 
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COMPUTER SIMULATION OF PLASMAS 
Velvin R. Watson 

NASA Ames Research Center, Moffett Field, California 


Recent advances in computers and cathode ray tube displays 
have made the computer an attractive tool for solving some problems 
in plasma dynamics. A recent investigation at Ames Research Center 
has shown that the computer is particularly useful for studying the 
dynamics of a eollisionless plasma using a computer simulation 
technique. An application of this technique is described in this 
paper, and several problems that can be investigated using this 
technique are listed. The further application of the computer to _ 
provide insight into the dynamics of a collisional plasma is 
currently being investigated. 

The computer simulation technique was used to study the plasma 
dynamics within a very high specific impulse accelerator. This 
simulation was accomplished by representing the many ions and 
electrons in the plasma by several thousand representative particles. 
These several thousand representative particles were assumed to 
constitute a good statistical sample of the actual ions and electrons 
in the plasmas. Starting with an initial distribution, the simul- 
taneous motion of the representative particles was calculated as a 
function of time. First the electric and magnetic fields were 
calculated from the distribution of particles; next the motion of the 
individual particles during a small increment in time was calculated 
using these fields with the equations of motion; and then the change 
in the distribution of particles caused by the motion of the particles 
during the small increment of time was calculated. This procedure was 
repeated for each new time increment and the field and particle 
positions were recorded for each time increment to provide a time 
history of the self-consistent fields and particle motion. 

The motions of the particles, the electric potential, and the 
magnetic field were displayed simultaneously on a cathode ray tube. 
These displays, which were recorded on motion picture film-, illustrated 
the coupling between the particle motions and the fields. In partic- 
ular, the displays illustrated how the plasma motion modified the 
applied electric potential, how the currents induced large magnetic 
fields, how the ions received their energy directly from the electric 
field, and how the ions were deflected toward the accelerator exhaust 
by the induced magnetic field. 

This computer simulation technique should be useful for solving a 
number of problems involving collisionless plasmas, including the 



following: (l) the flow of a collisionless plasma around "bodies, 

(2) the flow of plasma within the n Mather Focus" at the end of the 
cathode of a coaxial plasma gun, (3) propagation of waves in 
collisionless plasmas — particularly in inhomogeneous plasmas . 

An investigation of the techniques for applying the computer 
to solve problems in collisionless plasma dynamics indicated that 
one should not always proceed as far as possible with analytical 
mathematics prior to formulating the problem for the computer. 

The computer simulation technique described above employs the 
equation of motion for each of the several thousand particles. 

Instead of solving these several thousand equations one could have 
solved the Vlasov equation (which has been obtained from the 
particle equations of motion by further use of analytical mathe- 
matics) by a finite difference technique. Both techniques were 
used to solve for the voltage-current characteristics of a one- 
dimensional diode and the solutions were compared with experimental 
measurements of voltage and current. The computer simulation 
technique produced a much more accurate representation of the 
experimental measurements than the solution of the Vlasov equation 
by a finite-difference technique. 

The results of the computer simulation are the visual displays 
of the plasma motion rather than a closed form mathematical 
expression for some continuous function, such as a distribution 
function. The changes in the dynamics due to changes in parameters 
are therefore observed visually rather than interpreted from an 
analytical, expression for the distribution function. In a sense, 
the computer simulation is similar to a laboratory experiment where 
one changes the parameters and observes the results. Several major 
advantages of the "computer experiment" over the laboratory experi- 
ments are (l) measurements of the plasma properties can be made 
without perturbing the plasma with probes, (2) a microscopic as well 
as a macroscopic picture of the plasma is obtained, and (3) the 
plasma temperatures are not limited by availability of power supplies 
or by the technology of cooling at the plasma boundaries . 

The author believes there is a great potential for the "computer 
experiment" as a tool for solving problems in plasma dynamics and 
recommends that the computer facilities for performing these "computer 
experiments" be expanded to be commensurate with this potential. 


REPORTS 

1. Watson, Velvin R. : "A Collisionless Model for High Specific 

impulse Accelerators," AIAA 7th Electric Propulsion Conference, 
Williamsburg, Virginia, March 3-5, 19^9- 

2. Watson, Velvin R. : "Computer Simulation of a Plasma Accelerator, 

Stanford University Institute for Plasma Research Report No. 313, 
May 1969. 
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PULSED COAXIAL PLASMA GUP 
William C. A. Carlson 

NASA Ames Research Center, Moffett Field, California 


Experiments have heen conducted on a pulsed coaxial plasma gun 
powered with electrical energy from a 20 KV capacitor bank. 
Experiments to date have been conducted with between 3 and 20 KJ of 
electrical energy. The discharge is initiated by introducing gas 
into the evacuated, electrically energized gun. This mode of 
operation results in a distributed region of current flow which 
remains relatively fixed in space, rather than in a localized 
traveling current sheet which is typical for most similar guns. 

This later mode of operation has been described as a "snow plow 
model . " 

Several techniques have been used to determine the plasma 
velocity. These include streak photography, doppler shifting of 
radiation from the line, and impact pressure on a target. These 
experiments indicate that the velocity produced by the gun operating 
in the distributed or' deflagration mode is approximately an order of 
magnitude greater than that obtained when operation is in the 
traveling current sheet or detonation mode. 





NEW SPECTRAL LINES OF HIGHLY-IONIZED NEON 
FROM THE LANGLEY THETA PINCH 
Hans Hermans dorfer, Langley Research Center 


To help meet the 'Videspread and urgent demand for basic spectroscopic 
data (wavelengths, energy levels, transition probabilities, ionization 
potentials, and electronic binding energies) from a wide variety of fields 
of science and technology — plasma physics, solid-state physics, astro- 
physics, thermodynamics, spectrochemistry, laser research," the theta-pinch 
facility has been engaged in a program to produce and interpret new spectro- 
scopic data. The method and the results to date, including identification 
of a score of new neon lines, are described. 

In addition, the theta-pinch personnel have achieved the following, 
which are listed here but are not discussed in the paper: 

Put the new grazing- incidence UV spectrometer (10 <. X S. 1000 X) 
into proper operation 

Experimented with the ratio of working gas to carrier gas, e.g., 
NE: D , and with various fractions of the capacitor bank, to produce 
the correct temperature of the plasma and sufficient intensity of the 
spectral lines 

Received and further developed shutter from Dr. Schneider for 
controlling duration of exposure 

Attempted to produce 0^X11 Line @ 23.3 X — "satellite line due 
to inner-shell transition 

Diagnosed and eliminated cause of prefires and capacitor failures 

Published "laboratory observation of a 1 visible emission line of 
highly ionized argon at 4412 2 ." in Astrophysical Journal 

Published "Tables for the Calculation of Radial Multipole Matrix 
Elements by the Coulomb Approximation" in Astrophysical Journal 

Began effort to produce new spectra in isoelectronic sequence of 
B I to Ar XIV 


Collaborated on development of zone-plate camera with NRL for 
EUV photography of the Sun 

Designed and obtained apparatus for Thompson-scattering measurement 
of T 

e 

Designed and procured sliding-spark light source (currently used for 
U spectra at University of Florida) . 
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A METHOD FOR THE CALCULATION OF LARGE NUMBERS OF 
DIPOLE AND QUADRUPOLE TRANSITION PROBABILITIES 


L. P. Shomo and G. K. Oertel 
NASA, Langley Research Center 
Hampton, Va. 


The present status of our knowledge of transition probabilities or 
of f-values is poor compared to that of the wavelengths of atomic spectral 
lines. F-values are known only for a small number of the lines for which 
wavelengths are known, and the accuracy of the available values is usually 
poor. Although f-values as well as wavelengths are atomic constants 
obtainable in terms of integrals over wave functions, only the former involve 
the overlap between two different atomic states, so that cancellations often 
lead to radial matrix elements which are small compared to the uncertainties 
of the individual atomic wave functions. 


The calculation of f-values requires a knowledge of the wave functions 
of the two states involved in the transition. It is necessary to employ 
approximations to determine the wave functions for nonhydrogenic atoms or 
ions since Schrodinger's equation can be solved analytically only in the 
simple one-electron atom case. Several techniques that provide approximate 
wave functions have been developed. The most sophisticated of these tech- 
niques is the self-consistent-field method. The Coulomb approximation method 
is applicable to many lines of interest, and for those transitions where its 
assumptions are met, it gives f-values which are in close agreement with the 
self-consistent-field results. Even for relatively simple systems where the 
time required for self-consistent-field calculations is still relatively 
short, the Coulomb approximation method will result in a reduction of com- 
puting time by about three orders of magnitude. 


The Coulomb approximation formalism was generalized to the multipole 
case x to meet the need for knowledge of forbidden electric-multipole transition 
probabilities in astrophysics, and because large numbers of matrix elements, 
both dipole and quadrupole, are required for computations of the Stark 
broadening of isolated ion lines. ^ Tables for computing Coulomb approximation 
radial matrix elements were calculated-*- for dipole transitions s-p, p~d, d-f, 
f-g, g-h, and h-i, and for quadrupole transitions p-p, d-d, f-f, g-g, h-h, 
i-i, s-d, p— f, d-g, f-h, g-i, and h- j . A FORTRAN IV computer program was 
developed^ which takes the energy levels of a given element and stage of 
ionization and applies the appropriate selection rules and calculates dipole 
and quadrupole line strengths, transition probabilities, wavelengths, and 


* 

NASA Headquarters, Washington, D. C. 
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dipole f-values for all allowed transitions. The input is restricted to 
levels with no equivalent electrons and LS-coupling. Absorption oscillator 
strengths for Helium I calculated by this computer program are compared 
xrLth f-values calculated by other methods in Figure 1. 
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TABLE I. - ABSORPTION OSCILLATOR STRENGTHS FOR NEUTRAL HELIUM 


Transition 

Self- cons istent 
field 

Scaled Thomas 
- Fermi 

Coulomb 

approximation 

- 3 - 

Coulomb 

approximation 

4-+ 

BIB 

BB 

Sum rule 
modified 

Effective 

charge 

Wiese et. al. 
"best values 1 ' 


l'S - 2‘P 

0.2719 


0.2531 



0.3555 



0.2762 

SP 

3'P 

0.07203 


0.06825 



0.0722 



0.0734 

SP 

4'P 

0.02738 


0.02795 



0.0282 



0.0302 

W 

5'P 



0.01414 






0.0153 

LS 

6'P 

0.00793 


glJHISIIsssS 



0.0082 



0.00848 

LS 

2'S - 2'P 

0.3578 


0.342 ' 


IS 

0.3719 

■ H 

0.392 

0.370 


0.3764 

SP 

3'P ' 

0.1646 


0.242 



0.1431 


0.150 

0.156 

0.2948 

0.1514 

SP 

4'P 

0.0508 


0.0684 


0.04865 

0.04970 


0.062 

0.0503 

0.08728 

0.0507 

w 

5'P 



0.0296 


0.02218 

0.01936 

.»iB 

0.025 

0.0221* 

0.03818 

0.0221 


6'P 

0.01257* 


[0.0157_ 


0.01209 


BS 

0.012 

0.0124 

1 

0.01257 


2'P - 3'S 

0.04640 


0.0278 


0.04806 

0.04859 


0.04797 


0.02304 

0.0480 

w 

4'S 

0.00824+ 


0.00549 


0.008579 

0.008308 


0.00686 


0.004150 

0.00824 


5'S 



0.00211 


0.003213 

0.003041 




0.001519 

0.00308 

CA 

6'S 

0.00151* 


[0.00105 


0.001593 

0.001491 




0.0007502 

0.00151 


2'P - 3'D 

0.7253 


0.575 ‘ 


0.7104 

0.7006 

0.755 

0.7357 


0.6947 

0.711 

w 

4'D 

0.1205+ 


0.114 


0.1207 

0.1168 

0.118 

0.132 


0.1199 

0.1205 


5'D 

0.04308* 


0.0436 


0.04346 

0.04106 

0.0416 

0.04313 


0.04355 

0.04308 


6'D 



0.0218 


0.02105 

0.02054 

0.0199 



0.02115 

0.0213 

CA 

2 3 S - 2?P 

0.60021 

0.538 

0.5457 

0.5319 

0.542 

0.559 

0.535 

0.6686 

0.53908 

SP 

3 3 P 

0.05705 

0.0638 

0.05982 

0.06478 

0.0826 

0.052 

0.0768 

0.1876 

0.06446 

SP 

4 3 P 

0.02083 

0.0256 

0.02370 

0.02581 

0.0270 

0.031 

0.0232 

6.06502 

0.0231 

W 

5°p 


0.0124 

0.01143 

0.01022 

0.0123 

0.013 

0.0114* 

0.02999 

0.0114 


6°P 


0.00694 

0.006372 


0.00665 

0.007 

0.0061* 


0.0061 


2 s P -3 3 S 

0,07639 

0.0692 

0.06813 

0.07204 


0.07662 


0.02700 

0.693 

W 

4 3 S 

0.01159* 

0,0105 

0.01014 

0.008939 


0.00520 


0.004288 

0,01159 


5 3 S 


0.00376 

0.003581 

0-003159 




0.001488 

0,00365 

CA 

e’s 


0.00182 

0.001720 

0,001527 




0.0007089 

0.00176 

CA 

2?P - 3 a D 

0.6234 

0.611 

0.6168 

0.6089 

0.553 

0,6530 


0.6969 

0.609 

W 

4 3 D 

0.1229* 

0.123 

0.1244 

0.1192 

0.129 

0.0792 


0.1325 

0.1229 


5 3 D 

0.0466&+ 

0.0472 

0.04759 

0.04550 

0.0512 

0.04951 


0.04958 

0.04668 


6?D 


0.0236 

0.02374 


0.0260 



0.02446 

0,0215 

CA 


+■ Calculated by computer program in appendix A 

++ Calculated by linear interpolation in the orginal tobies of Bates and Damgaard (ref. 6) 

* Selected as "best value" by Wiese et. al. (ref. 9) 

sp Results of variational calculations of Schiff and Perkens (ref. 17) 

W Results of variational calculations of Weiss, A. W. presented by Wiese et. al. (ref. 9) 

LS Results of Low and Stewart quoted by Dalgarno and Stewart (ref. 19) 

CA Coulomb approximation results presented by Wiese et. al. (ref. 9) 

[ ] Results enclosed may be considered less accurate than the Coulomb approximation results (see appendix E) 
— Coefficient of fractional parentage used 
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HARTREE-FOCK CALCULATIONS OF ENERGY LEVELS AND TRANSITION 
PROBABILITIES OF ARGON I AND OF HIGHLY-IONIZED ATOMS 
IN THE BORON I ISOELECTRONIC SEQUENCE 

* 

Louis J. Shamey 
NASA, Langley Research Center 
Hampton. Va. 


Numerical Hartree-Fock calculations are being performed to obtain 
accurate nonrelativistic atomic wave functions. As an initial approxi- 
mation, the Hartree-Fock-S later (HFS) self-consistent calculations are 
performed, using the well-known programs of Herman and Skillman. This 
procedure does not distinguish between different terms within a configura- 
tion, thus resulting in only one set of energies and orbitals for a given 
configuration. These orbital energies, screening constants, and wave- 
function slopes are invaluable as input estimates for the more accurate 
(and sensitive) Hartree-Fock (HF) calculations, using the well-known programs 
of Charlotte Froese-Fischer. Because this Hartree-Fock calculation does 
distinguish between the different LS terms within a configuration, each 
term will be associated with its own set of orbital wave functions and 
energies; for this reason, the so-called average energy is no longer the 
same constant for the different term energies and thus should not he sub- 
tracted from the diagonal' electrostatic matrix elements, as is commonly done. 


These Hartree-Fock calculations are currently being applied to two 
theoretical spectroscopic problems, the work still being in progress. 

The HF results are used as the best available theoretical wavef unctions, 
from which all appropriate atomic parameters can he calculated. The^first 
application is a calculation of forbidden transitions in Argon I 3 p 4 p - 
3 p^. Electric quadrupole transitions take place bgtwe^n 3p S 0 a^d 
3p^ 4p Spin-orbit interactions couple the 3p 4p D^, and 

terms. Intermediate coupling calculations were performed to obtain tbe 
transition probabilities A(J-J’) of order 10 , which are relatively strong 
for forbidden transitions. Calculations are still in progress for obtaining 
magnetic dipole transition probabilities within this same transition array. 


The second application of these HF calculations is a systematic investi- 
gation of the spectrum, both energy levels and transition probabilities, of 
highly ionized atoms in the Boron X iso-electronic sequence, with particular 
emphasis on two atoms currently being studied ir* the laboratory: Ne VI, 

for which much is known observationally , and Argon XIV, for which little is 
known. Also of interest will be the appropriate stages of ionization for 
Mg, Al, Si, S, Ca, and Fe. The effects of configuration interaction and 
spin— orbit interaction are to be included in these calculations. This work 
is still in its early stages; preliminary results are encouraging. 


*NRC-NASA Resident Research Associate 
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LASER ENERGY ABSORPTION IN NON-UNIEOEM PLASMAS 

Loren C. Steihhauer and Harlow G. Ahlstrom 
University of Washington - Department of Aeronautics and Astronautics 


The advent of high power pulsed lasers, now makes it possible to 
heat a dense plasma to extremely high temperatures. Short, powerful 
laser pulses may heat the plasma to very high temperatures before 
dispersive mechanisms act to cool the plasma. One of the first questions 
that must be answered is how does the plasma react during the actual 
heating process, i.e. during the laser pulse. 

Under certain conditions of interest, some effects can be neglected: 
thermal conductivity, electron-ion equilibration, and plasma motion. 

This case for a general, non-uniform overdense plasma has been solved 
using the method of matched asymptotic expansions to obtain an analytic 
solution. 

Solutions have been calculated for all ranges of the absorptivity. 
These solutions show that the heating is either "simultaneous" when 
the absorptivity is small or due to a "heating wave" when the absorptivity 
is large. In the "simultaneous" case, the entire underdense region 
is heated simultaneously. A hot spot is produced near the critical 
density where the laser frequency equals the plasma frequency, and much 
of the laser radiation is reflected and lost. In the "heating wave" 
case, a heating wave propagates into the underdense region but never 
reaches the critical density. No hot spot is produced, and none of the 
laser energy is lost to reflection. 




J 0 /A ^joules /m 2 


SUMMARY OF PLASMA CONDITIONS. WHICH CAN BE ACHIEVED FOR A 
LINEAR DENSITY GRADIENT WITH A Nd + GLASS LASER 
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METHODS FOR ESTIMATING DYNAMIC MOVEMENT AND STRESS 
OP METAL PLATES WHICH CONDUCT CURRENT FROM AN 
UNDERDAMPED CAPACITOR BANK DISCHARGE 

M. D. Williams and J. A. Moore 
NASA. Langley Research Center 
Hampton, Va. 


Collector plates (metal plates with clamping through bolts) are used 
with many high current low inductance experiments. They are subjected to 
large pulse forces which must be structurally supported. The estimation 
of the movement and stress caused by these forces is a problem which occurs 
continually in the design of plasma devices. One solution to the problem 
has been to design the collector plates to support the peak undamped force 
under static conditions (the maximum conceivable load, neglecting resonances). 
The analysis and results presented here are for dynamic conditions. 

The analysis treats two theoretical models: (1) the mass-spring 
oscillator, and (2) the traveling wave model. 

The mass-spring oscillator model is valid for force pulse wavelengths 
much larger than the collector plate dimensions in the direction of the force. 
The mass of this model is driven by the force developed hy a damped ringing 
capacitor hank discharge. Mechanical damping and damping due to collector 
plate movement are neglected. The equation of motion is solved and the peak 
movement, normalized with respect to static displacement, is plotted for 
various electrical to mechanical frequency ratios, n, and bank damping ratios, 
<$, (fig. 1). In practice, this model describes collector plate movements 
when driven by capacitor banks with frequencies not greater than several kHz. 

For bank frequencies much larger than several kHz the force pulse wave- 
lengths are much smaller than the usual collector plate dimensions in the 
direction of the force and the traveling wave model is applicable. The 
equation of motion for the -traveling wave model is presented. The peak 
collector plate stress is localized in the first force pulse wavelength and 
is found to be a function of the bank damping ratio only. Peak stress normal- 
ized with respect to static stress is plotted as a function of damping ratio, 
6, (fig. 2). 

The results of this analysis exhibit several physically understandable 
phenomena and can be used to design and/or estimate the strength of collector 
plate systems. 




Log n 

Figure 1. - Normalized maximum movement of mass-spring model 
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ON THE MOTION OF VORTICITY IN MAGNETOHYDRODYNAMICS 

Robert C. Costen 
MSA, Langley Research Center 
Hampton, Va. 

A new equation has recently been derived for the motion of vorticity 
in a general MHD fluid, including the effects of viscosity, compressibility, 
nonhomogeneity, and nonconservative forces. The derivation of this equation 
was made possible by introducing a kinematic modification in the derivation 
of Kelvin's circulation theorem. In Kelvin's derivation, one obtains a 

d -> ~ 

formula for ff wn do, which is the time rate of change of the flux of 

vorticity w through a surface a of finite area, where a is taken to 
move with the fluid. If the fluid is inviscid, subject to conservative forces 
only, and if the density p is a constant or is a single valued function of 

^ _y 

the pressure p, one finds that — ff m*n do = 0 and concludes that under 

dt 

these conditions vorticity moves with the fluid. In our modified derivation, 

j 3 -v A 

one also finds an expression for — // wn do, but in this case the surface 

of integration a is taken to move - not with the fluid velocity v - but 
with" an arbitrary abstract velocity field U(x,y,z,t). The abstract velocity 

field U is then chosen so that -^7 ff orn do = 0 x?ithout imposing any 

restrictions on the fluid medium or the force fields. The velocity f^eld U so 
chosen may then be regarded as defining the motion of the vorticity w, and 
TJ is found to satisfy the equation 

(U - v) x to - - div 1? + ^ (v-v ) - - (p E + J x b) = V(J) 
p p s p e 

where p is the. fluid density, P is the pressure tensor including viscous 
terms, q is the source density of distributed fluid sources, v is the 
velocity of the source fluid, p E + J x B is the Lorentz force density, and 
4>(x,y,z,t) is an arbitrary scafer field. Note that U is not, in general, 
determined uniquely because of the arbitrary function <j>. Physically, this 
means- that for general distributions of vorticity there are arbitrarily many 
flux-preserving motions^ But in cases where a vortex tube is isolated, the 
displacement velocity U of the tube should be uniquely determined. The 
fluid source density q is included in this equation since the axis of a 
vortex is frequently ^ convenient place for injecting or withdrawing fluid. 

The pressure tensor P _^has been retained in full form to illustrate that a 
displacement velocity U may be assigned to the vorticity even in the 
presence of viscous diffusion. 

Little progress has been made to date in solving this nonlinear equation. 
But one may neglect viscosity and linearize the equation for the case of weak 
vorticity, weak fluid source density, and weak electromagnetic forces. This 
gives typically 
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X -y 

p (U - v ) X0 ) 1 + q (v - v ) VP - I x B = 0 

o o 1 \L o s p olo 


where q^, p^, and are first-order perturbation fields. Here 

p = cont. and the other zero order fluid fields satisfy the equation 
9v v o 2 

•7j“ + V r> ?= -V* where <£> o is the gravitational potential and 

<j> o is- an arbitrary function. The terms of the vortex motion equation all 

have the dimensions of force per unit volume. The first term p (U-v ) x m, 

o o 1 

represents the Magnus force density oh vortex filaments which do not move 

with, the unperturbed flow velocity v . The second term q n (v° - v S ) 

o 1 

represents the force density required to accelerate the source fluid from 

-J- P-1 

the source velocity v to the flow velocity v . The third term VP 

o p o 

: o 

a "bouyancy" force which acts upon density variations in the presence of the 

unperturbed pressure gradient. The equation states that the vorticity moves 

in such a manner that the Magnus force, source force, bouyant force, and 

electromagnetic force densities balance. 


is 


For application of this equation to an isolated vortex tube, it is con- 
venient to integrate it over the cross-sectional area of the tube. This gives 

-j- <VP > 

p < U - v > x r + <v - v > 0, 5— M. - I, x <B > = 0 

o olostLpllo 

o 

where is the circulation of the vortex tube, (1 is the mass pf source 

fluid being added to the vortex tube per unit length, M is the excess mass 
present in the tube per unit length, and 1^ is the total axial electric 
current in the tube. Brackets <> denote mean values over the cross section 
of the tube. The source term, the bouyancy term, and the electromagnetic 
term each make the motion of the vorfex tube deviate in some manner from the 
zero order fluid velocity v . The effect of the I x <B > term on vortex 
motion is well known, and our comments will be restricted ?o the other terms. 

One application of the source term is to the effect of the fluid boundary 
layer at a wall on the motion of a vortex whose axis intersects the wall. If 
fluid is sucked up from the boundary layer into the vortex as it moves (which 
will occur if the vortex is sufficiently strong to have a low pressure core) 
the boundary layer constitutes a fluid source with velocity v =0. The 
effect, as seen from solving the equation for this case, is to S cause the 
vortex to move to the right of the direction of the unperturbed flow v . 
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The bouyancy term applies to vortices which have a low density core. All 
vortices in a compressible fluid which have flow velocities approaching 
the speed of sound are naturally bouyant. Weaker vortices as well as 
vortices in an incompressible fluid may be made bouyant by injecting a 
lower density fluid into the core. The less dense fluid remains in the 
core of the vortex provided the radial pressure gradient in the vortex 
exceeds the other pressure gradients in the flow. One effect of bouyancy 
on a vortex whose axis lies in a horizontal plane is to cause the vortex 
to move laterally in the plane under the action of the downward pressure 
gradient due to gravity. This should cause the trailing vortices from the 
wing of an airplane* if made bouyant* to approach each other. Another effect 
of bouyancy is to cause two parallel vortices in near proximity to exert an 
attractive force on each other. This attractive force arises because each 
low density core is in the radial pressure gradient field of the other vortex. 
This, effect on vortices of like sense, which are revolving about each other, 
is to reduce their revolution rate. The effect on vortices of opposite 
sense and equal strength, which are drifting at constant separation, is to 
increase their drift velocity. The maximum changes in drift speed and 
revolution rate of vortex pairs due to bouyant attraction is on the order 
of 10%. Larger effects result from the downward pressure gradient caused 
by gravity. 
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LOW-ENERGY ION ATOM EXCITATION 
R. Novick 

Columbia Astrophysics Laboratory 
Columbia University 
New York, New York 10027 

ABSTRACT 

Experimental and theoretical studies have been made of the excitation 
cross sections of various neutral excited states of He in low-energy He**" - He 
collisions. Two striking results of the experimental observations are the 
large value of the cross section at low energies with thresholds slightly above 
the final state energies, and the strong oscillatory dependence of the cross 
sections- on projectile energy. These phenomena can be interpreted in terms 
of the molecular -potential energy curves of the He^ system. To understand 
the low-energy excitation it is noted that, in the diabatic approximation, 
tbe gerade ground term is repulsive and pseudocrosses all of the excited 
state terms. 2. These pseudocrossings occur at a value of internuclear separa- 
tion, Rj, of about 1 a.u. In the course of collisions with small enough 
impact parameters, the internuclear separation decreases and passes through 
Rj. As the system passes through the pseudocrossing, excitation occurs by the 
mechanism discussed by Landau, Zener, and Stuckelberg. 3 This accounts for the 
large cross section at low energies. Since the thresholds occur slightly above 
the final state energies, they provide a direct measure of the energy location 
of the pseudo crossing .4 Previously it was proposed that the oscillations were 
produced by an interference effect which arose from the double passage of the 
system through the crossing at Rj. However, it has now been shown 3 that any 
phase difference thus developed is strongly dependent on the impact parameter 
b; consequently, the integration over impact parameters required to obtain 
total cross sections averages out this interference effect. Such a phase 
difference might contribute to the oscillations reported in differential 
inelastic cross sections^ hut cannot explain the structure observed in the 
total cross-section data. 

In a recent paper, Rosenthal and Foley show that in addition to the 
crossing between the ground state and excited state curves at R-^, other 
crossings between tiro or more excited states exist at relatively large values 
of internuclear separation, R Q ~ 20a. u. In this model the excitation mechanism 
is essentially identical with than previously considered and hence the large 
cross sections at low energies, and the threshold energies, are explained as 
before. However, the oscillations are now considered to arise when two excited 
states populated at Rj develop a phase difference between Rj and R q before 
mixing again at R . Since the resulting phase difference is weakly "b "-dependent, 
it manifests itself in the form of oscillatory structure in the total cross 
section. Rosenthal and Foley show that their model correctly predicts the 
oscillatory structure in the excitation functions of the 3^S and 3o states, 
and, moreover, they predict that these oscillations should be in antiphase as 
is in fact observed. 


This work was supported by the National Aeronautics and Space 
Administration under Grant NGR 33-008-009. 
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New experimental studies have been made of the 2P excitation functions 
and of the velocity dependence of the cross sections. The latter studies were 
accomplished by using He^ and He^ in various combinations. These new data 
provide strong additional support for the theory of Rosenthal and Foley. 

The results presented here have stimulated considerable efforts to 
understand the detailed interactions occurring in low-energy ion-atom collisions. 
The theory recently developed by Rosenthal and Foley has been remarkably 
successful in leading to an understanding of those results. Some un an swered 
questions remain, however, and further theoretical efforts will be needed to 
answer them. Xt is clear that the full elucidation of the detailed structure 
of the excitation functions will provide a powerful test of our underst an ding 
of the dynamics of the unbound molecular state formed during these inelastic 
conditions. 
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Pig. 2. (a) Excitation cross section for the 4713-1 He x line in 
4+4 4+3 

He - He and He - He collisions plotted verses relative velocity 
at infinite internuclear separation, v(»). (B) Excitation cross 

section for the 4713-1 He I line in 4 He' f ‘ - 4 He and 4 He + - 3 He 
collisions plotted versus relative velocity in the final state. 
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PLASMA CHEMISTRY AND ION-MOLECULE INTERACTIONS 
by John V. Dugan, Jr., NASA Lewis Research Center 

Calculations of non-equilibrium number densities N e have been extended 
to argon plasmas for electron temperatures from 8,000 to 15,000 kJ 
The N e values are far below Saha at intermediate values of neutral 
density for optically thick plasmas. The results are non-hydrogenic 
for T e = 8,000 K; however, N e values are sensitive to the choice of 
model atom. 

Studies of an expanding hypersonic air stream have begun. The condi- 
tions consist o-f temperatures from 5^000 to 50,000 K and pressure 
ratios P/P 0 from 10"^ to 10 where P Q is atmospheric pressure. Pre- 
liminary calculations are for a nitrogen model with species N2> N, N2> 
N + , N ++ and e~. The relaxation times for electron energy equilibra- 
tion with heavy species are calculated from (elastic and inelastic) 
semi-class ica 1 collision coefficients. 

Numerical calculations of ion-molecule collision trajectories have been 
extended to including vibrational degrees of freedom on lifetimes of’ 
ion-molecule collision complexes. The nature and formation of such 
complexes has been studied via computer plots and computer-made motion 
pictures The nature of hindered dipole motion has been studied also. 
This behavior has chemical implications. 


REFERENCES 

1 J. R. Rose and J. V. Dugan, Jr., (unpublished). 

2 J. V. Dugan, Jr., R. B. Canright, Jr., R. W. Palmer and J. L. Magee, 
presented at VI international Conference on Electronic and Atomic 
Collisions, MIT Press, pp. 333 — 3^+0 . 



197 


TURBULENCE IN PLASMAS WITH COLLISIONS 
G. Sandri 

Aeronautical Research Associates of Princeton., Inc. 
50 Washington Road, Princeton, New Jersey 085^0 


1. Plasma turbulence and gasdynamic turbulence are analogous 
but differ In many respects. The major interaction among 
plasma particles is the electrostatic Coulomb force which is 
analogous to the intermolecular force in a neutral gas but 
has an infinite range. As a consequence, plasmas behave, 
under some conditions, like neutral gases but are much more 
unstable . 

2. It has been possible to develop a satisfactory theoretical 
description for dilute neutral gases which starts with the 
equations of motion of each of the particles in the presence 
of the forces exerted by all the others (Liouville equation) 
and leads through systematic expansions, first to a micro- 
scopic description (kinetic equation) and then to a macro- 
scopic description (fluid equations). Fluid-dynamic turbu- 
lence is described by macroscopic equations. A parallel 
theory is available for plasmas in weak magnetic fields if 
they are microstahle. This class of plasmas can be described 
by a well established kinetic equation which simplifies, for 
most practical instances of stable plasmas, (e = dielectric 
function ~ l) to the Landau-Pokker-Planck equation. On the 
basis of this equation, Spitzer and co-workers have calculated 
the transport properties of stable plasmas, in satisfactory 
agreement with the available observations. Furthermore, from 
the kinetic equation for stable plasmas, one can deduce the 
magnetohydrodynamic equations and, in principle, describe 

MED macro turbulence . 

3 . A 'Complication in the description of' plasmas arises from 
the fact that strong magnetic fields are often encountered 
in plasmas of interest. For "collisionless" plasmas subject 
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to strong magnetic fields, a theory has been developed (the 'Chew- 
Low-Goldberger theory) which is based on an expansion in powers of 
a parameter proportional to the Larmour radius. This theory shows 
a remarkable feature in contrast to the description that applies 
to neutral gases; namely, the kinetic and fluid descriptions can- 
not, be sharply separated. A further complication in describing 
plasmas is due to the fact that the condition of stability 
(e / 0) is seldom met in plasma work. For a collisionless system, 
linear microinstabilities have been investigated extensively and 
what is believed to be a fairly complete classification is now 
available (in fact, a beginning of nonlinear instability theory 
exists). The theory of MKD macroturbulenee has not been developed 
to a great extent, largely because of the belief that micro- 
instabilities are so damaging as to require the most immediate 
attention. A standard approach to the problem of microturbulenee 
is the "quasilinear" theory, which is designed for collisionless 
plasmas. We outline, below, a program for investigating micro- 
turbulenee in plasmas when collisions cannot be neglected. 

4. The first part of our program consists in investigating the 
behavior of turbulent magnetized plasmas by a method which is 
quite different from the quasilinear theory. Our approach 
consists in determining first the velocity distribution of 
charged particles in the presence of random electric and magnetic 
fields, and making then the velocity distribution function con- 
sistent with Maxwell's equations. To present, our work has been 
mainly concerned with the first phase of this approach. We find 
that it is possible to derive a kinetic description for particles 
subject to random magnetic fields from the Liouville equation for 
the entire system. The kinetic equation is satisfied by the 
ensemble average of the velocity distribution function of the 
particles. This equation can be shown to reduce to a Fokker- 
Planck equation if the magnetic field fluctuations are not very 
strong and the kinetic energy of the particles is large. This 
implies that the Fokker-Planck equation for magnetic Brownian 
motion is limited to situations that correspond to 
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average fields which are not very strong. When the mean field is 
very strong, or, alternatively, the kinetic energy of the particles 
is small, a considerably more complex description of the system 
than that afforded by the Fokker-Flanck equation is needed. We 
have not completed a derivation of a suitable equation in this 
case. We find, however, that it is possible to describe a number 
of situations in terms of equations representing both the pre- 
cession of the 'velocity vectors about the mean field and the 
''decay 1 ' of the distribution due to the interaction with the turb- 
ulent magnetic field. The precession frequency and the decay 
coefficient are 'somewhat complicated functions of the autocorre- 
lation of the magnetic field fluctuation. In general, the decay 
coefficient is negative, thus making the system unstable. 
Furthermore, the presence of turbulence makes the system of 
charged particles exhibit paramagnetic or diamagnetic behavior, 
depending on the structure of the autocorrelation function. The 
effect is exhibited by the shift in the precession frequency 
induced by the turbulence. 

5. The second part; of our program aims at developing a kinetic 
theory for unstable, spatially homogeneous plasmas. For 
example, a plasma with the "loss cone" Instability is envisaged. 

The techniques that we use stem from our analysis of stable 
plasmas. In this case, the two-particle velocity distribution 
function, as obtained by standard methods (adiabatic theory), 
has an infinite spatial range when the relative velocity vector 
is nearly parallel to the relative separation vector. We can 
show, however, that the two-particle correlation function has 
a finite spatial range if the effect of the decay of the velo- 
city distribution is properly taken into account. The method 
that has allowed us to develop a satisfactory two-particle corre- 
lation for stable plasmas is a local study of the phase space in 
the region where the adiabatic theory yields an infinite range. 
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The loss ccne Instability is analyzed with the same method., i.e.j 
in the vicinity of e = 0 where the two- particle correlation 
function diverges. The resulting two-particle distribution 
function will then be utilized to determine the evolution in 
time of the velocity distribution function which is initially 
unstable . 

6. We discuss the results available on the unstable homogeneous 
plasmas, as well as the work in progress on understanding the 
behavior of plasmas in the presence of fluctuating magnetic 
fields. We also outline methods for extending our approach 
to electrostatic fluctuations and for making the velocity 
distribution function consistent with Maxwell's equations. 
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HIGHER APPROXIMATIONS TO PLASMA TRANSPORT PROPERTIES 

W- E. Meador 

NASA, Langley Research Center 
Hampton, Va. 


The purpose of this study is to investigate the accuracies of the 
standard Chapman-Enskog and Grad 13-moment approximations in predicting 
transport phenomena in plasmas under a variety of conditions. Both methods 
are questionable when first-order perturbation limits are exceeded - a.g., 
when electron diffusion velocities are sufficiently large to invalidate 
linear flux-force relations. In addition, since the Grad approximation is 
specifically an expansion in the eigenfunctions of Boltzmann’s binary elastic 
collision operator for Maxwellian molecules, the reliability of that method 
must' be checked for more extreme potentials even to first order. Hence an 
investigation of the limitations of these techniques is imperative for a 
critical evaluation of much of the current literature. 

The method employed in the present research for determining these limi- 
tations involves the use of an approximate kinetic equation which is solved 
exactly and to which the Chapman-Enskog and Grad approximations can also be 
applied. Comparisons between the calculated results for various transport 
and other plasma properties thus provide valuable insight into the desired 
accuracies. The essential features are summarized in the following paragraphs. 

First-order (i.e. , nearly equilibrium) plasmas are found to be described 
quite accurately by the corresponding Chapman-Enskog theory, but not as reliably 
in some aspects by the 13-moment approximation. While the latter method predicts 
the diagonal transport coefficients of electrical and thermal conductivity to 
within one or two percent, errors on the order of 30 percent exist for such 
off-diagonal coefficients as thermal diffusion when the interaction potentials 
are very soft or very hard. Similar errors occur in the entropy density, the 
collisional production rate of entropy, and second-order contributions to the 
electron pressure tensor, all of which are corrected to nearly the exact re- 
sults by the addition of a few higher moments. Third-order contributions to 
the electron current density and the phase angles by which the current density 
and the heat flux lag an applied oscillating electric field are poorly predicted 
by the 13-moment approximation, the errors in these cases corresponding to 
factors up to four in the limits of hard and soft particles. 

All of the aforementioned errors ate explained by the inadequate descrip- 
tion of certain collision integrals when the 13-moment approximation is used 
to close out the macroscopic equations of change. However, a new phenomenon 
appears when a steady-state plasma becomes spatially inhomogeneous through 
second order by means of an applied electric field inducing transverse pressure 
and temperature gradients. The prediction of such gradients by the 13-moment 
approximations is extremely poor (even for Maxwellian molecules) because of 
the neglect of a new and important driving force for heat transfer which 
appears with the addition of the next nine moments. A similar phenomenon 
occurs in the viscous flow of neutral gases through circular pipes and pre- 
sents a somewhat different description of heat transfer from the picture of 
Navier and Stokes and the perturbation expansion of Chapman and Enskog. 
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The following publications have been generated by this research: 

(1) A Semiempirical Collision Model for Plasmas. NASA TR R-310, 1969. 

\ 

(2) A Critical Analysis of the Grad Approximation for Closing Out the 
Magnetohydro dynamic Equations for Plasmas. NASA TR R-325, 1969. 

(3) Third-Order Contributions to Electrical Conduction in Plasmas. 

NASA TN D-5477, 1969. 

(4) Induced Pressure and Temperature Gradients in a Plasma as a Result 
of Weak Electric Fields. NASA TN L-6877, 1969. 

(5) Fourth Approximation to Heat Transfer in Viscous Gases Flowing 
Through Circular Pipes, with Constant Temperature Gradients. NASA TR L-7018, 
1970. 
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NONLINEAR TRANSITION LAYER SOLUTIONS FOR PLASMAS 
USING GENERALIZED FUNCTIONS 

Robert C. Costen 
NASA, Langley Research Center 
Hampton, Va. 


Progress has been made in the application of generalized functions to 
obtain transition layer solutions of the Vlasov equations (Maxwell T s 
equations and the collisionless Boltzmann equation with electromagnetic 
forces) . The objective is to find analytic solutions for the fine structure 
of layers of rapid transition in a two— component plasma with an applied 
magnetic field for possible application to the magnetopause. It is hoped 
that this may also lead to solutions for the structure of collisionless shocks. 

Layers of rapid transition are well suited for mathematical description 
by^the use of generalized functions. For example, consider the representation 
F(r,w,t) H^(x) for a distribution function, where F is an analytic function 
of position r, velocity w, and time t, and H (x) is the Heaviside unit 
step function defined, for example, by the sequence 

H r (x) = | (1 + tanh nx) 

where n is the sequence index. This representation F(r,w,t)H (x) 
engenders a rapid variation in the distribution function near the plane x = 0 
(plasma-vacuum interface) and slow variations elsewhere. The sequence index 
n governs the thickness of the transition layer; as n -*• oo the thickness of 
the transition layer goes to zero. 

A formal method has been developed for finding solutions of the Vlasov 
equations where the distribution functions and electromagnetic fields are 
all generalized functions. The hyperbolic tangent representation of the step 
function (given above) is adhered to throughout because its deri vatives have 
the convenient form 


S> n (x) - 2n (H n (x) - H n 2 (x)) 

S '(x) = 4n 2 (H (x) - 3H 2 (x) + 2H 3 (x)) 
n n n • n 

where £ (x) is the Dirac delta function sequence. The acceleration term in 
the Boltzmann equation 

( jr + ) F + — (E + w X B) • — F « 0 

is nonlinear. Therefore the generalized function formalism used must embrace 
products. This is a departure from the usual generalized function theory, 
where products are not defined except in a convolution sense. The extension 
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to products is a simple one for the family of generalized functions treated 
above, for any product takes the form of a power of n multipolied by some 
polynomial in H (x) . (The consideration of whether such products are 
integrable in the limit n oo is not relevant in the physics of transition 
layers, because n, which governs the width of the transition layer, always 
remains finite.) Thus a class of generalized functions is defined which 
includes products of the members of the class. The basic element of the class 
is n^H^Cx), where the power is an integer and the power /3 is a 

non-negative integer. 


We seek solutions of the Vlasov equations from this class of generalized 
functions; that is, each field in the Vlasov equations is taken to have either 
the form 


B.ff.t) - ZL B <?>t) n*H^(x) 


n 


or the form 


ci = -A 0 = 0 
A oo 


* n (r,w,t) =2_- ZL *d /3 (? Zt) ^H^(x) 

o< = "A 0 =0 


where the coefficients and F^ are analytic functions. The infinite 

summation over fi is assumed to be uniformly convergent and allows the shape 
of the field in the transition layer to vary as any analytic function of 
H (x). Negative powers of n are included because they can become signifi- 
cant when multiplied by positive powers of n in the nonlinear term of the 
Boltzmann equation. By substituting series of this type for each of the 
fields in the Vlasov equations and applying an identity theorem for general- 
ized functions, one obtains a set of equations for the series coefficients. 
Any analytic solution of this derived set of equations constitutes a 
transition layer solution of the original Vlasov equations. 


A method which has proven to be successful in obtaining solutions of the 
derived set of equations is as follows: Pick a priori the highest allowable 

order (power of n) for each of the fields. The equation corresponding 
thereto, which is the highest order equation, will normally have the applied 
magnetic field (which is known) as one of the factors in the acceleration term 
of the Boltzmann equation. Thus the highest order equation is linear and can 
easily be solved for simple cases. One then proceeds to the equation corres- 
ponding to the next lower power of n. In the nonlinear term of this equation 
one of the factors is already known from solution of the previous higher 
order equation; hence this equation is linear also and amenable to solution 
in simple cases. By continuing to work downward in this fashion through 
equations of_decreasing order one obtains asymtotic solutions in the expansion 
parameter n ^ for each of the fields. The accuracy of these solutions 
increases as the transition layer becomes thinner (n increases) . 
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RESEARCH OH THE LOW DENSITY HOLLOW CATHODE 
Guiseppe Palumbo, Paul J. Wilbur, and William R. Miehelsen 
Colorado State University, Ft. Collins, Colo. 


A mercury vapor hollow cathode discharge is being investigated 
for application to electron bombardment thrusters. The hollow 
cathode is located at one end on the axis of an annular pole piece 
which carries magnetic flux to the main discharge of the thrusters. 
In conventional designs there is no magnetic field inside the pole 
piece. Preliminary experiments reported here indicate that substan- 
tial reduction in discharge power consumption can be attained by 
addition of a magnetic field internal to the cathode pole piece. 
Possible reasons for these results are discussed. 
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CONSTRICTED ARC PLASMA JETS 
Charles E. Shepard 

MSA Ames Research Center, Moffett Field, California 


Existing constricted-arc supersonic plasma Jets produce streams 
of very high enthalpy at relatively low densijy^in a variety of 
gases . Performance data for air and nitrogen 5 show that enthalpies 
corresponding to superorbital velocity can he produced at densities 
corresponding to flight at altitudes down to about 60 km. Although 
existing constricted-arc Jet facilities have proved very useful at 
Ames, there is a clear need for higher operating pressure in order 
to produce higher Reynolds numbers and to simulate flight at lower 
altitudes 

When the pressure is increased beyond the normal operating 
level, two related problems become apparent. The constrictor disks 
that make up the arc tube can fail either due to inadequate water 
cooling or to electrical breakdown (arc-over) between adjacent 
constrictor disks. For the 'Same total constrictor length, thinner 
disks reduce arc-over tendencies but are more difficult to cool . 

The performance of thinner, better cooled constrictor disks similar 
to those developed by Maecker has been studied both by means of 
computer calculations and by experiment. The results of the studies 
indicate that the maximum constrictor wall heat transfer rate can be 
increased from the present value of 3 kW/cm 2 to about 10 kW/cm 2 . 

Although water cooling is necessary for steady-state operation, 
heat capacity cooling offers some important advantages for short time 
operation. Not only are the parts cheaper and easier to make and 
modify, but also higher heat transfer rates can he maintained for 
pulsed operation than for steady-state water cooling, as shown in 
Figure 1. Two schemes for short time runs are being investigated. 
Millisecond operation from a LC lumped parameter discharge line has 
already been shown to be feasible in a contract study by General 
Dynamic s/Convair for Ames. Power levels of hundreds of megawatts for 
several milliseconds appear practical with about 10 6 Joule energy 
storage. Alternately, lower resistance inductors now being sought 
should permit test times of tens of milliseconds at many tens of 
megawatts power. 

Short time operation of a heat capacity cooled arc Jet from a 
conventional DC power supply may also be advantageous. Economy in 
construction and power supply costs for heat capacity cooling were 
important factors for a new 6 cm internal diameter constricted-arc 
Jet. This device is a prototype of a steady-state constricted arc 
Jet that will provide high radiative heat transfer rates to blunt 
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test bodies. The device consists of a 6 cm I.D. x 2.4 meter long 
constrictor tube followed by a divergent nozzle having an exit 
diameter of 10 cm. Thirty-six \ inch diameter tungsten rods form 
a heat capacity cooled multiple cathode at the upstream end of the 
constrictor tube. Ninety % inch diameter rods make up the multiple 
anode that is located along the nozzle surface. The design goals 
are a total enthalpy of ll6 MJ/kg and an impact pressure of 1.0 atm. 

Preliminary tests of the 6 cm diameter, heat capacity cooled 
constricted-arc jet at power levels up to l4 megawatts have produced 
total enthalpies in excess of 120 megajoules/kg and a test body 
impact pressure of about 0.4 atm. These tests indicate that useful 
heat transfer data can be obtained in the .2 to .5 second test 
intervals. As shown in Figure 2, calorimeter rake measurements 
indicate a nearly uniform radial heat transfer profile over about 
half of the exit diameter. 



Maximum Constrictor Wall Heat Transfer Rate, kW/cnr 


Limit based on 

melting at inner surface. 

T. = 1365°K 


Approximate limit for 
.water-cooling . 


Limit based on 
gas -seal temperature 


Test time, seconds 

Figure 1. Maximum constrictor heat transfer rate at inner wall as a function 
”■ of test time. 
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Figure 2. Approximate radial distribution of heat transfer rate based on time 
melt 1/8 inch diameter copper calorimeter rods. 


DEVELOPMENT PROGRESS ON A CROSSED 
MELD ROTARY- ARC IE A COEMEED VORTEX 
CONFIGURATION TOR ACCELERATIEG HXPERVELOCITY 
NOZZLE PLOWS 

Roger B. Stewart, NASA Langley Research Center 


During the past several years we have carried out experimental and theoretical 
studies of the flow "behavior in an arc heater employing an electromagnetieally 
driven confined vortex. These efforts have "been directed toward the develop- 
ment of a device that will produce high Reynolds number nozzle flows with 
gas velocities that approach reentry flight speeds. Several small (5 00 kw) 
arc driven devices have been operated extensively to determine the variation 
of performance and the ability to predict performance for different operating 
conditions and heater configurations. Early experiments are reported in 
reference 1 and a more recent range of performance is shown in figure 1 along 
with a schematic diagram of the device. At the high mass flow rate (10 gms/see) 
nearly 60 percent of the arc input power appears as net energy addition to 
the gas (thermal plus directed kinetic energy) . At the low mass flow rates 
the efficiency is 48 percent. Test results indicate that the overall 
performance rises significantly as the Hall parameter increases. This implies 
that substantial gains in performance can be obtained by operating at high 
magnetic field strengths. This is a unique advantage of these devices. 

The theoretical studies have been reported in reference 2. Closed form 
solutions for the azimuthal velocity, induced magnetic field and pressure 
distribution within the vortex have been obtained from the coupled set of 
electromagnetic and fluid dynamic equations within the MHD approximation. 
Numerical solutions to the non-linear energy equation have been obtained 
using a fifth degree polynomial for the variation of electrical conductivity 
with temperatureo Solutions to the energy equation provide the temperature 
distribution, Mach number distribution, and power required to heat and 
accelerate the vortex. Experimental measurements of pressure, velocity 
distributions inside the vortex chamber, and the power required confirmed 
the theoretical predictions. The experimental and theoretical studies 
have been applied to design and build a larger device with twice the arc 
gap and magnetic field strength, (B^ - k webers/meter 2 ). This device 

will soon be tested. In addition, the Structures Research Division, at 
Langley Research Center, is constructing a third and larger device and 
the Australian National University at Canberra is also contemplating the 
construction of s imi lar type of arc heater with an accelerating vortex. 

The ultimate performance of these devices is net known at present, however, 
experimental studies with these new devices will define new performance 
levels and furnish scaling data. Figure 2 lists typical performance. . 

In addition to the vortex studies we have carried out experiments with a 
device that incorporates a mechanically and electrically coupled linear 
accelerator section downstream of the circular channel, figure 3* Higher 
pressures and larger plasma interaction lengths are an advantage of this 
two-stage device over a single accelerating cycle. Steady state operation 
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of the vortex-linear device has shown a 5 8 percent efficiency for energy 
input to the gas <2 Considering the relatively low magnetic field strength 
(l. 1 wehers/meter ) this represents a considerable increase in performance 
over the single stage vortex device operating at the same field strength. 

Basic research studies applicable to both accelerators concerns the funda- 
mental behavior of high current arcs that are driven over cold metal surfaces 
by electromagnetic means. We have carried out an extensive series of 
experiments on arc behavior which are reported in references 3 a n d h. Of 
particular interest is the large increase in current sheet (or diffused arc) 
rotational velocities in the spiral cathode devices compared with concentric 
circular electrode geometries. An attempt to link this behavior with the 
radial displacement effect of the highly ionized region of the arc at the 
cathode has been made in reference 5* An order of magnitude argument 
indicates that significant additional "local" Lorentz forces may be 
present on the cathode sheath in spiral vortex geometries as compared 
with circular geometries. Figure k indicates the basic geometries studied 
and shews the results with velocity as a function of driving force for 
these different configurations. The cathode current sheet velocities 
for the spiral cathode are more than four times higher than the velocities 
for the circular cathode* The velocities in excess of 8000 meters/second 
are believed to be higher than any yet reported at the pressure levels 
shown. 


The behavior of several particular rotary arc devices has been studied 
both theoretically and exp er ime nt ally. The confined vortex flow that is 
heated and driven electromagnetically appears to offer higher simultaneous 
velocity and Reynolds number performance than is presently available 
with conventional arc heaters. Upcoming experiments with larger devices 
will define the hypervelocity nozzle flow conditions that can be achieved. 
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Conditions at 2.0 "by 3.8 cm 
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Figure 2. - Experimental conditions achieved in 500 kw vortex 

heater accelerator. 
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1 channel with a. nozzle .blow-down system. 


Stream static pressure, psia • . 

|Tot£L- enthalpy, Btu/lb . .• . <. . . ; ,>-» 
Weight flow, Ib/sec . '. . .. 

'Efficiency by heat balance,' percent . . 
i Applied 'magnetic field," -weber /met er^ 
l Arc -.current,-, amperes ..V .; , V ... h 
lArcV volt agey r volt s . Y. .• " Y Y-VV 

f* .“*•*. * * i 1 * { ■(!• ’* ' N 




3.1 

li,4oo 

*>.017 

‘58;o' ! 

*.14 


•-8C5 

•'44q 


jJl'gureVJj.V-" Vortex;. linear ; plasma’ device,’' 






215 


B aux.' weber meter2 

^ Cathode "A" 1.5 cm displacement 2.0 
^ Cathode "A" 1.25 cm displacement 1.0 
© Cathode n B", no displacement 65 and i 3 

& Cathode "A^", no displacement i.o 


Solid symbols denote tests with end wails 
Open symbols denote tests without end walls 



Electromagnetic driving force, BI, newtons/ meter 

Figure k. - Spiral and circular cathode experimental results* 
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COMPARISON OF THEORETICAL AND FLIGHT-MEASURED IONIZATION 
IN A BLUNT BODY REENTRY FLOW FIELD 

P. W. Huber, J. S. Evans, and C. J. Schexnayder 
NASA, Langley Research Center 
Hampton, Va. 


In order to define, circumvent, or alleviate the radio attenuation 
which occurs during spacecraft atmospheric entry one must first be able 
to describe the magnitude and distribution of electron concentration in 
the plasma flow near the antenna, and then to determine the interaction 
of the em field with this plasma. In the real-life situation, many factors 
such as reaction rates, ablation, separated flow., angle of attach, antenna 
configuration, etc. , combine to make the complete problem too complex to 
assess on a purely theoretical basis. As an alternative, one lays the 
groundwork using idealized theoretical models for various facets of the 
problem and then assesses the applicability of these models by comparison 
with data from the laboratory, or with flight measurements when the problem 
cannot be properly simulated on the ground. One purpose of the Langley 
Research Center project RAM program is to obtain reentry flight data which 
can be interpreted in terms of models for the plasma flow, for em propagation, 
and for alleviation of signal attenuation. Towards this end, RAM vehicles 
of hemisphere-9 0 cone configuration have reentered the earth’s atmosphere, 
at near-orbital velocity, at near-zero attitude, and with both ablating and 
nonablating heat shields. They were equipped with various onboard plasma 
diagnostic sensors, various communication systems, and a material injec- 
tion scheme to modify the plasma and reduce attenuation. ■ The vehicle’s 
nose radius is 0.5 feet. 

The purpose of this paper is to draw conclusions relative to the re- 
entry plasma ionization by comparison of results from theoretical models 
with the onboard diagnostic data and with the passive data obtained from 
interpretation of the signal attenuation. 

The theoretical models which are employed in the comparison include 
the following: (1) Pure air inviscid and merged boundary-layer- inviscid 

blunt-body flow-field models in which the various air species chemical 
kinetics reaction rates can be treated as independent parameters, (Ref. 1); 

(2) pure air boundary-layer correlations in which the influence of ambi- 
polar electron-ion diffusion and wall recombination are exhibited; 

(3) boundary-layer program in which the effects of alkali ablation contami- 
nants on the ionization profiles are exhibited, and (4) a nonhomo geneo us 
slab plane wave em plasma propagation model for determining signal attenua- 
tion due to the plasma (Ref. 2). 

The flight data used in the comparison includes onboard measurements 
of collected current from a rake of eight Langmuir probes extending 7 cm 
into the flow field at the aft part of the vehicle, onboard reflected power 
measurements from 15 reflectomers located at four station along the body 
(four different frequencies) , and onboard measurements of VSWR for three 
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of the communications antennas. Ground-based received signal strength 
measurements for different polarizations were made for the four communi- 
cations' antennas (four different frequencies) from a number of different 
range stations in radio contact with the spin stabilized reentry vehicle. 


The vehicle reentered at 25,100 ft/sec and its attitude was maintained 
at less than three degrees during the reentry data period. RAH C-l had 
a NARMCO ablating heat shield which contained more than 10^ ppm of sodium 
and potassium, while RAM C-2 had a beryllium and teflon heat shield which 


was free of alkali contaminants. 


Some of the conclusions which have resulted from the comparison of the 
theoretical with the experimental data are briefly summarized below , along 
with the basis upon which they are reached: 

-f. 

1. NO air reaction rate .- Computer experiments show that the peak 
profile ionization, and the C and X— band signal attenuation for pure air 
are quite sensitive to the N0 + rate in the altitude range above about 
120 Kft. C-band signal attenuation data in this altitude range for the 
flight with a beryllium heat shield followed by a teflon body coating shows 
good agreement with that computed using the NO^ rate recently measured in a 
shock tunnel, reference 3, while the X-band data agreement is not as good, 
but closer to this rate than the older rate of reference 4 which is slower 
by a factor of three. These results are shown in Fig. 1. 

2. Ablation ionization .- Theoretical ionization profiles and signal 
attenuations were computed fox cases where alkali material is contained in 
the boundary layer due to the NARMCO heat-shield ablation. Comparison of 
the signal attenuation data for the NARMCO flight, RAM C-l, with that of 
the beryllium-teflon flight, RAM C-2, and comparison of the theoretical 
model for pure air with the ablation impurity model both show a large effect 
of the ablation impurities. Furthermore, the air model using the fast NO 
rate is in reasonably good agreement with the RAM C-2 attenuation data, while 
the ablation model, although not agreeing quantitatively, shows the same 
qualitative trends as the RAM C-l data below altitudes of about 180 Kft. 
Comparison with the low altitude signal recovery data, where ablation rates 
are highest, is shown in figure 2. In figure 3, attenuation data from 

RAM C-l at higher altitudes is compared with the ablation model to suggest 
the very interesting probability that the impurity ionization disappears 
at about 200 Kft due to lack of ionization time. This is seen in the sudden 
drop of attenuation as compared with local equilibrium theory. 

3. Electron-ion diffusion .- Theoretical correlations are made for the 
RAM high altitude regime which show that the effects on peak profile ioniza- 
tion due to ambipolar electron-ion diffusion and wall recombination are very 
large at altitudes much above 230 Kft. Reflectometer measurements and VSWR 
indications of the peak electron concentration confirm this result, while the 
Langmuir probe profile measurements show, in addition to this large reduction, 
a pronounced distortion of the profile away from the wall as compared with a 
theoretical model with no electron-ion diffusion. 
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COMPARISON OF HIGH ALTITUDE ELECTRON CONCENTRATION 



FIGURE 4. 
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THEORETICAL ionization profiles in blunt body 
ATMOSPHERIC ENTRY FLOW FIELDS 

J. S. Evans, C. J. Schexnayder, and P. W. Huber 
NASA, Langley Research Center 
Hampton, Va. 


Techniques are described for the computation of electron concentra- 
tion profiles in the shock layer surrounding a blunt-nosed- vehicle during 
entry into Earth’s atmosphere at near-orbital speed and -at altitudes where 
ion diffusion effects can be neglected. Emphasis is on careful treatment 
of finite rate chemistry and boundary-layer effects. A streamtube approach 
is used to calculate gas composition and flow properties in both the invis- 
cid and viscous parts of the shock layer. Each streamtube (see fig. 1) 
begins immediately behind the shock front with free-stream composition and 
fully excited internal energy. Simultaneous integration of flow equations 
and chemical kinetic equations is carried out for streamlines defined by 
specified pressure distributions. The equations and their solution are 
modeled after reference 1. For the body shown (sphere-cone; = 1 ft., 
cone half angle = 9°) the boundary layer is too thin to have an appreciable 
effect on N e profiles below about 175,000 feet and nearly fills the shock 
layer at 230,000 feet. Thus, increasing numbers of streamlines are 
swallowed by the boundary layer as altitude increases. Streamlines are 
followed into the boundary layer on the basis of a mass flow criterion, 
but only gas composition is computed, since the flow properties are found 
from houndary-layer calculations. 

The species and the reactions are shown in figure 2. The atmosphere 
is assumed to consist of 0^ (21%) and ^ (79%). All possible atoms, 
diatomic molecules, and their singly charged positive ions are included, 
and, in addition, the negative ions 0“ and O^ - are also. The reactions 
are divided into three groups: neutral, charge production, and charge ex- * 

change. Since the charged species make up only a tiny fraction of the total 
gas, the thermodynamic properties of the mixture is for all practical pur- 
poses controlled by the neutral reactions. Creation and destruction of 
charges takes place only through the charge production reactions. The charge 
exchange group is important because it provides paths for rapid buildup 
and/or destruction of species which wou^_d otherwise be slow to change. For 
example, the rate of recombination of 0 by reaction 11 (3-body recombination) 
is negligible at 175,000 feet but the rate of exchange with NO to produce 
N0 + is rapid. Then N0 + combines with e . 

The importance of adequate provision for charge transfer is illustrated 
in figure 3, where, with no charge transfer reactions, O'*" and FT*" build 
up to- smaller initial values than they should and then decrease not at all 
during the expansion. A limited system of two charge transfer reactions, 
which connect O**' to 02 "*" and If*" to Ng"*", gives higher initial values 
but they are still essentially frozen during the expansion. In the latter 
case e~ remains quite high, because of the failure of the atomic ions to 
disappear. The effect of calculations like the above on N e profiles at 
1 5 6,000 ft, x/djj = 4, is shown in figure 4. 
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In figure 5 the body streamline values of y are plotted against 
altitude at several body stations and are compared with local ionic equi- 
librium values. (For local ionic equilibrium all ion concentrations are 
in equilibrium with local concentrations of the neutral species at the 
local temperature.) Since local ionic equilibrium represents a lower bound* 
electrons cannot recombine as completely at 100*000 feet as they can at 
175,000 feet. If the lower bound can be removed then further recombination 
can take place. .One way to lower the bound is by injecting water into the 
flow as shown on/?igure. It was calculated at altitude = 131,000 feet, 

W* = 1 (W* = mass flow of water/mass flow of air) by means of another stream- 
tube program (ref. 2) which has been developed to study the effects of water 
injection on N e profiles during atmospheric entry. This analysis accounts 
for the presence of accelerating and evaporating water droplets in the flow, 
as well as chemical reactions introduced by adding hydrogen. 
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FLOW FIELD AND TYPICAL STREAMLINE PATHS 



FIGURE 1. 


SPECIES AND REACTIONS FOR EARTH ENTRY 


NEUTRAN : 

1. 0g+M-*0+0+M 

2. Ng + M - N + N + M 

3. NO+M-N+O+M 

4. O + Ng - NO + N 

5. O + NO •* Og + N 


FIGURE 2. 


SPECIES 

02, Ng, O, N, NO, NO^j Og+, N 2 + , 0+, N*, e“. 


REACTIONS 
CHARGE PRODUCTION: 

6. O + O -* Og + + e” 17. 

7. N + N -* Ng + + e“ 18. 

8. N+ O'* NO + + e~ 19. 

9. Ng + e- -» Ng + + e“ + e“ 20. 

10. N + e" -* IT" + e~ + e“ 21. 

11. 0+ e~ -*0 + + e“ + e“ 22. 

12. Og + 6“ Og + 6“ + 6"” 23. 

13. NO + e~ - NO + + e- + e~ 24. 

14. O + O -» 0 + + O* 25. 

15. 0+ N- N* + O' 26. 

16. *0 + NO -* NO + + O" 27. 


28. 

29. 

30. 

31. 

32. 

33. 

34. 

35. 

36. 

37. 


O', Og“ 

CHARGE EXCHANGE: 
Ng + NO+ - Ng + + NO 
N + NO + - N + + NO 
O + NO + - 0+ + NO 
02 + N0 + - Og + + NO 
Ng + Og + “• Ng + + Og 
N + 02^ - N 4 - + 02 
0 + Og + - 0 + + 02 
Ng + 0 + - Ng + + O 
N +• 0 + -* H 1 * +‘0 
Ng + N+ - Ng + + N 
0' + M-*0 + e' + M 
Og" + M -* 02 +• e” + M 
Og 4- e~ -• 0~ + O 
O + N0 + - N+ + Og 
NO + N0 + - N2 + + 0 2 
O + N0 + - 02 + + N 
N + N0 + - Ng + + O 
N + N0 + - 0 + + N2 
NO + NO + -» Og + + Ng 
N2 + 0 + - N + + NO 
NO + 0 + - N+ + 02 


38. N + 02 + -* 0 + + NO 

39. O + N2 + - N+ + NO 

40. 02 + 0~ -» 02“ + O 
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STREAMLINE VARIATION OF CHARGED SPECIES 
(ALTITUDE = 156,000 ft) 



NO CHARGE TRANSFERS 
■TWO CHARGE TRANSFERS 



FIGURE 3. 


EFFECT OF CHARGE TRANSFER REACTIONS ON N e PROFILES 
(ALTITUDE = 156,000 ft; x/% = 4 . ) 





IONIZATION LEVEL AS A FUNCTION OF ALTITUDE FOR BODY STREAMLINE 



FIGURE 5. 


ALTITUDE, ft 



